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Preface 


This annual report contains information describing highlighted accomplishments for the 
past year from NASA's Aerodynamics Research and Technology Program. 

The Aerodynamics Research and Technology Program is conducted as part of the response 
to NASA’s charter to: (1) Preserve the role of the United States as a leader in aeronauti- 
cal science and technology, and the application thereof; (2) Improve the usefulness, perfor- 
mance, speed, safety, and efficiency of aeronautical vehicles; (3) Supervise and direct the 
scientific study of the problems of flight with a view to their practical solution; and (4) 
Ensure the timely provision of a proven technology base for a safe, efficient, and environ- 
mentally compatible air transportation system. 

The Aerodynamic Research and Technology Program includes both fundamental and 
applied research directed at the full spectrum of aerospace vehicles, from rotorcraft to 
planetary entry probes. The program encompasses analytical, computational, and experi- 
mental efforts conducted using the world’s best wind tunnel, computational, and flight 
research facilities at NASA’s Ames and Langley Research Centers. 

It is impossible to do justice to the breadth and depth of this comprehensive research 
program in a single document. The intent of this report is to render a balanced view of 
NASA’s Aerodynamics Program by presenting accomplishment highlights from fiscal year 
1991. Other documents and conferences provide a more extensive forum for detailed 
coverage of a given area. For brevity, we have not included all of the highlights that were 
submitted for consideration, and we apologize to those whose submissions were not 
incorporated. 

This report is arranged in chapters which outline the applied research and technology 
programs by vehicle type and the fundamental research by subject area. Each chapter 
includes an introduction by the appropriate Aerodynamics Division program manager, 
and each accomplishment highlight identifies the responsible researcher at the field 
centers. To facilitate communication, the address and phone numbers are also listed. 

During this year, the Aerodynamics Research and Technology Program was managed by 
the Aerodynamics Division, one of five divisions in the Office of Aeronautics and Space 
Technology (OAST) at NASA Headquarters. At the close of 1991 OAST was reorganized, 
and the five discipline divisions were replaced by separate aeronautics and space divi- 
sions in specific thrust areas. In line with this new organization, this will be the last 
report of this type covering the entire Aerodynamics Research and Technology Program. 

Mr. Louis J. Williams Dr. Kristin A. Hessenius 

Director Deputy Director 

Aerodynamics Division, Code RF 

National Aeronautics and Space Administration 

Washington, DC 20546 
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Chapter 1 

Subsonic Aerodynamics 

The objective of the Subsonic Aerodynamics Program is to provide the technology to improve the 
performance, efficiency, and economics of future transport aircraft. To accomplish this objective, 
the program is currently focusing on the areas of hybrid laminar flow control technology and high 
lift. 

A hybrid laminar flow control system has been installed on the wings of a Boeing 757 transport 
aircraft to investigate the effectiveness, maintainability, and cost of applying this technology to 
subsonic transports. 

High-lift studies include both wind tunnel tests and computational fluid dynamics (CFD) 
calculations. Two-dimensional wind tunnel tests that focused on three element airfoils have been 
conducted. Three-dimensional tests at high Reynolds Number and large scale are planned. CFD 
calculations are guiding the activity for both the two- and three-dimensional tests. 

Program Manager: Victor R. Corsiglia 

OAST/RF 

Washington, D.C. 20546 
(202) 453-2261 
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Figure 1-1. Hybrid Laminar Flow Control Flight Research 


2 



1-1 Hybrid Laminar Flow Control Flight Research 


Objective. To conduct flight research ex- 
periments with the Hybrid Laminar Flow 
Control (HLFC) system installed in the Boeing 
757 flight research aircraft. 

Approach. The program is jointly funded by 
the NASA, the U.S. Air Force, and Boeing 
Commercial Airplane Co. An HLFC system 
provides suction boundary layer control in the 
leading edge (to the front spar) to control the 
highly three-dimensional laminar boundary 
layer. Downstream of the front spar, the 
surface pressure distribution is used to stabi- 
lize the boundary layer flow and maintain 
laminar flow to the wing shock or aft pressure 
rise. The suction surface in the leading edge 
is a microperforated titanium skin with over 
19 million tiny, laser-drilled, closely spaced 
holes. 

Accomplishments. Last year the flight test 
experiments yielded momentous results that 
could advance significantly the application of 
this technology into the commercial transport 
fleet. The research showed that the suction 
system flow capacity, power requirements, 
weight, and complexity had been overesti- 
mated by more than 200%. 


Significance. The experiments of the past 
year have demonstrated that the suction re- 
quired to obtain large amounts of laminar 
flow over the wing chord (as much as 65% 
chord) is much less than had been previously 
believed, and amounts to one-third of that 
indicated by the early design calculations 
(shown in Figure A). In addition, pressure 
measurements made in the wake behind the 
trailing edge of the wing in the center of the 
17-foot test span (shown in Figure B) con- 
firmed a major reduction in wing drag over a 
wide variety of test conditions and showed 
that large amounts of laminar flow over the 
aircraft’s wing box area had been achieved as 
indicated by the early hot-film data. 

Status/Plans. Flight testing is continuing. 
Detailed experiments will be made to esti- 
mate the tolerance of extensive laminar flow 
to relaxation of the stringent laminar flow 
design requirements. 

Dal V. Maddalon 
Flight Applications Division 
Langley Research Center 
(804) 864-1909 


3 



CP 




1-2 \fortex Generators for High-Lift Airfoils 


Objective. To evaluate the effectiveness of 
very small vortex generators for controlling 
boundary layer separation over three-element 
high-lift airfoils. 

Approach. Small, vane-type vortex genera- 
tors were mounted on a Douglas three-ele- 
ment, two-dimensional model in landing con- 
figurations as part of a cooperative test pro- 
gram between NASA and Douglas Aircraft 
Company conducted in the Langley Low-Tur- 
bulence Pressure Tunnel (LTPT). The inves- 
tigation was conducted at near-flight condi- 
tions with chord Reynolds numbers of 5 mil- 
lion and 9 million and a Mach number of 0.2. 

Accomplishments. Separation control with 
vortex generators was investigated during 
the second and third quarter of FY 1991 on the 
Douglas high-lift model. Measurements in- 
clude lift, drag, surface pressure, wake pro- 
files, and fluctuating shear stress. Low-pro- 
file vortex generators mounted at 25% of the 
flap chord, and with a device heightof only 0.2 
inches at full scale, were effective in alleviat- 
ing flap separation. The resulting wake of the 
three-element airfoil was significantly nar- 
rowed, providing measured drag reductions 
of up to 50%, (see figure). In addition, lift 
coefficients at approach angles-of-attack (i.e. , 
8° to 12°) were increased by up to 10%. 

Significance. The results of this investiga- 
tion have indicated: (1) Separation allevia- 
tion on the flap can significantly improve the 
performance of the entire high-lift system. (2) 
A practical vortex-generator system that not 
only can produce effective separation allevia- 
tion but also produce low drag when deployed, 
no drag when retracted in cruise (stow inside 
the flap well), and simplicity and ease of 
maintenance. (3) Low-profile, flap mounted 
vortex generators can be integrated into a 
modern transport to maximize performance 
and eliminate buffeting. 


Status/Plans. More vortex generator tests 
are planned for 1992, as well as a vortex 
generator jet test on a NASA high-lift model. 

John C. Lin 

Fluid Mechanics Division 
Langley Research Center 
(804) 864-5556 
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DRAG POLARS FOR SOLID AND 
POROUS NACA 0012 AIRFOILS 


Mo, = 0.80 



SELF-ADAPTIVE AIRFOIL 
THROUGH POROSITY 


Porous NACA 0012 

M<o = 0.80, a « 1 .25° i Solid 

Mo, = 0.63, a = 2.00°/ equivalent 



Figure 1-3. Porosity for Transonic Airfoils With a View Toward Multipoint Design 
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1-3 Porous Transonic Airfoils for Multipoint Design 


Objective. To determine the feasibility of 
using porosity to achieve multipoint design 
for transonic airfoils. 

Approach. A computational pilot study was 
conducted with an in-house developed time- 
implicit upwind method for the nonconserva- 
tive Euler equations using floating shock fit- 
ting for accurate representations of shocks. 
The computational approach allowed for a 
faster and less-expensive assessment of influ- 
encing parameters, such as profile geometry, 
arrangement of porous surface patches, and 
degree of porosity, than an experimental study. 

Accomplishments. Solutions were com- 
puted for steady transonic flow over NAGA 
0012 and supercritical airfoils with solid as 
well as porous surfaces. The porous surfaces 
were applied to both upper and lower profile 
surfaces, and they extended over the nominal 
chord. Either connected or separated cavities 
were assumed to lie beneath the upper and 
lower surfaces. The porosity distribution is 
described by a modified sine wave with sev- 
eral amplitudes. The sinusoidal porosity dis- 
tribution was chosen to facilitate the numeri- 
cal approximation of flow through permeable 
surfaces. 

Applied to a NAC A 0012 airfoil, porosity gen- 
erally increases lift for a given angle-of-at- 
tack, in some instances by up to 65%. As 
indicated by the accompanying figure, the 
wave drag taken at constant lift for 
supercritical flow past a porous NACA 0012 
airfoil with separated cavities is up to one 
order of magnitude lower than for its solid 
counterpart. 

The accompanying figure also illustrates the 
new quality of a porous NACA 0012 airfoil 
being self-adaptive to dissimilar flow condi- 
tions. Using the calculated porous surface 
pressure distributions as target pressures, 
equivalent solid airfoil shapes were con- 


structed using a computational design tool. 
While the subcritical equivalent airfoil basi- 
cally retains the teardrop shape of the NACA 
0012 profile, the supercritical companion piece 
reveals a distinct flattening of the upper sur- 
face combined with hump-shaped closure to- 
ward the trailing edge. 

Making the surfaces of an already optimized 
supercritical airfoil permeable also broadens 
its operational speed and incidence range. 

Significance. Euler analyses indicated a 
potential of using passive venting techniques 
for restructuring the flow past airfoils such 
that they become self-adaptive to very dis- 
similar flow conditions. This property could 
open an alternate route for achieving 
multipoint design, that is, the design of air- 
foils which satisfies several, oftentimes con- 
flicting constraints. 

Status/Plans. The proposed concept will be 
experimentally verified in the 8-Foot Tran- 
sonic Pressure Tunnel at Langley Research 
Center. 

Peter M. Hartwich (ViGYAN Inc.) 

Applied Aerodynamics Division 
Langley Research Center 
(804) 864-2881 
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ORIGINAL PAGE 

SLACK AND WHITE PHOTOGRAPH 


Flight Results on a Subsonic Transport Flap System 


Pressure Distributions on 1st Flap Element 


Flap deflection = 40 deg 

Cl= 2.46 Airspeed = 101 kts (Indicated) 

Angle of attack = 10.21° RN, c = 9.96 x 10 s 



Tuft Flow Visualization of the Flap System 



Figure 1-4. Flight Results on a Subsonic Transport Flap System 
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1-4 Subsonic Transport High-Lift Flight Research - Phase I 


Objective. To obtain in-flight aerodynamic 
flow measurements on subsonic transport 
high-lift systems for correlation with wind- 
tunnel data and validation of CFD results to 
improve high-lift design methodology. 

Approach. As part of a multiphase research 
activity, Phase I flight experiments were con- 
ducted on a B737-100 (NASA 515) research 
aircraft to measure pressure distributions on 
the high-lift flap system of a transport con- 
figuration. Phase I activities included mea- 
suring pressure distributions at select 
span wise locations of the aircraft flap system. 
Pressure data were obtained in flight using 
belts of thin plastic tubing wrapped around 
the upper and lower surfaces of each of the 
trailing-edge flap elements at two spanwise 
locations. In addition, surface tufts were used 
to document flow-separation characteristics, 
and Preston tubes were used to measure sur- 
face shear-stress characteristics at selected 
locations. 

Accomplishments. Phase I flight tests on 
the NASA 515 research aircraft were com- 
pleted in May 1991. Steady-state test points 
(about 250 test conditions) were obtained for 
a range of altitudes up to 20,000 feet and 
airspeeds as low as the stick shaker speed. 
Chord Reynolds numbers ranged from 10 to 
20 million and Mach numbers ranged from 
0.16 to 0.40. Pressure data was obtained on 
each flap element of the triple-slotted flap 
system for deflections from 15° to 40°. Flow- 
separation characteristics on the flap system 
were documented with surface tufts, and sur- 
face shear-stress measurements were deter- 
mined from Preston tube readings. 

Significance. An aerodynamic database in- 
cluding pressure distributions and documen- 
tation of flow separation was established on a 
complex high-lift flap system in flight over a 
range of Reynolds and Mach numbers and 
angle-of-attack conditions. These data will 


contribute to the improved understanding of 
high-lift flow physics and will provide a good 
test case for correlation of wind tunnel data 
and validation of computational fluid dynamic 
results. 

Status/Plans. Additional flight tests (Phase 
IA) are planned for winter 1991/1992 to ob- 
tain full-chord pressure distributions along 
the wing and slat elements as well as the flap 
elements. Analysis of flight data will be made 
using two and three-dimensional computa- 
tional methods to help plan a follow-on Phase 
II flow physics flight experiment. Phase II 
activities will incorporate the use of replace- 
ment 737 slat and flap spare parts to house 
instrumentation for detailed surface pres- 
sure distributions, boundary layer, and wake 
flow measurements. 

Long Yip, Paul Vijgen 
Flight Research Branch 
Langley Research Center 
(804) 864-3866 
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Chapter 2 

Test Techniques and Instrumentation 


Technology is being provided for critical experimental research required to improve the measure- 
ment of the fundamental flow properties of fluids and the overall aerodynamic performance of 
aircraft components and configurations. The development of instrumentation and measurement 
techniques for real-time, flow diagnosis is being performed with emphasis on nonintrusive 
methods. These developments occur across the range of conditions from cryogenic to high 
temperature and from low subsonic to hypersonic speeds. 

One of the more excitinginstrumentation developments underway is the use of pressure sensitive 
luminescent paint on aerodynamic models. This will provide nonintrusive, inexpensive sensing 
of pressure in real time for wind tunnel and flight testing. 

Other research areas being addressed are (1) heavy gas wind tunnel testing to operate at high 
Reynolds numbers; (2) solid-state camera/optical image correction techniques to correlate flight 
test data of aircraft encounters with wing tip vortices generated by large aircraft; (3) nonintrusive 
Rayleigh-Raman multipoint measurements of gas density in hypersonic flow; (4) liquid crystal 
skin friction sensors capable of high-frequency response in high-speed flows; and(5) nonintrusive 
infrared thermography of global heat transfer rates or wind tunnel models in hypersonic flow. 

Program Manager: Gary Hicks 

OAST/RF 

Washington, DC 20546 
(202) 453-2830 
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Supersonic Oblique Wing Test 
NASA Ames 9X7 Unitary Wind Tunnel 
View of model lower surface - normal light 
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Figure 2-1. Luminescent Paint Sensor Development 
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2-1 Luminescent Paint Sensor Development 


Objective. To develop luminescent “paints” 
for the measurement of surface pressure and 
temperature fields in aerodynamic testing. 
Our ultimate goal is to develop this methodol- 
ogy into a practical tool such that it can 
replace conventional methods. 

Approach. This work is a cooperative re- 
search program between the Ames Fluid Me- 
chanics Laboratory and the University of 
Washington Chemistry Department. Candi- 
date luminescent molecule/coating formula- 
tions were identified, and trial coatings were 
made. The suitability of these trial coatings 
was assessed through static and wind tunnel 
tests. 

Accomplishments. In May 1991 the super- 
sonic oblique wing test was being conducted 
in the Ames 9 by 7-foot Unitary Plan Wind 
Tunnel. The pressure sensitive luminescent 
paint test was piggybacked onto the larger 
experiment, whose objective was to evaluate 
the performance and surface flow character- 
istics of a generic oblique wing model for Mach 
1.6 to 2.0. Figure A is a photo of the model in 
the test section showing the model’s lower 
surface. The lower surface of the wing was 
coated with the pressure sensitive paint. 
Because of the extreme test conditions, a low 
test section static pressure and a high total 
temperature, it was not possible to determine 
the absolute pressure level using the paint. 
Relative pressure change, however, could be 
seen. 

Figure B is a representative example of the 
data acquired. Shown is a map of the lower 
surface pressure field over the forward swept 
portion of the oblique wing. Conditions are 
notedin the figure. The unique field measure- 
ment capability of the paint method is evi- 
dent: the paint capturing the passage of the 
fuselage bow shock over the outboard portion 
of the wing and the passage over the wing 
inboard portion of a shock arising from reflec- 


tion of the wing tip shock off the fuselage. 
More subtle features such as the formation of 
streamwise vortices over the wing surface 
from the cross-flow and the subsequent sepa- 
ration can also be seen. 

Significance. The oblique wing test repre- 
sents the first use of the pressure paint under 
supersonic flow conditions. The flow condi- 
tions for this test were outside the previous 
operational experience with the pressure paint 
and were considered marginal at best. How- 
ever, the paint did provide qualitative infor- 
mation, From this experience we now have 
more information on the practical operating 
range of the pressure paint and the experi- 
mental methodology associated with its use 
in a large-scale supersonic facility. 

Status/Plans. Further development of the 
pressure sensitive paint to improve its char- 
acteristics continues. Near-term work is fo- 
cusing on the development of a temperature 
sensitive luminescent paint and a dual sensor 
paint that will allow simultaneous measure- 
ment of pressure and temperature. We are 
also pursuing the development of lumines- 
cent sensor coatings for unsteady aerody- 
namic applications. 

B.G, McLachlan, J. Bell, 

J. Schreiner 

Fluid Dynamics Research Branch 
Advanced Aerodynamic Concepts Branch 
Ames Research Center 
(415) 604-4142 
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Figure 2-2. Laser-Induced Fluorescence Instrumentation in the Ames 3.5-Foot Hypersonic 

Wind Tunnel 





2-2 Temperature and Density Measurements in Air Using Laser- 

Induced Fluorescence 


Objective. To provide nonintrusive mea- 
surements of flow-field temperature, density, 
and their fluctuations owing to turbulence 
with application capability in all air flows 
including hypersonic wind tunnel flows. 

Approach. Laser-induced fluorescence of 
oxygen is used in combination with Raleigh 
scattering from the same laser beam to pro- 
vide instantaneous and simultaneous mea- 
surements of temperature and density from 
each laser pulse. 

Accomplishments. Bench-top experiments 
have demonstrated that measurements can 
be made at the free-stream conditions of the 
Ames 3.5-foot hypersonic wind tunnel at 
Mach 10 with uncertainties of less than 2%. 
Initial measurements have been made in the 
3.5-ft Wind Tunnel at Mach 7. 


Significance. Nonintrusive measurements 
of flow-field properties and their fluctuations 
in hypersonic flows will provide significant 
new information for turbulent model de- 
velopment and the validation of 3-D numeri- 
cal simulation codes that is not obtainable by 
other means. 

Status/Plans. The instrumentation has been 
installed in the Ames 3.5-foot hypersonic wind 
tunnel and preliminary measurements have 
been demonstrated. Further refinement of 
the instrumentation is in progress. 

Robert L. McKenzie 
(415) 604-4749 
Douglas Fletcher 
(415)604-5244 

Experimental Fluid Dynamics Branch 
Ames Research Center 
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Figure 2-3. Application of Thermographic Phosphor Technique 
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2-3 Phosphor Thermography Technique 


Objective. To extend the capability of the 
relative-intensity, two-color thermographic 
phosphor technique to provide global, quanti- 
tative heat transfer data on models in hyper- 
sonic wind tunnels. 

Approach. Model surface temperature time 
history was accurately measured, pixel by 
pixel, prior to and during injection of the 
model through the nozzle boundary layer into 
the test core and following flow establishment 
of freestream flow about the model. One- 
dimensional conduction heat transfer equa- 
tions were used for each pixel on the model 
within the field-of-view of the camera (maxi- 
mum of 262,000 pixels for each time) to deter- 
mine the quantitative value of heat transfer 
rate. These values of heating were compared 
with benchmark values obtained with thin- 
film resistance gages for the same model ge- 
ometry at the same flow conditions and atti- 
tude. 

Accomplishments. Global heat transfer 
rates were measured on a ceramic “modified 
orbiter" model at Mach 10 in air using the 
relative-intensity, thermographic phosphor 
technique. Surface temperature time histo- 
ries were determined during the model injec- 
tion process, and newly developed software 
was used to reduce data to values of heat 
transfer rate. These results were compared 
with heating distributions measured with thin 
film gages on the same configuration at the 
same flow conditions (referred to as HALIS 
orbiter CFD code calibration study). Data 
from these two techniques along the model 
centerline and several span wise cuts showed 
good to excellent agreement. 

Significance. Phosphor thermography has 
evolved from a thermal mapping technique to 
a viable quantitative heat transfer measure- 


ment technique. Global heating distributions 
can be obtained minutes after a wind tunnel 
run, allowing researchers to optimize (i.e., 
tailor) investigation run by run. Compared 
with the thin-film gage technique (recognized 
as the most accurate method for measure- 
ment of heat transfer rate), phosphor ther- 
mography provides invaluable global infor- 
mation, and the time and cost to construct 
ceramic models are an order of magnitude 
less than the time and cost for thin-film mod- 
els. This phosphor technique has revolution- 
ized aerothermodynamic testing at LaRC and 
has generated considerable interest in the 
national/international hypersonic-testing 
community. 

Status/Plans. A next-generation phosphor 
system, which will provide tracking during 
the model injection process and accurate spa- 
tial definition including regions of high sur- 
face curvature, has been specified and is be- 
ing acquired. Advanced methods of applying 
phosphors are being examined and improved 
data reduction techniques are under develop- 
ment. Potential application to impulse facili- 
ties having run times less than 50 msec will be 
investigated. 

N. Ronald Merski 
Experimental Hypersonics Branch 
Langley Research Center 
(804) 864-7539 
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Holographic Interferometry 
Layout at EAST Reflected Shock Tunnel 



Figure 2-4. Density Field Measurements in a Nonequilibrium Expanding Flow Using 

Holographic Interferometry 
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2-4 Density Field Measurements in a Nonequilibrium Expanding 
Flow Using Holographic Interferometry 


Objective. To characterize the flow field of a 
two-dimensional nozzle using holographic in- 
terferometry as a first step towards the inves- 
tigation of nonequilibrium expanding flows. 

Approach. The 10-cm drive section of the 
Ames's electric arc-driven shock tube facility 
was converted into a reflected shock tunnel 
facility by installing a 2-D nozzle insert in the 
tube. A 3.2 km/sec shock in 100 Torr nitrogen 
reflected off the nozzle insert, providing a 
reservoir of compressed test gas with P5-1500 
psi and Ts=7000 K. A Nd:YAG laser doubled 
to 532 nm was used to produce holograms of 
the nonequilibrium expanding flow. 

Accomplishments. A diagnostics system 
capable of producing a snapshot (6-ns expo- 
sure time) single-plate, double-exposure ho- 
logram has been developed. Holograms of the 
expanding flow, including startup flow, at 
intervals of 9, 17, 20, 40, and 60 ps from the 
time of shock arrival have been recorded. 
Static pressures at 10 points in the wall of the 
nozzle have also been recorded. 

Significance. The relaxation time scale in 
expanding flows can be 1-2 orders of magni- 
tude lower than the time behind a normal 
shock. To develop a better understanding of 
this nonequilibrium in expanding flows, the 
vibrational populations during the relaxation 
will be measured. However, to accomplish 
this, the bulk density in the flow field must be 
known. The holograms provide the bulk den- 
sity data. 


Status/Plans. In addition to these holo- 
grams, synthetic holograms using the com- 
puted flow-field data will be generated. The 
experimentally obtained holograms then will 
be compared with the synthetic ones. Such an 
exercise will enable us to test our ability to 
compute the expanding flows. Also, the den- 
sity flow-field data from the holograms will be 
used in the next stage of the experimental 
program on the investigation of 
nonequilibrium expanding flows. 

Surendra P. Sharma, Scott A. Meyer, 

Walter D. Gillespie 
Aerothermodynamics Branch 
Ames Research Center 
(415)604-3432 
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2-5 Measurement Techniques for Hypersonic Flows 


Objective. To evaluate and develop experi- 
mental methods for use in scheduled Generic 
Hypersonics and N ASP-related programs. 

Approach. Two experimental techniques for 
hypersonic flows in helium were examined: 
(1) An exploratory test was made to deter- 
mine if infrared measurements could be used 
to map the two-dimensional temperature field 
on the surface behind a series of vortex gen- 
erators. Prior to the test, it was not known if 
the infrared imaging system would be sensi- 
tive enough to accurately determine the loca- 
tion of vortex-induced surface heating behind 
vortex generators. (2) Measurement of the 
constant-temperature hot-wire anemometer 
characteristics was made, and examination of 
hot-wire data reduction techniques in hyper- 
sonic flow was begun by a summer research 
associate, Professor Eric Spina of Syracuse 
University in collaboration with Catherine 
McGinley of Experimental Flow Physics 
Branch, Fluid Mechanics Division. 

Accomplishments. Despite the very low 
heat transfer rates on the surface of the base 
model in the High Reynolds Number Helium 
Tunnel, the wakes behind a series of vortex 
generators were visible, allowing the assess- 
ment of wake strength as a function of genera- 
tor angle and orientation. A frequency re- 
sponse of the hot-wire system was found to be 
adequate for measuring turbulence quanti- 
ties at Mach 11. 

Significance. Infrared measurements will 
enable a relatively simple determination of 
the geometry, size, and spacing of vortex gen- 
erators needed to trip a hypersonic boundary 
layer with minimum drag penalty. With a 
constant-temperature anemometer system, 
measuring mass-flow fluctuations in hyper- 
sonic boundary layers with adequate fre- 
quency response will be possible. These tests 
showed that a constant current system is not 
necessary. 


Status/Plans. Infrared tests will be contin- 
ued, and a parametric study of vortex genera- 
tor geometry will be made to determine the 
most effective configuration for hypersonic 
flow. Measurements of mass flux fluctuations 
in a hypersonic boundary layer will be made. 
Efforts to improve data reduction methods 
will continue. 


Ralph Watson 

Experimental Flow Physics Branch 
Langley Research Center 
( 804 ) 864-5723 . 
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Chapter 3 

Transition and Turbulence Physics 


The objective of the Transition and Turbulence Physics Program is to develop a fundamental 
understanding of flow structures relating to transition and turbulence and to incorporate these 
flow structures into sophisticated flow models for use with computational methods. Extensive 
efforts are being developed within the program to solve the full Navier-Stokes equations, which 
include transitioning and turbulent flows. However, the realities of today's supercomputer 
limitations in terms of memory and speed make it impractical to consider full Navier-Stokes 
solutions for all but the simplest problems. Practically speaking, research in transition and 
turbulence modeling with various levels of approximation is providing the bridge for the gap until 
the computer power is available. The understanding of the physical flow structures of 
transitioning and turbulence flows is being examined from two directions, computational and 
experimental. Extensive advancements in test instrumentation being developed within other 
discipline programs is enhancing transitional and turbulence experimental research. The 
development of quiet supersonic wind tunnels will provide an experimental tool for investigating 
supersonic viscous flows. Exciting advances are occurring in compressible boundary layer 
transition with the use of parabolized stability equations derived from the complete Navier- 
Stokes equations. The program supports the Center for Turbulence Research which assembles 
world renowned experts in the field to generate new concepts regarding turbulence. A similar 
effort is being established for fluid mechanics research emphasizing transition. 

Program Manager: Gary Hicks 

OAST/RF 

Washington, DC 20546 

(202) 453-2830 
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3-1 Compressible Boundary Layer Transition Using PSE 


Objective. To understand the mechanisms 
involved in compressible boundary layer tran- 
sition, specifically to provide a capability for 
boundary layer transition prediction in both 
“quiet” and “disturbed” environments. 

Approach. For convective instabihties such 
as T-S waves in boundary layers, the govern- 
ing PDEs are only weakly elliptic along the 
dominant flow direction. The original PDEs 
are reduced to a set of parabolized stability 
equations (PSE), which is parabolic along the 
dominant flow direction, allowing solution by 
single-sweep marching. Both nonparallel and 
nonlinear effects of growing boundary layers 
are studied up to the transition stage using 
this new PSE approach. 

Accomplishments. The PSE code for three- 
dimensional linear and nonlinear disturbances 
was developed for flat-plate geometry. Linear 
PSE calculations were performed for Mach 
numbers 0 to 4.5. The effect of boundary layer 
growth was shown to be important for both 
the first- and second-mode disturbances in 
supersonic boundary layers. Nonlinear PSE 
calculations show very good agreement with 
earlier temporal direct numerical simulations 
(DNS) at Mach 4.5 for the nonlinear evolution 
of a second-mode disturbance. The computa- 
tional time required is an order of magnitude 
less than that for DNS. 

Part (a) of the figure compares the growth 
rates of a linear disturbance as computed by 
the parallel PSE (solid line), the nonparallel 
PSE (dashed fine), and the multiscale DNS 
methods. The results indicate that nonparallel 
effects are more important for oblique (y = 
50°) waves than for 2-D (y = 0°) waves. A 
likely route to transition may consist of two 
oblique primary waves interacting with each 
other. Part (b) of the figure shows the oblique 
breakdown procedure for three different ini- 
tial amplitudes at the given frequency. Tran- 
sition is clearly located at the abrupt rise of 


wall shear. The results indicate that rapid 
growth of the secondary disturbance is trig- 
gered when the primary amplitude is suffi- 
ciently high. 

Significance. The PSE method can be used 
to study the nonparallel and nonlinear evolu- 
tions of disturbances in compressible bound- 
ary layers starting from the linear stage up to 
transition. Coupled with receptivity phenom- 
ena this approach offers a computationally 
viable means for studying and predicting the 
complex and intricate phenomenon of com- 
pressible boundary layer transition. 

Status/Plans. An axisymmetric version of 
the PSE will be developed and used for tran- 
sition studies in conical flows. The PSE theory 
will be extended to 3-D boundary layers. Work 
is also under way on formulating the leading 
edge receptivity problem using the PSE ap- 
proach. 

C. L. Chang 

Fluid Mechanics Division 
Langley Research Center 
(804) 864-5563 
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Figure 3-2. Large-Eddy Simulation of Compressible Isotropic Turbulence 
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3-2 Compressible Large-Eddy Simulation of Isotropic Turbulence 


Objective. To assess the subgrid-scale (SGS) 
model proposed for compressible large-eddy 
simulation (LES) by Speziale, Erlebacher, 
Zang, and Hussaini. 

Approach. The assessment was performed 
by conducting actual large-eddy simulations 
of compressible, isotropic turbulence and com- 
paring the results with direct numerical simu- 
lation (DNS) data. 

Accomplishments. An extensive new set of 
DNS data on 96^ grids was generated. These 
data covered a range of fluctuating Mach 
numbers (M) and fraction of compressible 
energy (x), the key parameters in compress- 
ible turbulence, according to a recent theory 
developed at ICASE. The proposed SGS model 
was utilized in LES for all of these cases. The 
figure compares the time evolution of the 
kinetic energy of decaying, isotropic turbu- 
lence from three different runs - a fine-grid 
(96^) DNS, a coarse-grid (32^) DNS, and a 
coarse-grid (32^) LES. As shown in the figure, 
coarse-grid LES is quite accurate on the total 
kinetic energy (E), but is somewhat less accu- 
rate for the compressible component of the 
kinetic energy (Ec). Moreover, the coarse- 
grid LES does substantially better than a 
coarse-grid DNS. Surprisingly, the model 
predictions for the compressible component 
sometime improve as the flow becomes more 
compressible. Some of the problems can be 
ameliorated by refining the constants in the 
SGS model. However, subtle changes in the 
model itself appear necessary to obtain more 
consistent predictions in the compressible 
regime. 

Significance. These were the first system- 
atic, compressible LES calculations that have 
been performed. Previous tests of the SGS 
model were based solely on DNS data and did 


not include actual large-eddy simulations. 
The results show that this first-cut SGS model 
for compressible turbulence performs better 
than anticipated. 

Status/Plans. The results of these first com- 
pressible LES calculations have sparked con- 
siderable interest within the Langley group. 
Several improvements to the compressible 
part of the SGS model have already been 
suggested. The improved model will be ap- 
plied to more challenging flows, such as com- 
pressible turbulence containing eddy 
shocklets, homogeneous turbulence in uni- 
form shear flow, and wall-bounded shear flows. 

Thomas A. Zang 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-2307 
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Figure 3-3. Effect of Compressibility Corrections on Vorticity Thickness Growth Rate 
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3-3 Compressible Turbulence Modeling for High-Speed 

Shear Layers 


Objective. To assess the performance of 
several compressibility corrections to turbu- 
lence models applied to high-speed shear lay- 
ers. 

Approach. The compressible Navier-Stokes 
equations were numerically solved using a 
two-equation turbulence model with and with- 
out compressibility corrections. Numerical 
results were compared with available experi- 
mental data to determine the most appropri- 
ate model. 

Accomplishments. Each of the compress- 
ibility corrections studied was developed to 
increase the dissipation in the kinetic energy 
of the turbulence in a shear layer. By lower- 
ing the eddy viscosity through a reduction in 
the turbulence energy, these models reduced 
the spread rate of free-shear layers for a wide 
variety of flow conditions as the convective 
Mach number increased. 

Significance. The vorticity thickness growth 
rate comparisons between computation and 
experiment indicate that the application com- 
pressibility corrections can significantly im- 
prove predictions of high-speed shear layer 
mixing. 


Status/Plans. This study is ongoing. Recent 
results were reported in June 1991 at the 
AIAA 22nd Fluid Dynamics, Plasma Dynam- 
ics and Lasers Conference. Further studies 
will include testing of additional compress- 
ibility corrections and detailed comparisons 
with shear layer experiments. 

John R. Viegas 
Fluid Dynamics Division 
Ames Research Center 
(415)604-5950 
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3-4 Development of a Mach 18 Quiet Helium Ttmnel 


Objective. To modify the Mach 20/40 Open 
Jet Helium Tunnel to operate as a Mach 18 
Quiet Helium Tunnel. 

Approach. Only the minimum essential 
modifications required to convert the existing 
conventional facility to a quiet flow tunnel 
were made. The settling chamber was outfit- 
ted with freestream turbulence and noise 
control devices. These consisted of a new inlet 
diffuser, two high-density woven wire porous 
plates, and two high open-area screens. The 
existing settling chamber was modified to 
accept a slotted nozzle and associated bound- 
ary layer bleed system piping. The down- 
stream portion of the existing Mach 20 conical 
nozzle was retained to eliminate the possibil- 
ity of nozzle wall transition due to the Gortler 
instability mode. The throat section, how- 
ever, was replaced with a highly polished unit 
with an annular sonic orifice upstream of the 
throat for removal of the settling chamber 
boundary layer. 

Accomplishments. Engineering design of 
the tunnel modifications was completed in FY 
1989. Manufacture and acquisition of the 
new hardware was completed dining FY 1990. 
Installation of the new components and asso- 
ciated piping has occupied most of FY 1991. 

Significance. Because of the high radiated 
noise environmentof conventional high-speed 
wind tunnels, transition experiments with 
known noise sensitivity or experiments where 
the noise sensitivity has not been established 
must be conducted under quiet flow condi- 
tions. Virtually all hypersonic transition tests 
fall in the latter category. If the modifications 
prove to be successful, this facility will offer 
the first hypersonic, quiet transition testing 
capability above Mach 6. 


projected for October 1991. Initial runs will 
use an existing, unpolished nozzle throat sec- 
tion and will focus on ensuring that the sys- 
tem helium flow is free of particles that may 
damage the highly polished nozzle finish. In 
the event that particles are found to be present, 
plans to install a filter upstream of the set- 
tling chamber are in preparation. 

Stephen P. Wilkinson 
Fluid Mechanics Division 
Langley Research Center 
(804)864-5733 


Status/Plans. At the present time, certifica- 
tion of the new high-pressure piping is 90% 
complete. Start-up and shake-down runs are 
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3-5 Direct Simulation of Compressible Homogeneous Shear Flow 


Objective. To increase the physical under- 
standing of the effects of compressibility on 
flow turbulence so that these effects can be 
modeled in flow computations. 

Approach. Homogeneous compressible tur- 
bulence was considered wherein a linear mean 
velocity field sustains the random turbulence 
field. A spectral collocation method, along 
with a third-order Runge-Kutta time advance- 
ment, was used to perform direct numerical 
simulations (DNS), which are highly resolved 
in space and time. The resulting database 
was analyzed for statistical and structural 
features of the flow. 

Accomplishments. Figure A shows the evo- 
lution of turbulent kinetic energy, K, as a 
function of time, 1, nondimensionalized by the 
shear rate, S, for several values of Mt, the 
Mach number of the turbulent velocity fluc- 
tuations. The figure shows that the level of 
turbulent kinetic energy decreases with in- 
creasing values of initial turbulent Mach num- 
ber, MtK>. Thus, compressibility decreases 
the growth rate of the turbulence intensities. 
The exact equation for the evolution of turbu- 
lent kinetic energy in homogeneous flows is 
d(pK) / dt = pP - pe$ - pec + p’d’, where the 
production, pP, and solenoidal dissipation, £s, 
are supplemented by two compressibility 
terms: compressible dissipation, pec, and pres- 
sure-dilatation, p’d’. Figure B shows that 
these compressibility terms alter the balance 
between production and dissipation by as 
much as 25% for the case with Mt,o = 0.4, 
suggesting that incompressible turbulence 
models are inappropriate. Physically, the 
augmented dissipation that is due to com- 
pressibility and the transfer of energy from 
the velocity field to the pressure field via the 
pressure-dilatation term contribute to the 
reduced growth of K A previous theoretical 
study had shown that the compressible dissi- 
pation obeys the scaling es/ec = M%. Figure C 
shows that for a variety of DNS cases, the 


ratio EctesM^t) approaches an equilibrium 
value of approximately 0.5 for large St. 

Significance. Simulations of homogeneous 
shear flow show that the growth of turbulent 
kinetic energy decreases with increasing Mach 
number - a phenomenon that is similar to the 
experimentally observed reduction of turbu- 
lent intensities in the supersonic shear layer. 
The numerically generated databases have 
exhibited great potential for developing mod- 
els for terms (e.g., the compressible dissipa- 
tion) important in high-speed flows. 

Status/Plans. Further analysis of the data- 
base will be carried out to improve the under- 
standing of compressibility effects and aug- 
ment our present turbulence modeling capa- 
bilities. 

Sutanu Sarkar 
ICASE 

Langley Research Center 
(804) 864-2194 
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FIGURE 2. Time development of the plane-averaged wall shear stress <r w > in Re = 

8000 transitional channel flow, * DNS (Zang et ah 1990); present results; 

LES (Piomelli & Zang 1990b). 


Figure 3-6. Dynamic Subgrid Scale Modeling and the New LES Program 
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3-6 Dynamic Subgrid Scale Modeling and the New LES Program 


Objective. Subgrid scale models used in 
large eddy simulations of turbulent flows have 
had several drawbacks: the inability to cor- 
rectly represent, with a single universal con- 
stant, different turbulent fields in rotating or 
sheared flows, the requiring of ad hoc damp- 
ing functions near solid walls, and the requir- 
ing of ad hoc intermittency functions to pre- 
dict the transition from laminar to turbulent 
flow. A new eddy viscosity model has been 
developed that alleviates many of these draw- 
backs. At present, our objective is to evaluate 
this model in a variety of incompressible, 
compressible, transitional, and fully devel- 
oped turbulent flows. 

Approach. In large eddy simulations (LES) 
the large-scale field is directly computed and 
the effects of small scales are modeled. In the 
dynamic mode approach, the coefficient of the 
subgrid scale eddy viscosity is a function of 
space and time, and hence varies in the differ- 
ent flows and flow regimes. The coefficient is 
computed dynamically dining the computa- 
tion rather than input a priori. This model 
was used in large eddy simulations of transi- 
tional and turbulent channel flow with heat 
transfer and in compressible homogeneous 
turbulence. 

Accomplishments. The concept of dynamic 
modeling was used to derive expressions of 
the subgrid scale eddy viscosity and turbulent 
Prandtl number. The results were in good 
agreement with the experimental data and 
direct numerical simulations, and were bet- 
ter than those of LES using conventional 
models that include ad hoc intermittence and 
damping functions. The model has performed 
remarkably well without any adjustments in 
very different flow situations. 

Significance. The idea of dynamic modeling 
is an altogether new approach in turbulence 
modeling. It has led to physically sound 
models for momentum and energy transport. 


Its inherent versatility should facilitate its 
use in applied CFD. 

Status/Plans. The dynamic model will be 
applied to complex flow configurations. These 
include flow over a cylinder, backward facing 
step, separated flat plate boundary layer, 
compressible flow transition on a flat plate, 
and flow over a concave wall. 

Parviz Moin 

Center for Turbulence Research 
Ames Research Center and Stanford Uni- 
versity 

(451) 723-9713 or 604-5127 
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3-7 Efficient Supersonic Wind Tunnel Drive System for Transi- 
tion Research at Mach 2.5 


Objective. To validate design principles for 
a new generation of low-cost Laminar Flow 
Supersonic Wind Tunnels (LFSWTs) using 
“Quiet” technology. The LFSWT is essential 
for boundary layer receptivity studies at high- 
speed civil transport (HSCT) Mach numbers. 

Approach. A 1/8 scale Proof ofConcept(PoC) 
model of the LFSWT was built to research 
unknowns for the full-scale design. The pri- 
mary purpose of PoC was to establish that the 
Fluid Mechanics Laboratory (FML) indraft 
compressor could support the desired LFSWT 
test envelope of M = 2.5 @ Re = 1.3 million per 
foot. The low compression ratio (~ 1 .8: 1) of the 
compressor precluded the use of any conven- 
tional tunnel drive system. For the past 19 
months, there has been an experimental re- 
search effort, supported by CFD studies, to 
develop a unique drive system for the LFSWT. 
In parallel, innovative design techniques for 
the “Quiet” settling chamber and supersonic 
nozzle are being investigated. 

Accomplishments. In February 1991, we 
achieved M = 2.5 in the PoC test section over 
the desired Re range. To achieve the lowest 
Re of 1 million per foot, a stagnation pressure 
(Po) of 5 psia, is required. A Po of 5 is less than 
the exit pressure (Pe) of 8 psia, which makes 
the compression ratio across the test section 
0.625:1, which is a major and unique accom- 
plishment. By comparison, a conventional 
supersonic tunnel requires a compression ra- 
tio greater than 2 :1 to run and an overpressure 
to start. The PoC has an advantage in that the 
starting and running test conditions are the 
same. Operation at such a low compression 
ratio is achieved through the use of a super- 
sonic diffuser (optimized by CFD support) at 
the end of the test section and dual ambient 
injector system. The primary injector is de- 
signed for an M = 2.4 exit flow and the second- 


ary for M = 2.0. The injectors pull conditioned 
air through the test section. The ratio of mass 
flow through the primary injector over that of 
the test section is an unprecedented 9.17:1. 

Significance. This low-budget project has 
demonstrated that we can operate a super- 
sonic wind tunnel with a nonspecializedindraft 
compressor. We now have the technology to 
design a drive system for the full-scale LFSWT. 
The LFSWT will enhance supersonic nozzle 
and transition research progress within NASA 
while impacting HSCT. 

Status/Plans. The LFSWT drive system is 
being designed, and the LFSWT design prin- 
ciples are being studied for the settling cham- 
ber and supersonic nozzle, which are critical 
to establishing quiet flow. 

J. Laub, S. Wolf, L. King, D. Reda 
Fluid Dynamics Research Branch 
Ames Research Center 
(415) 604-4136 
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Figure 3-8. Precursor Effect in Hypersonic Transition 















3-8 Numerical Simulation of Laminar Breakdown in Supersonic 

Transition 


Objective. To explain theoretically, the “rope- 
like structures” and the “precursor transition 
effect” that have been observed in supersonic 
transition experiments since the 1960s. The 
rope-like structures refer to the braided pat- 
terns that are visible near the boundary layer 
edge in schlieren photographs. The precursor 
transition effect refers to the origination of 
transitional disturbances near the boundary 
layer edge a considerable distance upstream 
of the first detection of transition by surface 
measurements. 

Approach. The approach was to conduct 
high-resolution direct numerical simulations 
(DNS) of laminar breakdown in supersonic 
boundary layers on cylinders and cones, uti- 
lizing linear stability theory (LST) and sec- 
ondary instability theory (SIT) to select cases 
that represent the most likely paths to transi- 
tion. 

Accomplishments. Two high-resolution 
temporal DNS of the strongly nonlinear, lami- 
nar breakdown stage of supersonic transition 
to turbulence were conducted: for Mach 4.5 
flow past a cylinder and for Mach 6.8 flow past 
a cone. The top frame of the accompanying 
figure shows a schlieren photograph that is 
typical of experimental results for transition 
on a cone in hypersonic flow. The essential 
features are displayed in the cartoon in the 
middle of the figure. The spatially recon- 
structed turbulence intensities extracted from 
the DNS of the Mach 6.8 cone case (bottom 
frame) display these key features. In particu- 
lar, the “precursor transition effect” has been 
captured, with even the spreading angles of 
the disturbance regions in good agreement 
with the experimental range: between 0.5° 
and 1°. Moreover, numerical schlieren flow- 
field visualizations from the DNS display 
“rope-like” structures are remarkably similar 
to those that have been observed experimen- 
tally. 

* 


Significance. A comparison of LST, SIT, 
and DNS results for the spreading angle of 
the disturbances indicates that the spreading 
rate cannot be attributed to linear instability 
alone (as has long been believed), nor even to 
secondary instability, but is, indeed, a strongly 
nonlinear effect. 

Status/Plans. The numerical databases al- 
ready generated will be analyzed in detail, 
with particular emphasis given to their impli- 
cations for transition modeling. 

C.D. Pruett 

Fluid Mechanics Division 
Langley Research Center 
(804) 864-6788 
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Figure 3-9. Receptivity of Low-Speed Boundary Layers 
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3-9 Receptivity of Low-Speed Boundary Layers 


Objective. To develop a research program to 
understand the essential issues underlying 
the receptivity process and to use this knowl- 
edge in developing new approaches to transi- 
tion prediction. Transition can be considered 
a three-stage process consisting of receptiv- 
ity, instability, and breakdown to turbulence. 
While most previous work has been concerned 
with the issue of flow stability, receptivity is 
now considered a critical first-order param- 
eter. 

Approach. A combined experimental/com- 
putational effort is under way to develop a 
receptivity database. Following previous work 
on receptivity to suction surfaces (see FY 1990 
Research Accomplishments), two new experi- 
ments are now under way. One experiment 
concerns the transition of a laminar, low- 
speed boundary layer in an adverse pressure 
gradient. The disturbance created by a puff 
periodically introduced through a small hole 
in the test surface is tracked with an auto- 
mated data acquisition system using phase- 
lock-averaging software. Pattern formation 
in the boundary layer - expressed as velocity 
vectors and vorticity contours - is examined 
using video animations. The second experi- 
ment concerns the receptivity to a roughness 
element and follows a theoretical treatment 
using triple deck theory. This experiment will 
study the nature of roughness-induced recep- 
tivity resulting from freestream acoustic dis- 
turbances. 

Accomplishments. The response of a low- 
speed boundary layer to a periodic jet at one 
phase interval is shown in the left panel of the 
accompanying figure. The minimum local 
boundary layer thickness is approximately 6 
mm. Animated flow sequences clearly show 
the generation of vortical flow structures and 
their downstream convection. The second 
experiment is currently being designed based 
on results of the computational visualization. 
The right panel in the figure shows the 


analytical streamlines induced by the “rough- 
ness element” and the down wind extent of 
this disturbance. 

Significance. This integrated computational/ 
experimental project is designed to under- 
stand the essential features of boundary layer 
receptivity. By collecting a detailed database 
concerning the receptivity process, enough 
knowledge will be gained to incorporate these 
effects in transition prediction schemes. The 
final result will be amplitude- and scale-de- 
pendent criteria for transition. 

Status/Plans. Current work is focused on 
low-speed, incompressible boundary layers. 
Further receptivity studies on suction holes 
are also planned. Once the basic physical 
processes are understood, emphasis will shift 
to high-speed flows. A special interest is in 
receptivity of supersonic boundary layer (1.5 
<L M ^ 2.5) to suction surfaces as are envi- 
sioned for the next-generation supersonic 
transport. 

S. Davis, J. Watmuff, M. Tadjfa, 

D. Reda 

Fluid Dynamics Research Branch 
Ames Research Center 
(415) 604-4197 
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Figure 3-10. Simulation of Homogeneous Turbulence on the Intel i860 Hypercube 
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3-10 Simulation of Homogeneous Turbulence on the 

Intel i860 Hypercube 


Objective. To develop programming and 
compiler techniques to exploit highly parallel 
computers for turbulence research. 

Approach. The three-phase approach was 
first to construct a highly optimizing compiler 
for the i860 microprocessor, then to develop 
and implement efficient inter-processor com- 
munication procedures, and finally to apply 
these to a program and compile an existing 
homogeneous turbulence algorithm. 

Accomplishments. First, the vectoral com- 
piler was ported to the i860 microprocessor 
and the hypercube multiprocessor environ- 
ment. A number of new compiler techniques 
were developed to exploit the high level of 
internal concurrency available in the i860. 
Techniques were also developed to alleviate 
the problem of inadequate raw memory band- 
width. The homogeneous turbulence code of 
Rogallo was then ported from the Cray Y-MP 
to the 128 processor NAS iPSC/860 Gamma 
Prototype. A highly efficient implementation 
of interprocessor communication was added 
to the code. A wall-clock speed of 1.64 Gigaflops 
was obtained for the largest mesh size (2563). 
Excellent agreement with experimental re- 
sults was obtained (see accompanying figure) 
at Reynolds numbers obtainable only at great 
expense on Cray computers. The code has 
also been ported to the 528 processor Caltech 
Delta Machine, which has a mesh rather than 
a hypercube architecture. 

Significance. It is recognized in the indus- 
try that highly parallel, distributed-memory 
computers are the future of supercomputing, 
and it is therefore crucial to begin learning 
how to implement our most computationally 
intensive applications on these machines. This 
project has shown that a typical turbulence 
simulation code can run on such a machine at 


a very high-performance level; in the case of 
128 40MHz i860’s, at a speed exceeding that 
of 10 Y-MP processors. Within the next 2- 
years, machines of the iPSC/860 type will 
have performance in the 50-Gigaflop range, 
far exceeding that of current supercomputers. 

Status/Plans. We will continue to refine the 
i860 Vectoral compiler, since this processor 
and its successors will be used in important 
new multiprocessors, port other turbulence 
codes to the i860 machines to be used as true 
production supercomputers, and port the com- 
piler and turbulence codes to other new mul- 
tiprocessor computers, such as the Stanford 
DASH machine. 

Alan Wray 

Computational Fluid Dynamics Branch 
NASA Ames Research Center 
(415) 604-6066 
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3-11 Compressible TUrbulence: Modeling Rapid Compression 


Objective. To develop new, and improve the 
old, models for compressible turbulent flows. 
A specific goal reported here is investigation 
of the rapid compression effect on turbulence. 
This effect is present in flows with turbulence/ 
shock interactions. 

Approach. In real flows, the individual 
mechanisms and processes of turbulence dy- 
namics cannot be studied in isolation. There- 
fore, the basis for the investigation of the 
rapid compression effect on turbulence has 
been the direct numerical simulations (DNS) 
of compression of homogeneous turbulence 
(Coleman and Mansour, 1991). The study led 
to a development of new turbulence closure 
models that represent the physics of rapid 
compression. 

Accomplishments. The study of the DNS of 
compressed turbulence yielded the following 
important findings: it was established that 
when nearly incompressible turbulence (with 
r.m.s. Mach numbers Mt«l) is rapidly com- 
pressed in one direction (ID), unexpectedly 
high levels of negative pressure-dilatation 
correlation are generated. The pressure-dila- 
tation term (lid) appears in the turbulence 
kinetic energy and its magnitude during the 
ID compression can become by 0(10) larger 
than the total dissipation (et); hence, fid can 
represent a significant loss of turbulent ki- 
netic energy (to pressure fluctuations). The 
striking aspect of this rapid compression 
mechanism is that it is most effective when 
Mt«l, and that it is inefficient when the 
compression is more isotropic, i.e., acting in 
all three directions. All these aspects have 
been included in the new model for pressure 
dilatation. The model-DNS comparison of the 
ratio Ild^t for two values of Mt = 0.02 and 0.44 
are shown in Figure A. Here, the abscissa Aot 
represents the compression time with 1/(1- 
Aot) being the total compression. The differ- 
ence in kinetic energy q2-levels (soley due to 
the pressure dilatation effect) is also pre- 
dicted by the model. 


Significance. In the flows of practical inter- 
est, the condition of directional rapid com- 
pression arises in flow configurations where 
turbulence passes through a shock. A typical 
example of this kind is a compression corner 
flow. The discussed mechanisms of the rapid 
compression effect on the pressure dilatation 
has so far not been considered in turbulence 
modeling. Its significance is apparent in a 
real flow configuration represented by a uni- 
form turbulent flow passing through a normal 
shock wave and shown in Figure B. Here, the 
DNS data (Lee 1991) of the turbulent kinetic 
energy (q2) response to the shock passage are 
compared with a turbulence model results, 
with and without the rapid contribution to the 
pressure-dilatation term; the abscissa xki 
represents the distance normal to the shock. 

Status/Plans. The new model is being tested 
in the compression corner flow configuration; 
rapid compression theory has been devel- 
oped; and modeling of turbulence subjected to 
rapid expansion is also being studied. 

Otto Zeman 

Center for Turbulence Research 
Stanford U/NASA Ames Research Center 
(415) 723-9596; (415) 604-4726 
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Figure 3-12. An Experiment to Guide Turbulence Modeling for Separated Flows 
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3-12 Separating Boundary Layer Experiment for 
Turbulence Modeling 


Objective. To evaluate turbulence models 
used in the prediction of separating flows in 
such applications as high lift airfoils. 

Approach. A computational study was un- 
dertaken to numerically calculate two cases 
of adverse pressure gradient: one with flow 
separation and the other attached. Solutions 
to the incompressible Navier-Stokes equa- 
tions were obtained with INS3D and various 
different turbulence models. Comparisons 
between numerical results and previously 
measured velocities and turbulent Reynolds 
stresses were performed to determine which 
is the most appropriate model to use in ad- 
verse pressure gradients 

Accomplishments. Various turbulence 
models such as the Baldwin-Lomax model, 
the J ohnson-King model, the Wilcox (k-omega) 
model, and the Baldwin-Barth model were 
incorporated into INS3D. Solutions were 
obtained with each of the models for both the 
separating case and the attached case of ad- 
verse pressure gradient. These solutions re- 
vealed that only those models that accurately 
predicted the Reynolds shear stress would 
accurately predict the pressure distribution. 
The numerical results as well as the ex- 
perimental measurements of velocities and 
Reynolds stresses were reported at the June 
1991 AIAA meeting. 

Significance. Accuracy of flow calculations 
involving strong adverse pressure gradient is 
inextricably linked to the accuracy of the 
turbulence model that is employed. Compu- 
tational Fluid Dynamics is being significantly 
advanced by the advent of the Johnson-King 
turbulence model, which was developed at 
NASA Ames Research Center and has been 
further validated in this study. Boeing Com- 
mercial Aircraft Company is actively and 
successfully using the Johnson-King model in 
the design of future aircraft. 


Status/Plans. The experimental database 
will be distributed to other researchers in the 
turbulence modeling community for use in 
model development. 

David M. Driver, Florian Menter 
Experimental Fluid Dynamics Branch 
Ames Research Center 
(415) 604-5396 
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Chapter 4 

Computational Methods and Validation 


The objective of the Computational Methods and Validation Program is to develop and apply 
advanced analysis and computational methods for solving complex, fluid dynamics problems and 
to perform detailed, benchmark experiments using redundant facilities and instrumentation to 
produce high-quality archival data sets to which computational fluid dynamics (CFD) solutions 
can be compared. Areas of interest include modeling turbulence and transition and computing 
complex flows (steady and unsteady, inviscid and viscous) over two-and three-dimensional 
geometries ranging in speed from zero to hypersonic and including such effects as mass inj ection 
and withdrawal. Additional objectives are to (1) demonstrate proof-of-concept computations for 
pioneering applications, (2) disseminate validated computer codes to the aerospace community 
and provide maintenance and consultation on their use, (3) develop innovative techniques for 
scientific visualization of flow-field solutions, and (4) provide an experimental database that is 
taken in the form and detail consistent with CFD modeling requirements and that has 
documented accuracy and limitations of the experimental data. 

The Computational Methods Program is focused on the present and future technology needs of 
the aerospace community. These needs include (1) developing faster and more efficient 
numerical algorithms to facilitate solutions of the full Navier-Stokes equations by large-eddy, 
simulations/small-scale turbulence modeling; (2) developing advanced geometric modeling and 
grid generation techniques for complex, three-dimensional configurations; (3) improving under- 
standing of the effects of grid characteristics on solutions accuracy, convergence, and stability; 
(4) enhancing computational capabilities through development and use of advanced computer 
architectures and expert systems concepts; and (5) developing improved methods for numerical 
simulation of aerothermodynamic flow phenomena associated with hypersonic cruise and 
maneuver vehicles, including real-gas chemistry. 

Experiments are designed for comparison with numerical CFD results in order to ( 1) understand 
flow physics, (2) develop physical models for CFD codes, (3) calibrate CFD codes, and (4) validate 
CFD codes. The experiments range in speed from subsonic to hypersonic and include a variety 
of configurations including generic, fighter/attack, subsonic transport, rotorcraft, ASTOVL, and 
propulsion systems. Work continues in developing high-quality data bases for several classes of 
flows, including (1) high-and low-aspect ratio wings in subsonic and transonic flows; (2) simple 
3-D turbulent flows, including time histories; (3) flow fields about aircraft components; (4) 
propulsive lift flow interactions in ground effect; and (5) unsteady flow interaction in rotor flow 
fields. 

Program Manager: Pamela F. Richardson 

OAST/RF 

Washington, DC 20546 

(202) 453-9857 
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Figure 4-1, Direct Numerical Simulation of Transition and Turbulence in a Spatially 

Evolving Boundary Layer 
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4-1 Direct Numerical Simulation of Transition and Turbulence in 
a Spatially Evolving Boundary Layer 


Objective. To develop a high-order-accu- 
rate, finite-difference method to compute com- 
pressible, transitional, and turbulent flow 
and to demonstrate the capability of the 
method by computing transition to turbu- 
lence on a fiat plate. A study of the physics of 
transition and the creation of a data base for 
turbulence and transition modeling are an 
integral part of this effort. 

Approach. The unsteady, compressible for- 
mulation of the Navier-Stokes equations are 
solved using a high-order-accurate, upwind- 
biased, finite-difference method. This method 
is an extension of an earlier high-order-accu- 
rate method developed to compute incom- 
pressible, turbulent flows. The new code that 
has been developed uses a zonal methodology 
to efficiently process the number of grid points 
used for the computation. 

Accomplishments. The high-order-accu- 
rate method was successfully used to compute 
transition to turbulence on a flat plate in a 
high freestream disturbance environment. 
The numerical results agree qualitatively with 
the available experimental data. Flow visual- 
ization of the computed flow indicated that 
the transition region was characterized by 
detached shear layers and pairs of counter- 
rotating stream wise vortices. The results 
indicate that the essential features of the 
transition process have been captured in the 
computation. Figure A shows instantaneous 
span wise vorticity contours just above the flat 
plate and Figure B shows the computed skin 
friction along the plate compared with experi- 
mental data. 

Significance. This computation is the first 
“transition to turbulence” simulation. It indi- 
cates that computing transition to turbulence 
in certain spatially evolving boundary layers 
to a reasonable degree of accuracy is possible 
on currently available supercomputers (pro- 
vided there is judicious use of zonal methods 


and high-order-accurate schemes). The finite- 
difference method developed in this study can 
be extended in a straightforward manner to 
curvilinear grids and thus enable direct simu- 
lations of transitional and turbulent flow over 
general geometries. 

Status/Plans. The method has been ex- 
tended to curvilinear grids. A new code to 
compute compressible, turbulent/transitional 
flow over a turbine airfoil is being developed. 

Man Mohan Rai 
Fluid Dynamics Division 
Ames Research Center 
(415) 604-4499 
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INRIA II CASE 3.4 


Mach Contours 


4-2 Multigrid Algorithm for Hypersonic Viscous Flows 


Objective. To develop a multigrid algorithm 
for the efficient solution of viscous flows at 
hypersonic Mach numbers. 

Approach. A Full Approximation Scheme 
(FAS) multigrid was developed based on an 
implicit upwind-biased algorithm for the com- 
pressible Navier-Stokes equations. The algo- 
rithm uses upwind differencing for the con- 
vective and pressure terms and central 
differencing for the viscous terms. The im- 
plicit equations were modified so that the 
primitive variables, instead of the conserved 
variables, were updated during the iteration. 
The time advancement and correction stages 
of the FAS algorithm were modified to main- 
tain positivity of the thermodynamic vari- 
ables, density and pressure. The FAS algo- 
rithm was modified to include multiple itera- 
tions on coarser meshes in order to speed up 
the upstream transfer of information through 
the subsonic part of the boundary layer. 

Accomplishments. The convergence rates 
for a number of hypersonic strong-interaction 
viscous flows were improved with the multi- 
grid scheme. The test cases were those of the 
second Institut National de Recherche en 
Informatique et en Automatique (INRIA) 
Workshop on Hypersonic Flows for Reentry 
Problems. An application to case 3.4, corre- 
sponding to a laminar strong-interaction flow 
over a compression ramp at a Mach number of 
11.68, is shown in the accompanying figure. 
The Mach contours show the leading-edge 
shock and separation-induced shock at the 
compression; the enlarged view shows the 
separated flow in the region of the compres- 
sion corner. The variation of separation ex- 
tent with computer time is shown for the 
multigrid and single grid schemes. An order- 
of-magnitude improvement in efficiency oc- 
curs with the multigrid scheme. 


Significance. The multigrid method, previ- 
ously shown to be highly successful for sub- 
sonic and low supersonic flows, has been ex- 
tended to viscous hypersonic flow. An order- 
of-magnitude increase in efficiency obtained 
with the scheme allowed grid refinement stud- 
ies to be completed for all test cases of the 
INRIA code-validation workshop. 

Status/Plans. Further applications of the 
algorithm are continuing, including blunt 
bodies and three-dimensional configurations. 

James L. Thomas, David H. Rudy 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-2146 
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Figure 4-3* Unstructured and Adaptive Multigrid for the Three-Dimensional 

Euler Equations 
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4-3 Unstructured and Adaptive Multigrid for the Three- 
Dimensional Euler Equations 


Objective. To develop an accurate and effi- 
cient method for computing steady-state com- 
pressible flow about complex three-dimen- 
sional configurations. 

Approach. The steady-state 3-D Euler equa- 
tions are solved on an unstructured tetrahe- 
dral mesh by a Galerkin finite-element 
scheme. The flow variables are stored at the 
vertices of the mesh, and an edge-based data- 
structure is employed to minimize memory 
requirements. An unstructured multigrid 
technique is used to accelerate convergence to 
steady-state. This procedure operates on a 
sequence of non-nested coarse and fine meshes, 
and the patterns for interpolation between 
the various meshes of the sequence are deter- 
mined in a preprocessing step, using an effi- 
cient search algorithm. This strategy is com- 
bined with an adaptive meshing approach, 
where new, finer meshes are automatically 
generated as the solution evolves. 

Accomplishments. The present methodol- 
ogy is capable of providing accurate and effi- 
cient solutions for steady-state inviscid 3-D 
flows without incurring large memory over- 
heads. The figure illustrates the computation 
of transonic flow over an ONERA M6 wing. A 
multigrid sequence of four meshes has been 
employed, with the last two being adaptive 
meshes. The final solution was obtained in 
100 multigrid cycles on the finest grid, which 
required 35 minutes of CPU time on a single 
Cray Y-MP processor, and 22 million words of 
memory. 

Significance. The ability to accurately and 
efficiently predict compressible flows over 
complex 3-D geometries is of particular im- 
portance to the aircraft industry. The use of 
adaptive meshingenables unprecedented reso- 
lution of highly localized phenomena, while 
the multigrid strategy maintains the overall 
efficiency of the scheme. 


Status/Plans. Future work will center on 
the inclusion of more adequate surface model- 
ing and grid generation techniques, as well as 
the extension of the present work to viscous 
turbulent flow cases in three dimensions. This 
will require the implementation of a multiple 
field-equation turbulence model. 

Dimitri J. Mavriplis 

Institute for Computer Applications in 

Science and Engineering (ICASE) 

Langley Research Center 
(804) 864-2213 
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4-4 S3D - An Interactive Surface Grid Generation Tool 


Objective. To develop a surface grid genera- 
tor with geometry modeling and surface 
gridding capabilities in an integrated, inter- 
active environment, with the goal of vastly 
reducing turnaround time. 

Approach. The scope of surface grid genera- 
tion requirements were defined, and areas 
that have not been adequately addressed or 
treated by existing software were identified. 
The initial development effort of S3D has 
been directed at filling or enhancing those 
areas. A prototype S3D with primitive inter- 
face was developed and applied on a variety of 
complex geometries. These applications were 
helpful in identifying problem areas and per- 
mitted a greater focus on such issues as user- 
friendliness and ease of operation. The S3D 
code has since been expanded and enhanced 
with a user interface incorporating the latest 
in workstation technology. The code is now 
ready for hands-on application by CFD ana- 
lysts. 

Accomplishments. S3D can start from a 
geometry definition based either on a 
discretized curve set or a rectangular point 
set. Prominent geometric features such as 
discontinuities and high curvatures are eas- 
ily preserved and deficiencies in geometry 
data easily removed in transforming from 
geometry definition to surface definition. At 
the heart of S3D is the analysis of curves and 
surfaces by robust and widely applicable piece- 
wise cubic and bi-cubic interpolation tech- 
niques made possible by highly reliable curve 
and surface fitting schemes. Point redis- 
tributions are accomplished, through smart 
and user-friendly interfaces, by the use of 
hyperbolic function-based, two-sided stretch- 
ing functions which have been shown to re- 
duce the spacing-induced truncation error. 
Some of the more advanced features of S3D 
include surface-surface intersections, opti- 
mized surface domain decomposition and 
rec'omposition, and automated propagation of 


edge distributions to surrounding patches. In 
the area of database management, S3D per- 
mits a variety of data formats including that 
of PLOT3D for ease of data transfer. 

Significance. The efficiency and ease with 
which S3D can be used to handle commonly 
encountered tasks in surface grid generation 
will help CFD analysts drastically reduce the 
time required in gridding complex surface 
geometries. 

Status/Plans. A workshop on S3D is planned 
for this October. Future plans include adding 
options to handle some of the common data 
exchange formats under Interim Graphics 
Exchange Standard (IGES) and capabilities 
to design and analyze geometries based on 
non-uniform rational B-spline (NURBS). 

Raymond Luh 
Fluid Dynamics Branch 
Ames Research Center 
(415) 604-4494 
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4-5 Accelerated Preparation of Grid Generation Input Data 


Objective. To greatly reduce the user time 
required to generate three-dimensional mul- 
tiple-block computational volume grids. 

Approach. The approach is to use the 
graphical capabilities of a powerful scientific 
workstation to accelerate the process of col- 
lecting and formatting the input data re- 
quired by a grid generation program. This 
new software tool leads the user through the 
process of specifying input data by asking 
questions in sequence, error-trapping the re- 
sponses, and then formatting those data for 
the grid generator. Users having more expe- 
rience can view all the input data in a random- 
access fashion, enter or modify it, and then 
format it with automatic error-trapping. A 
suite of useful utilities completes the software 
system. 

Accomplishments. Preliminary tests indi- 
cate that a reduction by a factor of ten has 
been achieved in the number of man-hours 
required to generate 3-D volume grids. 

Significance. The generation of a suitable 
grid for computational analysis of a realistic 
aerodynamic configuration can take as long 
as 6 months, frequently dwarfing all other 
aspects of the computational fluid dynamics 
(CFD) process combined. This problem be- 
comes even worse when the physical system 
to be modeled undergoes change requiring 
modification to the grid, such as deflecting 
control surfaces, rotating nozzles, or 
aeroelastic effects. Thus, any tool that greatly 
reduces the user time required to generate 
grids is a very significant advance. 

Status/Plans. The software system to effect 
the above, called 3DPREP, is complete. An 
improved program to view 3-D computational 
grids, called 3DECANT, is complete. The 
elliptic volume grid generator, called 
3DGRAPE, is also complete. 

v 


The next step is to combine the three modules 
to produce a truly integrated interactive 3-D 
volume grid generator. But this integration 
task is more significant than one might ex- 
pect. The graphical user interface must run 
on a workstation, but the elliptic grid genera- 
tor must run on a supercomputer. Thus, the 
integrated system must run over a network- 
connecting workstation and supercomputer 
running concurrently. This formidable pro- 
gram design and coding task is now well 
under way. 

The integrated system will further accelerate 
the process by enabling the user to view the 
iterative grid generation process as it is tak- 
ing place. This will give the user a better 
understanding of the grid generation process, 
and the ability to improve a grid generation 
effort that is working or to abort one that is 
not. 

Reese L. Sorenson 
Fluid Dynamics Division 
Ames Research Center 
(415) 604-4471 
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Figure 4-6. Numerical Simulation of the YAV-8B Harrier VSRA in Ground Effect 
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4-6 Numerical Simulation of the YAV-8B Harrier VSRA in 

Ground Effect 


Objective. To develop a validated powered- 
lift flow-field analysis capability useful in the 
design of advanced powered-lift aircraft. 

Approach. The flow field surrounding the 
YAV-8B Harrier in low- velocity jet-borne flight 
has been simulated through the numerical 
solution of the Reynolds-Averaged Navier- 
Stokes equations. This simulation also uses 
an engine model to specify inlet and nozzle 
mass flows and nozzle temperatures based on 
throttle inputs and inlet temperature and 
pressure. Over 2.8 million grid points in 18 
grids are used to resolve the flow field using a 
diagonalized approximately factored algo- 
rithm and the chimera overset grid method. 

Accomplishments. Simulation conditions 
of 30 knots forward speed at 30 feet above 
ground level with exhaust nozzles rotated 81° 
below horizontal have been modeled. This 
computation was performed on a Cray Y-MP 
8/128 using 70 single processor CPU hours 
and running at a rate of 130 MFLOPS (million 
floating point operations per second). Im- 
proved, spatially varying time step methods 
for powered-lift flow fields have allowed a 
factor of six increase in the global time step 
and a corresponding decrease in the cost of 
computation. Software has also been devel- 
oped! for the analysis of infrared visualiza- 
tions which will be used in the validation of 
the aircraft simulation capability. 

Significance. Ground effect flow fields about 
powered-lift aircraft, such as the Harrier, are 
highly complex. Analysis is generally resis- 
tant to small-scale wind tunnel techniques 
due to a wide variety of scaling effects. This 
leads to expensive full-scale powered wind 
tunnel tests to quantify the effects of ground 
proximity on vehicle performance. Some of 
the hazards of operating in the ground envi- 
ronment include hot-gas ingestion, foreign 
object damage, and the “suckdown” effect. 
This simulation demonstrates a promising 


and potentially money saving new technique 
for performance prediction of powered-lift 
designs in the ground environment. 

Status/Plans. As a High-Performance Com- 
puting and Communications Program 
(HPCCP) Computational Aero Sciences (CAS) 
Grand Challenge, the powered-lift project is 
developing numerical analysis tools for com- 
plete powered-lift aircraft on massively par- 
allel computers. Simulation methods for aero- 
dynamics, propulsion, and controls will be 
combined. Additional flight tests are sched- 
uled for the collection of infrared imaging 
validation data. Following validation, it is 
expected that an accurate and neccessary 
design and analysis tool for powered-lift air- 
craft will result. 

William R. Van Dalsem, 

Fluid Dynamics Division 
Ames Research Center 
(415) 604-4469 
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Figure 4-7. High-Speed Civil Transport Navier-Stokes Configuration 
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4-7 High-Speed Civil Transport Navier-Stokes Computations 


Objective. To develop an efficient numeri- 
cal procedure for computing viscous flow over 
a high-speed civil transport configuration. 

Approach. A multistage Runge-Kutta time- 
stepping scheme with multigrid acceleration 
technique was employed for computing the 
steady-state solutions of the thin-layer Navier- 
Stokes equations. The accuracy of the nu- 
merical solutions was enhanced by the use of 
a matrix dissipation model. 

Accomplishments. A multigrid-based fi- 
nite volume numerical scheme developed for 
computing high Reynolds number viscous 
flows over aircraft components at transonic 
speeds has been extended to accommodate 
supersonic flows. The resulting numerical 
code, known as TLNS3D, has been applied to 
compute viscous flow over a generic high- 
speed civil transport model at the design 
cruise Mach number M = 3. The angle-of- 
attack varied from 0° to 8° in these computa- 
tions to correspond to an in-house experimen- 
tal investigation, and each test condition re- 
quired approximately 2.5 hours on a Cray-2 
computer. The computed lift and drag coeffi- 
cients are compared with the experimental 
data in the adjoining figure. The computed 
solutions are in excellent agreement with the 
data over the entire angle-of-attack range, 
even for the drag coefficient, which is domi- 
nated by the viscous component at lower lift 
values. 

Significance. Navier-Stokes solutions for 
the flow over configurations of practical inter- 
est can be obtained routinely using the multi- 
grid-based TLNS3D code. Because of the 
improved efficiency of this code, parametric 
studies to examine the performance of differ- 
ent configurations are feasible without re- 
sorting to wind tunnel testing in early phases 
of vehicle development. 


Status/Plans. The method described here 
will be generalized to accommodate block- 
structured grids so that flow over complex 
aircraft can be calculated. 

Veer N. Vatsa 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-2236 
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Figure 4-8. CFD-Based Hypersonic Wind Tunnel Nozzle Design 
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4-8 CFD-Based Aerodynamics Design of Hypersonic Wind Tunnel 

Nozzles 


Objective. There was an immediate need for 
a new design procedure for high Mach num- 
ber hypersonic wind-tunnel nozzles, for which 
the classical design procedures (using method- 
of-characteristics wi th boundary layer correc- 
tion) have failed to produce the high-quality 
flow fields necessary for fundamental flow 
experiments and computational fluid dynam- 
ics (CFD) validation studies. The classical 
design procedure for hypersonic nozzles has 
been shown to break down at Mach numbers 
greater than 8 and when the boundary layer 
covers a large percentage of the test section 
diameter. 

Approach. To develop a procedure that uni- 
fies the best of classical design practices, 
CFD, and optimization procedures, an effi- 
cient CFD code, which solves the parabolized 
Navier-Stokes (PNS) equations using an ex- 
plicit upwind algorithm, was coupled to a 
nonlinear least-squares (LS) optimization 
procedure. A LS problem was formulated to 
minimize the difference between the com- 
puted flow field and the objective function, 
and consisted of the centerline Mach number 
distribution and the exit Mach number and 
crossflow velocity profiles. The aerodynamic 
lines of the nozzle were defined using cubic 
splines, the slopes of which are optimized 
with the design procedure. The thick bound- 
ary layer typical of hypersonic nozzles is ac- 
counted for by incorporating the solution of 
the PNS equations in the design procedure. 

Accomplishments. Initial programming and 
validation have been completed on a com- 
puter code that exercises a LS/PNS optimiza- 
tion procedure for designing hypersonic wind 
tunnel nozzles. The program couples an opti- 
mum design procedure based on the LS mini- 
mization problem to an explicit upwind PNS 
CFD code. In the accompanying figure, a 
Mach 15 helium nozzle designed using the 
new CFD-based procedure is compared with a 
nozzle designed with the method-of-charac- 


teristics with boundary layer (MOC/BL) cor- 
rection methods. The CFD-based nozzle de- 
sign reduced the fluctuations of static pres- 
sure to less than + 2% for the uniform core 
region. This is a significant improvement 
over the classically designed nozzle, which 
had fluctuations of static pressure greater 
than ±30% for the core region. 

Significance. The new procedure for design- 
ing hypersonic wind-tunnel nozzle contours 
demonstrates that a CFD computer code can 
be efficiently coupled to an optimization algo- 
rithm. The effect of changes in the nozzle wall 
parameters is evaluated by computing the 
nozzle flow field using the PNS equations. 
This new procedure provides a method for 
designing hypersonic nozzles of high Mach 
numbers in which classical procedures break 
down. 

Status/Plans. Work is under way to apply 
the new procedure to design other wind-tun- 
nel nozzles and to extend it for use in the 
design of hypersonic inlets. 

John J. Korte 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-6920 
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4-9 Turbulent Flow Over a Backward Facing Step 


Objective. To develop a detailed database 
of the turbulent flow over a backward facing 
step for the validation of direct simulations, 
large-eddy simulations, and turbulence mod- 
els. 

Approach. State-of-the-art instrumenta- 
tion was used to acquire the mean velocities 
and turbulence quantities in the turbulent 
flow. The primary instrumentation was a 
three-component laser velocimeter, which al- 
lowed the measurement of all three velocity 
components and the full Reynolds stress ten- 
sor. Test conditions and geometries were 
defined through interaction with the compu- 
tational fluid dynamics (CFD) developers. 

Accomplishments. During FY 1990, the 
initial investigation of the flow over the step 
was performed. The test conditions for that 
experiment were chosen to be suitable for 
direct simulations. The detail and quality of 
the data generated interest in the turbulence 
modeling community, which requested a 
higher Reynolds number case. A new facility 
has been constructed that allows a wide range 
of Reynolds numbers and expansion ratios. 
Currently, the flow at a Reynolds number of 
47,000 based on step height with an expan- 
sion ratio of 1:2 is being investigated. 


Significance. The data being acquired in 
this investigation are unique in that they 
document the full Reynolds stress tensor. 
Upon completion of the experiment, this data- 
base will represent the most complete statis- 
tical database for the validation of CFD meth- 
ods. 

Status/Plans. The acquisition of the flow- 
field data is currently under way and will 
continue throughout the next fiscal year. 

Scott O. Kjelgaard 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-2160 
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4-10 Near-Wall Laser Velocimeter Measurements 


Objective. To increase near-wall Reynolds 
stress database for three-dimensional turbu- 
lent boundary layer flows with significant 
turning and separation. 

Approach. A new Laser Doppler Velocime- 
ter (LDV) instrument was developed capable 
of obtaining near-wall measurements in a 
three-dimensional boundary layer. Measure- 
ments were performed in a 3-D turbulent 
boundary layer flow apparatus, located at 
Stanford University for which previous mea- 
surements of the Reynolds stresses were ob- 
tained with a hot-wire anemometer. 

Accomplishments. Turbulence measure- 
ments were acquired in a 3-D boundary layer 
at a distance of 0.02 mm (0.001 in.) from the 
surface. Conventional LDV approaches are 
limited to a distance of about 1 mm from the 
surface. The technique has undergone fur- 
ther improvements, such as the development 
of a single channel multifrequency signal pro- 
cessing method for measuring the cross corre- 
lations. 3-D Navier-Stokes calculations (on a 
coarse grid) have been performed for the 3-D 
wedge flow experiment. 

Significance. The accurate prediction of 
skin friction and heat transfer is critically 
dependent on the modeling of the Reynolds 
stresses within the viscous sublayer. The 
present LDV development will allow the mea- 
surement of these stresses within the 
sublayers even under highly 3-D conditions. 
Such data will aid in the development and 
validation of improved turbulence closure 
models. 


Status/Plans. Further develop the multifre- 
quency signal processing system to obtain 
velocity components -vw and -uv Reynolds 
shear stresses. Fiber-optics are being devel- 
oped to facilitate 3-D measurements. 

Dennis A. Johnson 
Fluid Dynamics Division 
Ames Research Center 
(415) 604-5399 
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Chapter 5 

Numerical Aerodynamics Simulation (NAS) 


The Numerical Aerodynamic Simulation (NAS) program provides readily accessible 
supercomputing capability to the United States’ top aeronautical researchers in Government, 
industry, and academia. The NAS program also includes research and technology development 
to ensure application of emerging technologies to computational fluid dynamic and other 
computational sciences. Specifically, the current research involves the enhancement of user 
interfaces and software and hardware technology for parallel computer architectures. 

The objectives of NAS are (1) to maintain a pathfinding role in providing leading-edge 
supercomputing capabilities to NASA, DoD, and other Government agencies, industry, and 
universities as a critical element for continued leadership in computational aeronautics and 
related fields; (2) to stimulate the development of state-of-the-art, large-scale, computer systems 
and advanced computational tools for pioneering research and development; and (3) to provide 
a strong research tool for OAST. 

To maintain the lead in large-scale computing capability, NAS is implementing a strategy of 
installing, at the earliest possible opportunity, the most powerful high-speed processor (HSP) 
available. NAS maintains at least two HSPs, one of which is fully operational and represents 
more mature technology, and the other, which is a higher performance prototype or early 
production model. The current HSPs are a Cray-2 installed in January 1988 and a Cray Y-MP 
installed in November 1988. The Cray Y-MP is the first computer to sustain a computation rate 
of a billion floating point operations per second (GFLOP). Current plans are to replace the Cray- 
2 with a new processor (HSP-3) in 1992. 

NAS was the first supercomputing facility to install a standard operating system (UNIX) and 
communication software on all systems. UNIX offers the flexibility of both batch and interactive 
computing and provides a common user interface on all user-visible subsystems. NAS is 
currently upgrading its remote network nationwide from NASnet, a switch-based, high-perfor- 
mance communication network, to Aeronet, a router-based system, allowing researchers at 
remote locations to have similar interactive capability as local users at the NAS facility. 

The vision for the NAS program is to provide by the year 2000 an operational computing system 
capable of simulating an entire aerospace vehicle system within a computing time range fttwn one 
to several hours. It is estimated that a computing time rate of one trillion floating point 
operations per second (TFLOPS) is required to accomplish this goal. 

Program Manager: Pamela F. Richardson 

OAST/RF 

Washington, DC 20546 

(202) 453-9857 
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5-1 Parallel Computers in the NAS Program 


Objective. To evaluate the performance 
capabilities of parallel computers running 
computational fluid dynamics (CFD) codes. 

Approach. Abase level of performance sta- 
tistics was established and algorithms that 
increase the performance of parallel comput- 
ers were developed. 

Accomplishments. The NAS Program has 
acquired two parallel computers, a Connec- 
tion Machine-2 (CM-2) from Thinking Ma- 
chines, Inc. (TMI) and an iPSC/860 system 
from Intel. The CM-2 has 32,000 1-bit serial 
processors and 1000 64-bit Weitek floating 
point processors with a SIMD (single instruc- 
tion, multiple data) architecture. The peak 
speed is 14 GFLOPS (giga [billion] floating 
point operations per second). Scalar code is 
executed on either Sun or VAX front-end 
computers, which broadcast instructions to 
the CM-2 processors. Interprocessor commu- 
nication occurs three ways: general router, 
compiled, and NEWS (North, East, West, 
South, that is, nearest neighbor). The CM-2 
has 4 gigabytes of memory, and a 25 gigabyte 
Datavault for on-line storage as well as a 
high-performance parallel interface (HiPPI) 
to UltraNet. The CM-2 was upgraded in 
January 1991 to take advantage of TMI’s new 
slice- wise architecture, which presents a pro- 
gramming model of the CM-2 based on 1000 
floating point processors, instead of the 32,000 
bit serial processors. Such a model is more 
relevant to the CFD work done at NAS. TMI 
has recently provided a time-sharing operat- 
ing system and network queing system. 

The iPSC/860 system was installed in Janu- 
ary 1990 as part of Intel’s gamma prototype 
project with DARPA. This system has 128 
nodes, each with an Intel i860 64-bit micro- 
processor and 8 megabytes of memory. It has 
a 40 megahertz RISC chip rated at 60 
MFLOPS (M=mega [million]). Peak perfor- 
mance is thus 7.5 GFLOPS. Total memory 


available is 1 gigabyte. The iPSC/860 is a 
MIMD (multiple instruction, multiple data) 
architecture. Interprocessor communication 
is through message passing. The intercon- 
nect network is a hypercube. The iPSC/860 
has a concurrent file system of 10 gigabytes 
supporting a peak transfer rate of 10 mega- 
bytes/second. At present, all communication 
is done through a small (PC 80386) front-end 
computer, which also handles partitioning of 
the nodes and loading of jobs for execution. 

Significance. The highest performance at 
NAS for CFD-related applications on parallel 
computers is 1.6 GFLOPS using a 32-bit iso- 
tropic turbulence code. This code used an 
assembly language FFT (Fast Fourier Trans- 
form) and a special compiler, Vectoral. The 
problem size was 256^ nodes. A similar code 
runs at about 160 MFLOPS on a single Cray 
Y-MP processor. This code is achieving a 10- 
fold performance-to-price increase on the Intel 
compared to the Cray Y-MP. 

Status/Plans. NAS will continue to explore 
the potentials of both parallel systems. 

T. A. Lasinski 
NAS Systems Division 
Ames Research Center 
(415) 604-4405 
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Performance of NAS Pseudo CFD Applications (7/24/91) 


Algorithm 
64 3 grid 

Computer 

Number of 
Processors 

Time/Iter 

(sec) 

MFLOPS 

(YMP) 

Efficiency 

(actual/peak) 

LU 

BBN 

62 

12.13 

36 

1/17 


CM2 

32K 

2.84 

150 

1/86 


iPSC 

128 

2.34 

182 

1/42 


YMP 

1 

2.10 

203 

3/5 


YMP 

8 

0.35 

1218 

<1/2 

SP 

BBN 

62 

3.71 

95 

1/7 


CM2 

32K 

11.63 

25 

1/521 


iPSC 

64 

2.42 

122 

1/31 


YMP 

1 

1.49 

198 

>1/2 


YMP 

8 

0.26 

1335 

1/2 

BT 

BBN 

62 

11.83 

77 

1/8 


CM2 

32K 

24.25 

37 

1/352 


iPSC 

64 

4.54 

199 

1/19 


YMP 

1 

4.22 

214 

3/5 


YMP 

8 

0.64 

1411 

>1/2 


LU: Lower Upper Triangular Solve 
SP: Scalar Pentadiagonal 
BT: Block Tridiagonal 


Figure 5-2. Performance of NAS Pseudo CFD Applications 
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5-2 NAS Parallel Benchmarks 


Objective. To develop a new methodology for 
the performance evaluation of highly parallel 
computers running computational fluid dy- 
namics (CFD) applications. 

Approach. The principal distinguishing fea- 
ture of the new benchmarks is their “pencil 
and paper'’ specification — all details of these 
benchmarks are specified only algorithmically. 
In this way many of the difficulties associated 
with conventional benchmarking approaches 
on highly parallel systems are avoided. The 
benchmarks consist of a set of kernels, the 
“parallel kernels,” and a simulated applica- 
tion benchmark. Together they mimic the 
computation and data movement character- 
istics of large-scale CFD applications. 

Accomplishments. A set of benchmark ker- 
nels and pseudo applications has been se- 
lected and completely specified. Sample se- 
rial implementations for the benchmarks have 
been developed. Parallel versions of the bench- 
marks have been implemented on the 8-pro- 
cessor Cray Y-MP, the 128-processor Intel 
iPSC/860, and the 32,000-processor Connec- 
tion Machine CM-2 (from Thinking Machines). 
The benchmarks have been distributed to a 
large number of computer vendors and aca- 
demic researchers interested in performance 
evaluation. 

Significance. The NAS Parallel Benchmarks 
represent the first comprehensive effort to 
assess the performance of highly parallel com- 
puters. They present a new methodology that 
allows the suitability of parallel supercom- 
puters for aerosciences applications to be 
evaluated. In addition, the benchmarks are 
an invaluable tool for selecting computers for 
the High Performance Computing and Com- 
munication Program (HPCCP) testbeds and 
will ensure the selection of computers suit- 
able for the HPCCP “Grand Challenge” prob- 
lems. 


Status/Plans. Implementation on other par- 
allel computers will continue in the next year. 
The inclusion of additional kemels/applica- 
tions into future releases of the benchmarks 
is under investigation. 

Horst D. Simon 

Computer Sciences Corporation 
NAS Systems Division 
Ames Research Center 
(415) 604-4322 
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Figure 5-3. Virtual Wind Tunnel 
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5-3 Virtual Wind Tunnel 


Objective. To explore the effectiveness of 
virtual environment technology for flow visu- 
alization. 

Approach. Virtual environment technology 
provides a new approach to user interfaces in 
computer software. This approach involved 
integrating a variety of input and display 
devices to give the user the illusion of being 
immersed in an interactive computer-gener- 
ated environment. The computer-generated 
scene was displayed in stereo creating the 
illusion of depth, and was rendered from a 
point of view that tracks the user’s head 
movements. The user had an input device, 
typically an instrumented glove, through 
which objects appeared to be directly manipu- 
lated in the computer-generated environment. 

Accomplishments. A virtual environment 
was implemented for exploring numerically 
generated three-dimensional unsteady flow 
fields. A variety of interesting techniques 
were used for visualization and for navigation 
through the flow. A boom-mounted six de- 
gree-of-freedom head position-sensitive ste- 
reo CRT system was used for viewing. AVPL 
Dataglove™ Model II, which incorporates a 
Polhemus 3 Space™ tracker, was used to 
sense the user’s hand position, orientation, 
and finger joint angles. The user’s gestural 
movements were interpreted as commands 
for injecting various tracers (e.g., “smoke”) 
into the virtual flow field. An eight-processor 
Silicon Graphics computer system was used 
for computation and rendering. 

Significance. We have shown that a high- 
performance graphics workstation can be used 
to visualize three-dimensional flow fields in a 
virtual environment at acceptable frame rates; 
virtual environments are an effective tool for 
the rapid exploration of three-dimensional 
flow fields, both steady and unsteady; and the 
system may be applicable to the visualization 
of other vector fields. 


Status/Plans. Future plans include distrib- 
uting the software so that computation and 
data management is done on a supercomputer 
and upgrading to a higher resolution head- 
tracking display system. 

Creon Levit 
NAS Systems Division 
Ames Research Center 
(415) 604-4403 
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5-4 Vector Field Topology Visualization Software 


Objective. To create simple visualizations 
that capture the important properties of vec- 
tor fields. 

Approach. A new module, called Topology, 
has been created for the Flow Analysis Soft- 
ware Toolkit (FAST) graphics software pack- 
age to find and display a portion of the topol- 
ogy of vector fields. FAST is a multiprocess 
software environment that consists of a col- 
lection of separate modules that can be run 
simultaneously, sharing data and allowing 
the user to create, load, analyze, visualize, 
animate, and record computational fluid dy- 
namics (CFD) data. FAST was created at 
Ames Research Center, and beta versions of it 
are currently being tested at various U.S. 
sites. 

Accomplishments. Topology, the new FAST 
graphics software module, finds, classifies, 
and displays critical points, i.e., where a vec- 
tor field vanishes. The software integrates 
curves through the field from initial positions 
near critical points on the invariant mani- 
folds of the linearized field near the points. By 
carefully choosing the curves to integrate, the 
user can find the vortex cores of all vortices 
containing critical points and/or examine the 
skin friction surface topology to determine 
the location of separation lines and other 
features of interest. The software has very 
flexible mechanisms for exploring other topo- 
logical features as well. 

Significance. Users of Topology can quickly 
find and display some of the important fea- 
tures of vector fields. In particular, separa- 
tion lines and most vortex cores may be found 
quickly, easily, and, in many cases, automati- 
cally. Since the software functions properly 
on the large, multi-zoned, blanked (masked) 
grids typical of state-of-the-art CFD, scien- 
tists can gain additional insight into the re- 
sults of their computations. 


Status/Plans. We plan to extend the soft- 
ware to improve visualization of vortices and 
to integrate topological surfaces from the 
neighborhood of critical points. We will inves- 
tigate implementing Topology in other visual- 
ization environments, particularly the FLORA 
(FLOw Real-Time Analysis) software, also 
developed at Ames, which enables a user to 
interactively generate particle traces and sur- 
faces, and EXPLORER, developed by Silicon 
Graphics, Inc. 

A1 Globus 

Computer Sciences Corporation 
NAS Systems Division 
Ames Research Center 
( 415 ) 604-4404 
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Figure 5-5. Distributed and Cooperative Visualization of Unsteady Fluid Flow 
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5-5 Distributed and Cooperative Processing of Unsteady 

Fluid Flows 


Objective. To utilize distributed processing 
and cooperative processing to provide a pow- 
erful tool for visualizing computational fluid 
dynamics (CFD) data. 

Approach. Distributed processing provides 
the user with the combined capabilities of a 
graphics workstation and a supercomputer, 
under a single application. The supercom- 
puter's capabilities ofhigh performance, large 
memory, large disk storage, and fast disk 
access are necessary to contain, access, and 
calculate the large data sets endemic to un- 
steady flow analyses. The high-speed graph- 
ics capability of the workstation presents high- 
resolution animations of features extracted 
from large data sets. In addition, the human 
interface environment provided by the work- 
station allows for controlling the analytical 
tools of the visualization system. 

With cooperative processing, visualization is 
used as a communication medium. The para- 
digm is illustrated by the analogy of one 
scientist outlining ideas on a chalkboard for 
another scientist. Equal accessibility to the 
chalkboard by both scientists provides a plat- 
form for the exchange of ideas. The images 
produced on graphics workstations can be the 
platform. In combination with distributed 
processing, scientists may use this platform 
at geographically distant locations. 

Accomplishments. Tempus Fugit/Interview 
is a CFD visualization application. Tempus 
Fugit (“time flies”) interactively creates im- 
ages animated over time from large data sets 
representing unsteady fluid flows. The com- 
panion program, Interview, presents simul- 
taneously, identical and animated images to 
another workstation in a shared environment 
with Tempus Fugit. Interview attaches to the 
supercomputer process created by Tempus 
Fugit and transforms the process to a shared 
environment. Identical images are presented 
to both workstations simultaneously. Thus, a 


cooperative environment is created via a net- 
work connection to the supercomputer and is 
uninhibited by geographical proximity. 

Significance. The initial implementation of 
Tempus Fugit/Interview is a proof-of-the-con- 
cept capability providing interactive visual- 
ization of unsteady fluid flows by using dis- 
tributed processing. Interview’s prototype 
WYSIWIS (“what you see is what I see”) 
interface provides a base for further research 
into group communication via a visualization 
medium. 

Status/Plans. A prototype of Tempus Fugit/ 
Interview has been completed. Additional 
visualization techniques and cooperative in- 
teraction methods will be added as opera- 
tional experience is gained. 

Michael J. Gerald-Yamasaki 
NAS Systems Division 
Ames Research Center 
(415) 604-4412 
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Chapter 6 
Rotorcraft 


The objective of the Rotorcraft Program is to provide the technologies for helicopters and other 
rotorbome aircraft to achieve quiet, low-vibration operation with increased performance, agility, 
maneuverability, and stability, all with acceptable handling qualities. Much of the work is done 
in conjunction with the U.S. Army and the FAA and in cooperative programs with industry. 

Improved analysis can now handle many local aerodynamic phenomena, but the transonic, 
unsteady, complex wake interaction flow requires another generation of codes. Therefore, much 
of the rotorcraft program still is empirically based, with validation of analysis as a goal. This 
approach requires new test techniques, upgraded test facilities, carefully instrumented models, 
sophisticated simulations, and increased use of the latest generation of computers. New ideas 
also are part of the program, since the complexity of rotorcraft makes it a fertile area for 
innovative approaches. 

The heart of the aerodynamic portion of the program is airloads research. Small scale, pressure- 
tapped blade data have become available with the full-scale results due in 1992. Complementary 
efforts are under way in component, interference, and wake testing and prediction. These 
databases are used in acoustic and vibration research, which constitutes more than half of the 
program resources. The airload models also will be used for rotor state control and higher 
harmonic control to suppress noise or vibration or to enhance maneuverability. 

Flight dynamics research in the Aerodynamics Division is based on unsteady rotor dynamics, 
with emphasis on simulation and eventually flight tests using a variable stability UH-60. New 
flight test capability is planned, with research challenges in higher-frequency control for rotor 
state control to achieve high agility. 

Higher-speed rotorcraft activities are now focused on commercial tiltrotor technology. Certifi- 
cation issues are being addressed on the simulator and in noise testing and prediction. Also under 
way are improvements in performance, interior noise, vibration, and stability. The importance 
of these research opportunities for civil applications have been reported under contract to Boeing 
Commercial Airplane Company. 

Program Manager: Mr. George Unger 

OAST/RF 

Washington, DC 20546 
(202) 453-5420 
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Figure 6-1. Large-Scale Tiltrotor Performance Program 
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6-1 Large-Scale Tiltrotor Performance Program 


Objective. To acquire rotor performance 
and loads data with and without the wing 
installed to determine the influence of the 
wing on rotor operation, provide rotor perfor- 
mance data up to and beyond the aircraft 
operating envelope, and evaluate acoustics 
for terminal area operations. 

Approach. A large-scale test program was 
developed for the 40- by 80 foot wind tunnel. 
The Prop Test Rig (PTR) was used to evaluate 
the forward flight performance for three dif- 
ferent large-scale rotor blade sets including 
the V-22, XV-15 Metal, and XV-15 ATB (Ad- 
vanced Technology Blades). The required 
measurements include acoustic as well as 
dynamic and performance data. Measure- 
ments made with and without the wing would 
be used to quantify effects of the wing on the 
rotor performance and acoustics. 

Accomplishments. Testing was initiated, 
and isolated rotor data up to 240 knots was 
acquired on the 0.658-scale V-22 rotor prior to 
incident. The PTR recovery plan was initi- 
ated to rebuild/strengthen hardware to pro- 
vide full test envelopes in both 40 by 80 and 80 
by 120 test configurations. 

Significance. The data obtained during 
testing will be used to provide experimental 
validation of the V-22 tiltrotor aircraft cruise 
performance. Future acoustic data will pro- 
vide critical information into the feasibility of 
tiltrotor terminal air operations. 


Status/Plans. The Ames Prop Test Rig will 
be rebuilt for both 40 by 80 cruise perfor- 
mance testing and 80 by 120 transition test- 
ing for performance and acoustics. The test- 
ing will use existing government-owned 25 ft 
diameter XV-15 rotor hardware. 

Jeffrey Light 

Rotorcraft Aeromechanics Branch 
Ames Research Center 
(415) 604-4881 
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SYMMETRIC FLIGHT MODE 
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Figure 6-2. Ultrotor Aeroelastic Stability Control 
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6-2 Tiltrotor Aeroelastic Stability Control 


Objective. To investigate the use of active 
controls to delay tiltrotor whirl-flutter type 
instability. 

Approach. The comprehensive rotorcraft 
analysis code CAMRAD/JA was used to ob- 
tain a set of linear differential equations, 
which describe the motion of a tiltrotor air- 
craft at various speeds. The hub motions from 
wing/body vibration are a standard input to 
CAMRAD/JA. The hub-motion is calculated 
using a separate structural analysis program. 
The CAMRAD/JA output consists of the open- 
loop system matrices which describe the air- 
craft motion in the state-variable domain. 
The matrices formed an input to a separate 
program, which performed the closed-loop, 
active control calculations. Additional input 
consists of the sensor model and the active 
controls feedback gain factors. The program 
performed an eigen-value analysis to deter- 
mine the flutter stability for both the open- 
and closed-loop systems. Time response cal- 
culation can be performed to estimate the 
magnitude of the required active control in- 
put for closed-loop stability. A separate util- 
ity was used to combine aircraft trim informa- 
tion from CAMRAD/JA with wing vibration 
data from the structural analysis program to 
define a sensor model, which is based upon 
physically measurable signals. 

Accomplishments. The use of active con- 
trols to alleviate whirl-flutter was investi- 
gated for a cantilevered, high-wing tiltrotor 
model and an advanced joined- wing tiltrotor 
aircraft. Sensor models, defining the feed- 
back of either pure state-variables or physi- 
cally measurable wing accelerations to the 
rotor cyclic pitch were evaluated. The analy- 
sis showed that for both tiltrotor aircraft con- 
figurations, the feedback of accelerations, rep- 
resenting the wing torsion and chordwise 
motions, to the longitudinal cyclic pitch could 
increase the flutter velocity for the symmetric 
flight mode. The figure shows that for the 


cantilevered-wing tiltrotor configuration the 
wing-chord instability could be increased from 
280 knots to above 300 knots by feeding back 
vertical and horizontal accelerations mea- 
sured near the wing tip to the longitudinal 
cyclic pitch. 

Significance. The results of the investiga- 
tion have shown that (1) active controls, using 
the rotor cyclic pitch can be used to delay the 
occurrence of whirl flutter instability; (2) an 
active cyclic pitch control level of 0. 010°-0.012°, 
which equates to a 10-15 pound active control 
force applied at the hub, was required to 
stabilize the tiltrotor; and (3) contamination 
of the sensor output signal with a 25% noise 
level did not adversely affect the closed-loop 
stability. 

Status/Plans. The whirl-flutter alleviation 
investigation will be extended to include the 
anti-symmetric flight mode and for the com- 
plete aircraft in which the symmetric and 
anti-symmetric flight modes are coupled 
through the active control system. Optimiza- 
tion of the feedback system for a given sensor 
model will also be studied. 

Johannes M. van Aken 
Rotorcraft Aeromechanics Branch 
Ames Research Center 
(415) 604-6668 
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6-3 Tiltrotor Download Reduction 


Objective. To understand the physics of 
tiltrotor download and develop download re- 
duction concepts to minimize the adverse per- 
formance impact of the tiltrotor’s wake on the 
aircraft’s wing in hovering, vertical, and tran- 
sition flight. 

Approach. An experimental and analytical 
program has been pursued to quantify and 
demonstrate download reduction concepts. A 
model 7-foot diameter V-22 rotor system was 
used in hover at the Outdoor Aerodynamic 
Research Facility (OARF). Full potential and 
Navier Stokes models of tiltrotor hovering 
flow fields have been successfully developed 
to characterize and analyze the physics of the 
unsteady, vortical flow fields. 

Accomplishments. A successful hover test 
was conducted at the OARF in FY 1991. 
Results included a demonstration of down- 
load reduction with upper surface wing blow- 
ing. Shadowgraphs of the wake/wing interac- 
tion were also acquired. Several papers were 
presented on three tiltrotor wake analysis 
codes. Two steady codes, a full-span potential 
flow code and a semi-span Navier-Stokes code, 
were developed. Wing pressures within the 
rotor’s wake as well as outside the wake along 
the wing were predicted and compared. A 
new two-dimensional vortical wake interac- 
tion code using a velodty/vorticity formula- 
tion on a staggered grid was reported. 

Significance. Results from the experimen- 
tal program are directly applicable to aircraft 
configurations and result in a 10 to 20% re- 
duction in download on current generation 
tiltrotor designs. The shadowgraph wake 
flows quantify the unsteady flow field, and 
the wing pressure database will be used to 
validate new 3-D analysis codes. The compu- 
tational methods are the first successful analy- 
ses of tiltrotor flow fields in hover, and are 
available to industry for wing design on the 
next generation of tiltrotor aircraft. 


The 2-D velocity/vorticity formulation results 
are the first-ever calculations to accurately 
predict the distributed airfoil pressures on 
both the upper and lower wing surfaces. The 
staggered grid method provides for the direct 
solution of the incompressible Navier-Stokes 
equations by means of a fully coupled implicit 
technique. This unique method was demon- 
strated to accurately predict the mean lift and 
drag and the periodic loading for several 2-D 
bodies including an NACA 0012 airfoil with 
and without a deflected flap at -90° incidence. 

Status/Plans. Plans are to report the results 
of the experimental download program and to 
conduct a second test evaluation of other 
download reduction concepts. The prediction 
methodology will be expanded to include 3-D 
steady and unsteady analysis , turbulent flow, 
and rotor/fuselage interaction. 

Paul M. Stremel 

Rotorcraft Aeromechanics Branch 
Ames Research Center 
(415) 604-4563 
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Figure 6-4. Directivity Characteristics of the XV- 15 in Hover 




6-4 Tiltrotor Hover Acoustics 


Objective. To perform full-scale flight tests 
for validation of ROTONET and to determine 
methods of minimizing rotorcraft noise. 

Approach. A cooperative program with Ames 
Research Center has been defined to evaluate 
the acoustic impact in hover of the Advanced 
Technology Blade system on the XV-15 
tiltrotor aircraft. A database for validation of 
ROTONET is an objective of this program. 

Accomplishments. A hover test evaluation 
of the XV-15 tiltrotor aircraft was conducted 
at Ames Research Center. The hover test 
program was conducted in December 1990. 

Significance. The results of the hover test 
program have shown that: (1) A significant 
reduction of in-plane noise was achieved by 
reduced rotational speed, as expected. (2) The 
longitudinal directivity measurements iden- 
tified a strong impulsive character propagat- 
ing aft of the vehicle. This character has been 
identified before and related to the “fountain 
effect.” Triangulation is currently in progress 
to identify the location of the impulsive source 
(see figure). 


Status/Plans. The acoustic data from the 
hover test program is being prepared as two 
articles for publication. AHS Technical Spe- 
cialist Meeting Rotorcraft Acoustics and 
Fluid Dynamics. Follow-on test programs are 
planned for level flight and terminal area 
descent/ascent conditions: a quick program 
scheduled for August 1991 and a more com- 
prehensive program scheduled for early 1993. 


Danny R. Hoad 
Applied Acoustics Branch 
Langley Research Center 
(804) 864-5060 

Brent Wellman, Martin Maisel 
Flight Experiments Branch 
Ames Research Center 
(415)604-6573/6372 
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Figure 6-5. Civil Tiltrotor Technology 
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6-5 Civil Tiltrotor Technology 


Objective. To determine the technical, eco- 
nomic, and environmental factors that are 
critical to the success of commercial, passen- 
ger-carrying tiltrotor aircraft in the national 
transportation system. 

Approach. A follow-on study was performed 
with Boeing Commercial Airplane Group 
teamed with Bell Helicopter Textron and 
Boeing Helicopter to identify highest payoff 
technologies, economic evaluations, market 
size and “cost/price/ market” loop closure, 
civil requirements, and the potential changes 
needed to the V-22 for civil application. Op- 
erational analyses were conducted for steep 
approaches to vertiports, and actions required 
to achieve a national civil tiltrotor transpor- 
tation system were defined. 

Accomplishments. The study found antici- 
pated demand for a fully designed commercial 
passenger carrying tiltrotor is forecast to be 
2625 vehicles by the year 2000, growing to 
4925 in the year 2010. The market value of 
2625 aircraft is $32 to $42 billion, with half 
this being a domestic market and the remain- 
der a potential worldwide export market. The 
study found that the tiltrotor has promise as 
a means to relieve airport congestion signifi- 
cantly. A projected airline schedule of airline 
service in the northeast corridor indicated 
that as replacement commuter aircraft a 
tiltrotor network could open up 1000 airport 
slots per day (1/3 of 1989 daily slots). 

The study found that the cost/price/market 
loop can be closed. A fully designed commer- 
cial passenger-carrying tiltrotor was found to 
be different from the V-22 and must overcome 
three barrier technology challenges: (1) low 
external noise, (2) human factors-based pilot 
controls for commercial flight, and (3) termi- 
nal area low-speed control and navigation, 
allowing safe, precise, steep approaches. The 
study recommends increased research to over- 

v 


come the technology barriers and to develop 
enhancing technology to encourage civil ap- 
plications. 

Significance. The United States holds a 
commanding lead in tiltrotor technology. This 
study identified the opportunities, challenges, 
and actions required to develop a civil tiltrotor 
aircraft that is both market responsive and 
effectively integrated in the national trans- 
portation system utilizing its unique capa- 
bility. The development of the recommended 
actions will reduce the risk to industry. 

Status/Plans. The Phase II Study is com- 
plete and a final report was released in Febru- 
ary. NASA is working with FAA to help 
determine Pre-Terminal Operating Procedure 
Criteria for safe, steep approaches to urban 
vertiports. Barrier technologies and other 
high payoff technologies are being assessed 
by NASA for future R & T consideration and 
the Advanced Tiltrotor Transport Technology 
new initiative. 

John Zuk 

Civil Technology Office 
Ames Research Center 
(415) 604-6568 
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XV-15 APPROACH 



Figure 6-6. Measured Contours of the 65 Day-Night Noise Level (DNL) for the XV-15 
Tiltrotor Aircraft With 3 Blades, Corrected for 4 Blades, for Takeoff and Approach Conditions 
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6-6 Theoretical Determination of Noise Reduction With 
Increasing Blade Number for the XV-15 Tiltrotor Aircraft 


Objective. To predict the change in the 65 
day-night noise level (DNL) contour area for a 
tiltrotor type aircraft as a result of increasing 
the number of rotor blades per rotor. The 
flight conditions of particular interest were 
takeoff and approach. 

Approach. Contour corrections were ob- 
tained theoretically by comparing ROTONET 
acoustic predictions made for the baseline 
XV-15 tiltrotor system (three blades per ro- 
tor) with predictions made for the XV-15 modi- 
fied to four blades per rotor. The resultant 
correction factor was applied to modify a mea- 
sured XV-15 ground contour. 

The predictions were calculated using 
ROTONET Phase I and included propagated 
tone and broadband noise. Effective per- 
ceived noise level (EPNL) values were calcu- 
lated at five rows under the ground track: 63° 
and 45° port of the ground track; on the 
ground track; and 63° and 45° starboard of the 
ground track. All predictions were made for 
the XV-15 in a level flyover at 90 knots air- 
speed, 85° nacelle tilt, 95 knots groundspeed, 
250 feet altitude, and gross weight of 13,000 
pounds. Wing off-loading of rotor thrust was 
accounted for. Noise values from the port and 
starboard rotors in isolation were summed to 
give the combined rotor noise signature. The 
difference between the predicted baseline 
EPNL and extra blade EPNL was averaged 
over the five observer locations. This average 
was assumed to be a noise reduction attrib- 
uted to a unit increase in the number of blades 
per rotor. 

Accomplishments. The resultant noise re- 
duction due to one additional blade per rotor 
was 6.2 EPNdB. This correction, when ap- 
plied to the measured XV-15 ground contours 
(see figures), corresponds to a shrinking of the 
65 DNL contour surface area by a factor of 
0.23 for both takeoff and approach conditions. 


Significance. This study has demonstrated 
a capability to use the ROTONET system 
noise prediction code to predict trends and 
magnitudes of resulting noise reduction from 
changing rotor system design parameters. 

Status/Plans. Future plans will examine 
the predicted noise results to explain the 
differences in relative magnitudes of the noise 
corrections and to apply the approach to other 
combinations of rotor head design and operat- 
ing parameters. 

Robert A. Golub 
Acoustics Division 
Langley Research Center 
(804) 864-5281 
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Figure 1. XV-15 collective actuator loads trends with configu ration 
mo difi cations - comparison of predicted and measured loads at low 
speed. 
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Figure 2. XV-15 collective actuator loads correlation using a modified 
version of CAMRAD/JA -Helicopter Mode, MSL, 35 to 80 KTS 


Figure 6-7. XV-15/ATB Flight Investigations: Advanced Technology Blades Loads, 

Stability, and Performance 
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6-7 XV-16/ATB Flight Investigations: Advanced Technology 
Blades Loads, Stability, and Performance 


Objective. To investigate advanced technol- 
ogy applied to tiltrotor blades for its impact on 
the performance, dynamics, loads, and 
aeroelastic stability of the advanced Technol- 
ogy Blades. These data provide a configura- 
tion perturbation for the development and 
validation of analytical prediction codes. 

Approach. Advanced aerodynamics and 
mixed matrix composite materials are com- 
bined in the design of the Advanced Technol- 
ogy Blades (ATB) for XV-15 Tiltrotor Re- 
search Aircraft. The highly twisted blades 
are configured with a compound planform, 
thin tips, and an increased solidity to improve 
hover lift capability without degrading per- 
formance at high speeds in the airplane mode. 
Flight tests of the ATB will support the devel- 
opment of a broad range of analytical method- 
ologies. 

Accomplishments. Initial flight testing re- 
vealed unanticipated high oscillatory loads in 
the XV-15 rotor control system. Subsequent 
analyses identified potential control system 
and blade modifications to reduce loads and 
improve structural dynamic stability. To vali- 
date the analyses, variations of tip twist, 
blade chordwise, e.g., control stiffness, and 
blade sweep have been examined in flight. 
The most significant factor in the reduction of 
loads tested to date is the increase of the cyclic 
swashplate stiffness, accomplished by a tem- 
porary modification, which displaced the natu- 
ral frequency of the system sufficiently far 
from the rotor 3-per-rev excitation to enable 
evaluation of the ATB’s performance and loads 
in all operational modes. Also, improvements 
to the Stability and Control Augmentation 
System (SCAS) have been installed for flight 
evaluations. 

Significance. The measured XV-15/ATB 
control loads flight investigation revealed limi- 
tations and deficiencies in current analytical 
and design methods. The flight data provides 


an important configuration perturbation (from 
the original XV-15 metal blades) to develop 
and validate predictive codes. 

Status/Plans. A scheduled 200-hour mqjor 
transmission inspection and overhaul will 
follow the completion of the FY 1991 XV-15/ 
ATB flight tests. A vibration survey will be 
conducted immediately following the flight 
activity to document the airframe resonant 
frequencies in support of the dynamics and 
load analyses. Control system and blade 
modifications to increase stability and reduce 
loads will be performed during the inspection 
period. A full spectrum of tiltrotor flight 
experiments in support of the advanced 
tiltrotor transport technology activity mil 
commence in late FY 1992 after the inspec- 
tion/overhaul effort is completed. 

Brent Wellman, Martin Maisel 
Flight Experiments Branch 
Ames Research Center 
(415) 604-6573/6372 
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Figure A XV-15/ATB collective actuator loads correlation using 
a modified version of CAMRAD/JA • Helicopter Mode, 35 to 80 KTS 



Azimuth (deg) 

Figure B XV-15/ATB Normal bending moment correlation using a modified 
version of CAMRAD/JA showing variation in correlation for different 
wake models • Helicopter Mode, MSL, 35 KTS. 


Figure 6-8. Analysis of Tiltrotor Phenomena Utilizing CANRAD/JA 


98 





6-8 Analysis of Tiltrotor Phenomena Utilizing CAMRAD/JA 


Objective. To evaluate the accuracies of the 
state-of-the-art comprehensive rotorcraft 
analysis codes, such as CAMRAD/JA (Com- 
prehensive Analytical Model of Rotorcraft 
Aerodynamics and Dynamics/Johnson Aero- 
nautics), in the prediction of rotor system 
performance, loads, and stability and to iden- 
tify and address sources of inaccuracies to 
enhance the analysis and design of advanced 
rotor systems. 

Approach. Extensive correlation efforts iden- 
tified the limitations of the current state-of- 
the-art in comprehensive rotorcraft analysis 
codes. Improvements in the comprehensive 
codes will be judged by respective improve- 
ments in correlation. Advanced high-speed 
tiltrotor blade optimization efforts also have 
been initiated using a rigorous two-point op- 
timization of the XV-15/ATB (Advanced Tech- 
nology Blades) rotor system. The validity of 
the results of the optimization effort will be 
strengthened by improvements in correlation 
of the current state-of-the art comprehensive 
codes. 

Accomplishments. Extensive correlation 
efforts have been performed and reported 
using predictions generated from a modified 
version of CAMRAD/JA and XV-15/ATB rotor 
correlated with control system loads flight 
data, as shown in Figure A. This effort iden- 
tified the importance of accounting for control 
system flexibility. The modified version of 
CAMRAD/JA reflected the flexibility, and this 
resulted in improved correlation. This modi- 
fied version of CAMRAD/JA is currently be- 
ing used to generate results for the advanced 
high-speed tiltrotor blade optimization effort. 
The predicted performance and stability gains 
resulting from this effort will reflect the im- 
proved accuracies obtained with the modified 
version of CAMRAD/JA. 


Extensive correlation results have been ob- 
tained utilizing CAMRAD to predict airloads 
for the AH-1G and correlating them with AH- 
1G airload measurements. 

A major effort was also undertaken to up- 
grade the UH-60 CAMRAD model through 
improved input and input source documenta- 
tion as requested by UH-60 PEER Review. 

Significance. Improvements in correlation 
will strengthen the validity of the optimiza- 
tion work and lead to a realistic rotor blade 
design for advanced helicopters and high- 
speed tiltrotor configurations. Improvements 
in performance and aeroelastic stability must 
be achieved to enhance the feasibility and 
various measures of aircraft efficiency for 
rotorcraft at high speeds. 

Status/Plans. Further aeroelastic stability 
and airloads correlation efforts are planned 
for the XV-15/ATBs, the AH-IG Cobra, and 
the UH-60 Phase II program to ensure that 
inaccuracies are identified and comprehen- 
sive rotorcraft analysis codes are improved 
and reflected in the advanced high-speed 
tiltrotor blade optimization effort. 

J. J. Totah, J. F. Madden III 
Rotorcraft Technology Branch 
Ames Research Center 
(415) 604-4126 
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Figure 6-9. Comparison of Measured and Predicted Acoustic Signal for a Model Helicopter 
at a Strong Blade-Vortex Interaction Noise Condition 
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6-9 Blade-Vortex Interaction Noise Prediction Validation 


Objective. To understand rotorcraft noise 
source mechanisms and to develop and vali- 
date noise prediction codes, specifically to 
validate a Langley-developed noise predic- 
tion code for blade-vortex interaction (BVI) 
noise using measured rotor blade surface pres- 
sures. 

Approach. Far-field noise and blade pres- 
sure data obtained from a model rotor test 
were employed to vabdate the rotor noise 
prediction code, WOPWOP, for BVI noise. 
The test was a cooperative effort between 
Sikorsky Aircraft, the United Technologies 
Research Center (UTRC), NASA Langley and 
Ames Research Centers, and the U.S. Army 
Aeroflightdynamics Directorate. The model 
rotor was a 1/6 scale, 4-bladed, swept-tip 
design. The rotor blades were instrumented 
on both the upper and lower surfaces with 176 
recess-mounted pressure transducers. This 
large number of transducers, combined with a 
high data digitization rate, resulted in a high 
resolution aerodynamic database. The blade 
pressure resolution was sufficient to capture 
short time events, such as BVI. Noise predic- 
tions were made using the measured blade 
pressures as input into the noise prediction 
code WOPWOP. The predicted noise time 
histories were then compared with the mea- 
sured time histories. 

Accomplishments. Measured and predicted 
acoustic time histories have been compared 
for a number of rotor operating conditions and 
measurement locations. A comparison be- 
tween predicted and measured results for a 
BVI case at a moderate flight speed of 85 
knots is presented in the figure. The top 
figure illustrates the microphone position 
under the flight path, and the lower figures 
show the measured and predicted acoustic 
results. The results are presented for one- 
rotor period and clearly show the impulsive 
BVI noise from each of the four blades. The 
corhpaiison is good, allowing for blade-to- 


blade variability in the measured data. The 
maj or features of the measured data are fairly 
well reproduced in both shape and magni- 
tude. 

Significance. The rotor noise prediction 
code WOPWOP can predict BVI noise given 
accurate high resolution blade loads data as 
input. This is significant because noise pre- 
dictions made with predicted blade loads have 
not been in good agreement with measured 
noise data. The recently obtained acoustic 
and blade pressure data will add immensely 
to the understanding and development of pre- 
dictions for both rotor blade loads and acous- 
tics. 

Status/Plans. Noise predictions for a variety 
of operating conditions will be made and com- 
pared with measured data. In addition, noise 
predictions from predicted blade loads will be 
made and compared with measured results. 
The predicted and measured blade loads will 
be compared and deficiencies in the predicted 
blade loads will be identified. This will allow 
for the development of an improved blade 
loads prediction and thus a better noise pre- 
diction capability. 

Casey L. Burley 
Acoustics Division 
Langley Research Center 
(804) 864-3659 
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6-10 High Resolution Rotor Blade Loads for Blade-'Vbrtex 
Interaction Noise Prediction 


Objective. At present, the missing element 
in providing Langley with the first fully ana- 
lytical prediction capability for helicopter main 
rotor noise is the specification of detailed 
blade loading. A study was undertaken to 
provide this detail, through a unique recoding 
and post-process modeling of the rotor perfor- 
mance code Comprehensive Analytical Model 
of Rotorcraft Aerodynamics and Dynamics 
(CAMRAD). 

Approach. The accurate prediction of impul- 
sive blade-vortex interaction (BVI) noise, as 
well as the prediction of low frequency loading 
noise when BVI occurs, requires accurate 
blade loading definition at azimuthal steps of 
1° or less, as well as small radial blade incre- 
ments. The well-known and often used per- 
formance code CAMRAD has sophisticated 
free-wake and aeroelastic blade modeling, 
but it has a 15° azimuthal resolution limita- 
tion. The approach taken to obtain a high 
resolution version of CAMRAD was to first 
increase the code resolution version of 
CAMRAD to the intrinsic limit (found to be 
10° azimuth). Post-processing codes were 
developed to determine blade positions and 
tip vortex strengths and trajectories at very 
small incremental steps of 1°. At each incre- 
mental step, the induced velocities and blade 
sectional loadings were calculated. 

Accomplishments. The CAMRAD code 
modifications and the post-processing code 
for the free wake and first order unsteady 
aerodynamics have been completed. An ap- 
plication of the code is shown in the figure, 
where the rotor disk loading, the lift encoun- 
tered by the blades as they rotate, is plotted. 
Compared with the 10° resolution results, the 
1° high-resolution results offer an order of 
magnitude increase in the level of detail and 
understanding of the blade-wake interactions. 
The loading reveals the tip vortex positions 
and the locations where they pass through the 
rotor disk. The high-resolution code has been 


applied successfully in support of an experi- 
mental effort to reduce BVI noise through the 
use of higher harmonic control (HHC) of blade 
pitch. After incorporating HHC capability 
into the code, the effect of HHC on BVI noise 
levels was accurately predicted. 

Significance. The effort has produced a 
more sophisticated and versatile prediction 
capability for the harmonic loading and im- 
pulsive BVI noise. The long experience of 
government and industry researchers in the 
use of CAMRAD for aerodynamic and 
aeroelastic calculations makes a high resolu- 
tion version particularly attractive. 

Status/Plans. Efforts are under way to com- 
plete the loads definition by installing a near 
wake model, as well as an advanced unsteady 
aerodynamics model. Detailed noise directiv- 
ity predictions are to be obtained by using the 
high-resolution loading calculations in 
Langley’s WOPWOP rotor noise code, and 
these will be compared with existing acoustic 
data sets. The code will continue to be applied 
in support of the noise reduction research on 
higher harmonic control. 

Thomas F. Brooks 
Acoustics Division 
Langley Research Center 
(804)864-3634 
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6-11 Ray-Acoustics Approach to Fuselage Scattering of 

Rotor Noise 


Objective. To predict the acoustic field of 
installed rotors and propellers. The specific 
objective of the present research is to develop 
a numerical method for the aircraft fuselage 
scattering of noise from a rotor. 

Approach. A numerical implementation of 
ray acoustics theory in conjunction with a 
paraxial ray approximation is used to model 
the scattering of an impulsive noise source by 
a fuselage. This is a two-point ray tracing 
problem, where both the source position and 
observer location are known in space. This is 
modeled as a boundary value problem that 
requires an iterative solution, and is very 
expensive to solve computationally. In order 
to improve numerical efficiency, paraxial ray 
approximation (PRA) is used. This is devel- 
oped from a Gaussian beam method (GBM), 
where Gaussian beams are high-frequency 
asymptotic time-harmonic solutions to the 
wave equation, and are concentrated close to 
rays. A GBM solution requires no iteration, 
but does involve a superposition of a large 
number of beams to compute the noise at a 
given observer location. The PRA evaluates 
the time-travel ray path and amplitudes both 
along the ray and in its vicinity. So a single 
ray can be used to obtain the solution at the 
observer, providing it passes sufficiently close 
to the observer position. 

Accomplishments. An example of the fuse- 
lage scattering result is given in the figure for 
rotor impulsive noise source modeled as a 
point load rotating in a plane above the scat- 
tering body. The geometry of the model is 
shown in the figure. The point loading rotates 
at angular velocity ft, on a disk of radius a, at 
a height H above the body. A spherical body 
of radius b is depicted in the figure. Numeri- 
cal results are also shown in the figure for 
three different scattering bodies at an ob- 
server angle 0 of 45°, b/a=0.5, H/a=0.5, and a 
source Mach number (ft a/speed of sound) of 
0.7. The result for a free-field case, where 


there is no fuselage, and a single impulse is 
seen. For the case of a sphere, the initial 
impulse is received cleanly, and is followed by 
a smaller, broader pulse due to scattering 
from the sphere. Where the sphere is re- 
placed by an ellipsoid, a larger amplitude 
scattered pulse is seen. These results are 
highly dependent on observer position in ad- 
dition to the shape of the scattering body. 

Significance. Acoustic scattering from wing 
and fuselage surfaces for full aircraft configu- 
rations is always a concern in experimental 
aeroacoustics. In addition, model supports, 
microphone struts, and other wind tunnel 
hardware can cause reflection and scattering 
problems. This technique can be used to help 
understand measured acoustic results, and 
especially to explain differences between mea- 
sured results and theoretical predictions, as 
predictions usually do not include scattering 
surfaces. 

Status/Plans. A “user-friendly” version of 
the numerical implementation of this theory 
is being developed. This will be a useful 
diagnostic tool to the experimental acousti- 
cian, who will be able to define the specific 
geometrical setup of a particular test. 

Michael A. Marcolini 
Acoustics Division 
Langley Research Center 
(804) 864-3629 
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Figure 6-12. Comparison of Mid-Frequency Noise Contours for the BO-105 Rotor 
in Normal Flight and With Higher Harmonic Control 
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6-12 Rotor Impulsive Noise Reduction Using Higher 

Harmonic Control 


Objective. Impulsive blade-vortex interac- 
tion (BVT) noise is one of the most objection- 
able types of helicopter noise. A previous 
rotor noise experiment in the Langley Tran- 
sonic Dynamics Tunnel (TDT) proved that 
higher harmonic control (HHC) of blade pitch 
could reduce overall sound power. This ex- 
periment was undertaken to obtain the effect 
of HHC on the directivity of BVI noise. 

Approach. The test, conducted in the ane- 
choic German-Dutch Wind Tunnel (DNW), 
was a cooperative effort between NASA Lan- 
gley, the German aerospace research estab- 
lishment DLR, and the European helicopter 
companies MBB and Aerospatiale. A large 
database was acquired for a 40% scale model 
of the MBB BO- 105 main rotor. Noise and 
vibration measurements were made for a 
range of operating conditions where higher 
harmonic blade pitch was superimposed on 
the normal collective and cyclic trim pitch. 
HHC control modes of 3 per rev (3P), 4P, and 
5P at various amplitudes and phases were 
tested, as well as special mixed 3P/4P/5P 
modes. The noise was measured over a large 
measurement plane underneath the model 
rotor using a traversing in-flow microphone 
array. 

Accomplishments. It was found that the 
BVI noise directivity could be significantly 
modified by using HHC. The figure presents 
BVI noise contour plots over the measure- 
ment plane for a particular descent flight 
condition, advance ratioof 0.15 and an equiva- 
lent descent angle of 6°, which is considered 
typical for full-scale helicopters. On the left, 
a baseline case is shown where HHC is not 
used. The levels given are mid-frequency 
values and represent the BVI portion of the 
total noise. The plot to the right is the noise 
contour for a particular HHC pitch schedule 
which causes a significant 6 dB reduction in 
the advancing side lobe levels, although the 
retreating side increases by 2 dB. The rela- 


tive levels and directivity can be modified by 
changing the area of the rotor disk in which 
the HHC is employed. 

Significance. This database is the first 
available on BVI noise radiation and directiv- 
ity using HHC. The results support the key 
conclusion of the previous Langley test that 
HHC can be used to reduce BVI noise for 
descending flight conditions where BVI noise 
is most important. These results are doubly 
significant because of the two different rotor 
types used in the investigations. The present 
rotor is hingeless with dynamically scaled 
blades, while the rotor tested in the TDT was 
fully articulated with much stififer blades. 

Status/Plans. Analysis is continuing to es- 
tablish more fully the BVI noise reduction 
benefits versus operational and HHC param- 
eter variations. In addition, the vibration and 
low frequency noise data will be used to help 
establish the practicality of particular con- 
trols. The results will be compared with 
acoustic and vibratory predictions based on 
an improved high-resolution loading version 
of the CAMRAD rotor performance program 
currently under development at Langley. 

Thomas F. Brooks 
Acoustics Division 
Langley Research Center 
(804) 864-3634 
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Figure 6-13. Joint Army/Boeing Research Program 
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6-13 Cooperative Army/NASA/Boeing Pressure-Instrumented 

Rotor Program 


Objective. To (1) investigate the basic for- 
mation of the rotor wake directly behind the 
advancingblade, (2) validate a circulation box 
technique for measuring rotor loads, and (3) 
capture a blade-vortex interaction occurrence 
in the rotor wake. 

Approach. A wind tunnel test was con- 
ducted in the 14- by 22-foot subsonic tunnel 
using the Boeing Single Rotor Helicopter 
(SRH) drive system, a pressure- instrumented 
model of the Boeing 360 rotor and a two- 
component laser velocimeter. Measurements 
of the dynamic blade surface pressures and 
the flow-field velocities were obtained for sev- 
eral flight conditions. 

Accomplishments. Measurements of the 
blade surface pressures using the 66 trans- 
ducers installed in the four rotor blades and 
velocity measurements using the laser were 
made for nine flight conditions ranging from 
high-speed forward flight to descent condi- 
tions. Velocity data were obtained for over 
1600 flow-field measurement locations. In all 
cases, average pressure and velocity data as 
well as time-dependent data were acquired. 

Significance. These data will be used to 
correlate surface pressure measurements with 
wake velocity measurements to determine 
the nature of the forcing function for rotor 
blade airloads. These data also provide a link 
between the acoustic measurements obtained 
in the German Dutch Wind Tunnel (DNW) 
and the flow-field phenomena which gener- 
ates rotor noise. The data identifying the 
formation of the rotor tip vortex and wake will 
be used to improve analytical predictions of 
rotor wakes. The velocity data around the 
blade in a box will be used to validate the 
circulation measurement technique, and, if 
successful, will lead to the capability to mea- 


sure blade loading without expensive and 
troublesome pressure transducers in the ro- 
tor blades. 

Status/Plans. The pressure data are being 
digitized and reduced at Boeing Helicopters. 
The laser velocimeter data will be processed 
at Langley. A joint publication is anticipated 
next year. Correlation with analysis will be 
done both at Boeing with their proprietary B- 
65 code and at Langley using CAMRAD-JA 
with FPR coupling. 

Susan L. Althoff 
Subsonic Aerodynamics Branch 
Langley Research Center 
(804) 864-5059 
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6-14 Rotor/Fuselage Aerodynamic Interactions Program 


Objective. To quantify rotor/fuselage aero- 
dynamic and acoustic interactions through 
wind tunnel testing and computation. The 
program results will assist in the design of 
rotorcraft systems by taking advantage of 
mutually beneficial aerodynamic interactions. 

Approach. To achieve the program objec- 
tive, a full-scale test of a Bell 412 rotor system 
with a simple axisymmetric fuselage was con- 
ducted in the 40- by 80-foot wind tunnel at 
Ames. Small-scale tests with similar fuselage 
shapes have been conducted in the past. 

Accomplishments. Data from the 40- by 80- 
Foot Wind Tunnel test has been reduced and 
analyzed. A paper based on the test results 
was presented at the 1991 American Helicop- 
ter Society Forum. 

Significance. Results from the wind tunnel 
test revealed that the magnitude of the aero- 
dynamic interactions on the simple fuselage 
were small. The importance of rotor shaft and 
hub interactions was also shown. In addition, 
comparisons between full-scale and small- 
scale results demonstrated the importance of 
model configuration and installation differ- 
ences. For example, the figure shows a com- 
parison of the rotor effect on the upper surface 
pressures for full- and small-scale data. In 
the figure, (ACp)R represents the difference 
in pressure coefficient that is due to the rotor. 
The large differences between the full-scale 
and 1/6-scale results are attributed mainly to 
the lack of a rotor mast in the 1/6-scale test. 
The reason for the offset between the full- 
scale and 0. 15-scale data is unclear, although 
potential causes include the unsealed hub 
and significantly different model mounting 
system used in the 0.15-scale test. 


Status/Plans. Plans are to complete the 
documentation of the full-scale wind tunnel 
test results. 

Thomas R. Norman and 
Rotorcraft Aeromechanics Branch 
Ames Research Center 
(415) 604-6653 
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Helicopter Fuselage 
(source panel model) 


Rotor Wake 
(vortex lattice model) 






RWF code output (LaRC-Berry) 


Figure 6-15. Rotor-Wake-Fuselage Code Development 
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6-15 Rotor-Wake-Fuselage Code Development 


Objective. To predict the aerodynamic envi- 
ronment of a realistic helicopter configura- 
tion which accounts for the interaction of the 
rotor wake and the helicopter fuselage. 

Approach. This program has been an in- 
house development effort. Existing code ele- 
ments for computing the wake from a single 
rotor blade and the influence of a nonlifting 
source panel have been integrated into a code 
that computes the interacting flow field of 
multiple rotor blades with cyclic pitch, their 
shed wakes, and a nonlifting fuselage. The 
resulting computer code Rotor-Wake-Fuse- 
lage (RWF) has been under preliminary de- 
velopment since 1986. 

Accomplishments. During FY 1991, pre- 
liminary copies of the implementing code 
(RWF) have been provided to two requesting 
U.S. helicopter companies. A version of the 


code is under modification for production op- 
eration on the Langley Supercomputer Net- 
work System (SNS). The code is being used to 
assess fuselage effects on rotor inflow. 

Significance. The computer code RWF is a 
valuable tool for investigating of realistic fu- 
selage effects on helicopter configurations. 

Status/Plans. After validation of the SNS 
production version of the RWF code, the code 
will be made available to interested U.S. us- 
ers. A low level of continuing development is 
planned to respond to requirements of the 
target user community. 

John D. Berry 

Subsonic Aerodynamics Branch 
Langley Research Center 
(804) 864-5090 
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Figure 6-16. University of Maryland/Ames Research Center Bearingless Rotor 

Aeroelastic Analysis 
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6-16 Five-Bladed Bearingless Rotor Program 


Objective. To document the performance, 
loads, and stability of five-bladed bearingless 
rotors, and to develop and validate compre- 
hensive analyses to accurately predict and 
design new rotor configurations. 

Approach. Testing of the McDonnell Dou- 
glas Helicopter Company (MDHC) Advanced 
Bearingless Rotor (MDAR) and the Sikorsky 
Five-Bladed Bearingless Main Rotor (SBMR) 
was conducted in the 40- by 80-foot wind 
t unne l through the flight envelope up to 200 
knots. A comprehensive database on rotor 
performance, blade, flexbeam and torque tube 
loads, and aeroelastic stability was acquired. 
The results of these test programs will assist 
in the development and validation of the Uni- 
versity of Maryland Advahced Rotorcraft 
Code (UMARC). 

Accomplishments. A joint cooperative pro- 
gram was established for the MDAR tests. 
The first test was scheduled for the fourth 
quarter of CY 1991 in the 40 x 80 on the 
MDHC test stand. The second phase of the 
program with a pressure instrumented rotor 
has been initiated. A joint cooperative pro- 
gram was established with Sikorsky for test- 
ing of the SBMR in the 40 x 80 on the Rotor 
Test Apparatus in the second quarter of CY 
1992. The UMARC Code has been docu- 
mented and correlation is in progress for the 
SBMR stability data. The MDHC HARP rotor 
has been modeled and the MDAR modeling is 
awaiting data release from MDHC. 

Significance. The correlation efforts with 
UMARC and experimental stability data have 
been encouraging. Correlation with the SBMR 
hover stability data has been good. Correla- 
tion with a four-bladed model scale bearingless 


rotor wind tunnel test data for shaft fixed and 
shaft free (see figure) have shown good agree- 
ment. 

Status/Plans. Wind tunnel test planning is 
continuing for both the MDAR and SBMR. 
Development and validation of UMARC con- 
tinues as data becomes available. 

Steve Jacklin 

Rotorcraft Aeromechanics Branch 
Ames Research Center 
(415) 604-4567 
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Figure ® Calculated Mach contours for AH-1G blade at a high speed 
flight condition where transonic flow is encountered. 


Figure 6-17. Rotor Airloads Correlation Using CFD/Lifting-Line Methods 
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6-17 Rotor Airloads Correlation Using CFD/Lifting-Line Methods 


Objective. To assess the predicting capabil- 
ity of current state-of-the-art rotorcraft aero- 
dynamic codes on modem rotors by compar- 
ing results with existing flight test data. 

Approach. The comprehensive rotorcraft 
code, CAMRAD/JA, was used to model the 
aerodynamic and dynamics of the complete 
rotorcraft configuration. To obtain the de- 
tailed blade airloads, the Full Potential Rotor 
Code (FPR) was coupled with CAMRAD/JA in 
an iterative loop. The final solution was the 
predicted airloads on the rotor system which 
were then compared with flight measure- 
ments at different operating conditions . 

Accomplishments. A detailed airloads cor- 
relation effort on an AH-1G Cobra helicopter 
has been completed. Computational methods 
were used to calculate the tjlade airloads and 
give insight into the nature of blade vortex 
interaction loads at low advance ratios. Fig- 
ure A shows the correlation between calcu- 
lated measured blade section lift during blade 
vortex interaction encounters. As shown in 
the figure, the section lift was calculated for 
several vortex core sizes. The effect was 
minimal for this case, implying the vortex did 
not pass near the blade. At high advance 
ratios, the same approach was used to predict 
the transonic compressibility effects that are 
present on the advancing side of the rotor 
disk. The use of computational fluid dynam- 
ics (CFD) allowed both a quantitative and 
qualitative analysis of the flow phenomena in 
a very efficient manner. Figure B shows 
predicted Mach contours at a high advance 
ratio where transonic flow is present. Also, 
computational fluid dynamics methods have 
been used to predict the rotor airloads on the 
UH-60 Blackhawk and to compare them with 
existing experimental data from the German- 
Dutch Wind Tunnel (DNW). This CFD model 
is being developed in preparation for Phase II 
of the UH-60 Airloads Program, currently 
under way at Ames. 


Significance. The accurate prediction of 
rotor airloads is one of the biggest challenges 
in the field of theoretical aerodynamics and 
one which remains to be solved. Current 
rotorcraft aerodynamic codes make use of 
many approximations and assumptions that 
ultimately must be replaced by more accurate 
methods. To determine the validity of such 
methods and improve upon them, correlation 
work with experimental data is needed. Such 
improved codes will help in the design of 
better rotor systems. 

Status/Plans. CFD methods are currently 
being used to model the UH-60 Blackhawk 
rotor system from hover to high advance ra- 
tios . Preliminary results have been compared 
with model wind tunnel data at different 
flight conditions. These methods will help in 
the prediction of rotor airloads during Phase 
II of the UH-60 Airloads Program. 

Francisco J. Hernandez 
Rotorcraft Technology Branch 
Ames Research Center 
(415) 604-1322 


117 




Ill;;-- ; : ;|1 

mmmamm 

5 : - .: V ;: :: l--,-;- -v 

.. ■ . . ' ■ 


: ; ; ; 





6-18 Rotor Flow Field Investigation 


Objective. To provide a comprehensive da- 
tabase documenting the velocity flow field 
generated by a thrusting helicopter rotor in 
forward flight. This database is intended to 
provide a means to validate the accuracy of 
the velocity and rotor load predictions of rotor 
computational methods. 

Approach. A wind tunnel test was con- 
ducted in the 14- by 22-foot subsonic tunnel 
using the 2-Meter Rotor Test System and a 
two-component laser velocimeter to measure 
the velocity induced by the rotor above the 
plane of the rotor blade tips. Measurements 
were also made in two planes through the 
separated flow region behind the hub. A 
second laser was used to generate a light 
sheet to provide flow visualization of the vor- 
tex structure in the rotor wake. 

Accomplishments. Measurements of the 
rotor-induced velocity were made at 180 loca- 
tions in planes above the rotor blade tips at a 
height of 0.50 chord for an advance ratio (ratio 
of freestream velocity to rotor tip speed) of 
0.23. Measurements with and without blades 
were made in the region of separated flow 
behind the rotor hub at 15 longitudinal sta- 
tions, with 15 points at each station. In all 
cases, average velocity data as well as time- 
dependent data were acquired. Flow visual- 
ization studies were conducted for the rotor 
wake region above the rotor disk and down- 
stream of the rotor for advance ratios of 0.15, 
0.23, and 0.30. 

Significance. The rotor inflow measure- 
ments, which were obtained during this test 
program, demonstrate the first production 
use of the newly developed Frequency Do- 
main Processors (FDPs) for the laser veloci- 
meter data acquisition system. The hub sepa- 
ration measurements provide the basis to 
define the extent of the separated flow behind 
the rotor hub for computational rotor-fuse- 


lage interaction code validation. The flow 
visualization work provides the beginning of 
a visual database of the rotor vortex wake 
structure. 

Status/Plans: The data obtained in this 
program will be used in-house for correlation 
with rotor performance codes and rotor-fuse- 
lage interaction codes. The flow visualization 
video data will be digitized and enhanced for 
further analysis by image processing soft- 
ware. 

Susan L. Althoff, Joe W. Elliott 
Subsonic Aerodynamics Branch 
Langley Research Center 
(804) 864-5059 
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Chapter 7 

Fighter/Attack Aircraft 


The objective of the Fighter/Attack Aircraft Program is to provide the required technologies for 
fighter/attack aircraft to achieve efficient, sustained supersonic cruise and maneuver perfor- 
mance; efficient store carriage and deployment at supersonic speeds; increased agility at 
subsonic speeds and acceptable handling qualities at extreme angles-of-attack; and short takeoff 
and vertical landing (STOVL) operation. 


Improved prediction methods and/or experimental techniques are now available for high-lift 
aerodynamics, propulsion integration, weapon carriage, supersonic store cavity and separation 
aeroacoustics, integrated flight controls, and systems design for fighter/attack aircraft. Wind 


thrust vectoring for propulsive control at extreme angles-of-attack. In addition, powered lift 
systems consistent with the operation of advanced STOVL aircraft have been identified, and 
configuration studies have been conducted to assess the impact of integrating these systems with 
supersonic airframe designs. 


The Fighter/Attack Aircraft Program is focused on (1) CFD modeling and validation of cavity 
effects and near-field trajectory simulation for weapons launch during maneuver and modeling 
of 3-D flow fields for vehicles at high angles-of-attack with separated flow; (2) development of 
innovative propulsive and aerodynamic control systems concepts to provide increased vehicle 
control at angles-of-attack near and beyond maximum lift; (3) improved understanding of control 
requirements for control throughout the relevant flight envelope; (4) large-scale, ground-based 
testing of a STOVL fighter concept using ejector lift/vectored thrust to define critical transition 
aerodynamics; (5) studies of STOVL ground effects including inlet reingestion and ground 
erosion; and (6) integration studies to quantify the impact of emerging technologies and define 
the critical research areas for supersonic STOVL. 

Program Manager: Benjamin Neumann 

OAST/RF 

Washington, DC 20546 
(202) 453-2795 
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F/A-1 8 Aircraft 


Angle of Attack, degrees 


Comparison of small-scale and full-scale tail buffet frequency versus angle of attack. Experimental 
results are for fluctuating pressure measurements at a single point on the vertical tail while the 
CFD point is for the unsteady root bending moment on the tail. 
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Figure 7-1. Tail Buffet Frequency vs Angle-of-Attack 







7-1 Full-Scale Tfest of F/A-18 at High Angles-of- Attack 


Objective. To examine the aerodynamics of 
high angle-of-attack flight of fighter aircraft. 
Particular emphasis is placed on understand- 
ing the LEX-vortex (Leading Edge Extension) 
induced tail buffet and on developing a pneu- 
matic forebody flow control concept for lateral 
control at high angle-of-attack. The test has 
also provided important data for computa- 
tional fluid dynamics (CFD) validation and 
resolution of scaling issues, which have be- 
come apparent in comparing small-scale wind- 
tunnel data with flight results. 

Approach. A retired Navy F/A-18A has been 
converted into a wind-tunnel model for test- 
ing in the 80- by 120-foot wind tunnel. The 
model is heavily instrumented to measure 
time-averaged forebody pressure distributions 
and the unsteady pressure distribution on 
one of the vertical tails. In addition, the nose 
of the aircraft has been modified with the 
addition of jet-type and slot-type blowing for 
flow control at very high angle-of-attack. 

Accomplishments. Preliminary data from 
the wind-tunnel tests matched very well with 
flight data for overall forces and moments as 
well as for the radome and forebody pressure 
distributions. The unsteady vertical tail pres- 
sure data will be used to address wind tunnel 
scaling issues, to identify flow structures and 
dominant frequencies in the burst vortex, and 
to validate CFD analyses. Preliminary analy- 
sis of these data indicates that the frequency 
content of the fin pressure measured in the 
full-scale wind-tunnel test is close to that 
observed in several small-scale experiments. 
Similar pressure measurements are not yet 
available from the flight experiment. 

The pneumatic forebody flow control has been 
shown to be effective for practical blowing slot 
lengths and blowing rates. At a given angle- 
of-attack, the yawing moment is linear over a 
wide range of blowing mass flow coefficient. 
At a 50° angle-of-attack, the pneumatic sys- 


tem can provide more yawing moment than a 
full rudder deflection does at low angles-of- 
attack. 

Significance. The pressure distributions 
measured on the full-scale model will provide 
a rich source of CFD validation data, which 
will help to improve the accuracy of 3-DNavier- 
Stokes computations. Current CFD efforts 
are ongoing to compute the unsteady burst- 
vortex induced pressures on the F-18 tail. 
The wind tunnel measurements will also be 
used in the development and validation of 
CFD for unsteady flows. 

It has been demonstrated that pneumatic 
forebody flow control can generate control 
yawing moments of sufficient magnitude to 
provide directional control at extreme angles- 
of- attack where rudder control is no longer 
effective. These full-scale data are critical in 
order to define the planned flight experiment 
of pneumatic forebody flow control system on 
the High Alpha Research Vehicle, HARV. 

Status/Plans. A second wind-tunnel test 
will be conducted in 1992 to optimize the 
forebody flow control system in preparation 
for flight testing, extend the surface pressure 
database to include the wings for use in CFD 
validation and correlation with flight, and to 
refine the measurement of unsteady pres- 
sures on the vertical tails of the aircraft. 

James C. Ross 

Fixed-Wing Aerodynamics Branch 
Ames Research Center 
(415) 604-6722 
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Figure 7-2. Tail Buffet Research 
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7-2 Tail Buffet Research 


Objective. To study the effects of vortical 
flow fields on tail buffeting of modern high- 
performance aircraft, and explore methods of 
reducing buffet levels while maintaining the 
favorable effects of vortex flow on the high 
angle-of-attack aerodynamic characteristics. 

Approach. A cooperative program was es- 
tablished between NASA, the Navy, and 
McDonnell Aircraft Company to conduct low- 
speed wind tunnel tests on a 0. 16-scale F/A-18 
model with dynamically scaled and instru- 
mented vertical tails. This research provided 
data on the unsteady aerodynamic loadings 
and surface pressures of the tails. Changes to 
the configuration aimed at modifying the vor- 
tical flow field were tested to investigate the 
effect on both buffeting response levels on the 
tails and the overall aerodynamic character- 
istics of the configuration. Additional tests 
were conducted to obtain detailed aerody- 
namic information on the impact of the modi- 
fications. 

Accomplishments. Testing in the Langley 
30- by 60-foot tunnel provided a large data- 
base of dynamic response characteristics of 
the vertical tails and steady-state total force 
and moment data for several configurations. 
Correlation of measured tail buffet levels with 
previous, limited industry testing of the 
baseline F/A-18 configuration was very good. 
The study showed that repositioning and/or 
weakening of the primary vortex system can 
reduce tail buffet levels but can also have a 
large impact on the overall lift and stability 
characteristics. In addition, substantial 
amounts of tail buffeting were seen due to 
fuselage and wing wake flow, in the absence of 
a vortex system. Follow-on aerodynamic-only 
testing indicated that many of the stability 
effects were caused by changes in interactions 
between the vortex system and tail surfaces. 

Significance. The results of this test indi- 
cate' that relatively minor changes in an 


aircraft’s configuration can have a significant 
impact on its tail buffet characteristics. Thus, 
consideration of tail buffeting should be ad- 
dressed early in the design stage of an air- 
craft, when configuration changes are still 
possible. Two issues to consider are: (1) re- 
gions of stalled or separated wake flow can be 
a m^jor contributor to tail buffet; and (2 ) 
aerodynamic/stability impact can be due more 
to interactions of the vortex flow field with 
other aircraft surfaces than to changes in the 
vortex flow field itself. 

Also, successful correlation with previous data 
using scaling parameters suggests that the 
same approach may be used to develop wind 
tunnel-to-flight scaling methods. 

Status/Plans. Follow-on investigations of 
various innovative concepts to alleviate tail 
buffeting are planned. Correlation of results 
with research programs of other organiza- 
tions is planned through a workshop on tail 
buffet research in August, 1991, for Aero- 
nautics Technology Technical Panel HTP-5 
members of The Technical Coordination Pro- 
gram. Results from current and future wind 
tunnel tests will be correlated with airplane 
flight test data for development and valida- 
tion of scaling methods. 

Gautam H. Shah 
Flight Dynamics Branch 
Langley Research Center 
(804) 864-1163 
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7-3 Development of High-Angle-of-Attack Nose-Down Pitch Con- 
trol Requirements for Relaxed Static Stability Combat Aircraft 


Objective. To develop high-angle-of-attack 
nose-down control requirements that can be 
used as design guidelines early in the devel- 
opment cycle of advanced combat aircraft 
incorporating relaxed static stability in pitch. 

Approach. A cooperative program was es- 
tablished between NASA and the U.S. Navy 
to perform the research. A systematic piloted 
simulation study was conducted on the Lan- 
gley Differential Maneuvering Simulator to 
develop and assess candidate guidelines. The 
simulation process involved specific assess- 
ment maneuvers and pilot rating approaches 
developed for this application. A key element 
of the investigation was a parametric study 
involving systematic variations of critical 
parameters affecting nose-down control capa- 
bility, particularly the minimum level avail- 
able at high angles-of-attack (Cm*, see ac- 
companying figure). Extensive involvement 
of Navy researchers and pilots permitted ac- 
celeration of the research process. 

Accomplishments. A comprehensive data- 
base of piloted simulation results was gener- 
ated to enable the assessment and develop- 
ment of potential guidelines. The results 
were based on analysis of flight motions and 
pilot comments and ratings obtained during 
the parametric study. These results show a 
strong correlation of pilot rating to the air- 
plane short-term pitch acceleration and pitch 
rate buildup in response to nose-down com- 
mands applied at high-angle-of-attack condi- 
tions. Using these data, candidate design 
guidelines and flight demonstration require- 
ments were defined. Additional simulation 
tests were conducted to verify that these de- 
sign guidelines could be used to predict pilot 
opinion of the aircraft response. 

Significance. The use of relaxed static sta- 
bility in pitch for performance enhancement 
is a key design approach for advanced combat 
aircraft. Successful application of this con- 


cept requires design guidelines for defining 
high-angle-of-attack nose-down control re- 
quirements early in the development process. 
The output from this study is being used by 
the Navy to define these requirements for 
design and flight testing of the F-18 E/F. The 
Air Force is also interested in this study, both 
in terms of potential application to Advanced 
Tactical Fighter and for development of hand- 
book criteria. 

Status/Plans. A preliminary set of guide- 
lines and flight demonstration requirements 
has been developed based on the results ob- 
tained to date. Simulation tests will be con- 
ducted to verify that the design guidelines 
are valid for complex maneuvering and to 
prepare for full-scale flight tests. A limited 
flight evaluation/validation program on atest 
F-18 airplane will be conducted by the Navy 
in September 1991. Detailed flight testing of 
proposed guidelines using the NASA High 
Alpha Research Vehicle (HARV) is planned. 

Marilyn E. Ogbum, John V. Foster 
Flight Dynamics Branch 
Langley Research Center 
(804)864-1136 
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7-4 Investigation of the Tumbling Phenomenon in Aerodynamic 

Configurations 


Objective. To determine the effect of plan- 
form geometry and mass distribution on the 
tumbling characteristics of wings. 

Approach. Static and dynamic force tests 
were conducted in the 30- by 60-foot tunnel 
and free-tumbling tests yere performed in 
the 20-foot Vertical Spin Tunnel using a se- 
ries of generic wing models. Each set consists 
of twelve models with four different leading- 
edge sweep angles and three different trail- 
ing-edge sweep angles. Aerodynamic data 
through 360° angle-of-attack was used to de- 
velop a three-degree-of-freedom simulation 
for calculating tumbling behavior. 

Free-tumbling tests were conducted over a 
range of center-of-gravity and inertial condi- 
tions. A free-to-pitch apparatus is being de- 
veloped and instrumented to enable continu- 
ous rotation testing of a dynamically scaled 
model for single degree-of-freedom analysis. 

Accomplishments. Free-tumbling results 
have shown that higher aspect ratio wings 
generally have a tumble/no tumble boundary 
at a more positive static margin than lower 
aspect ratio wings. This boundary can be 
shifted forward for the low-aspect ratio wings 
by changing the mass distribution. Initial 
results of the simulation show tumbling be- 
haviors similar to those observed in free- 
tumbling tests. 

Significance. As fighter design trends drive 
toward more tumble-susceptible configura- 
tions, the information generated from this 
research will provide design guidelines to 
assist in configuration development. 


Status/Plans. The aerodynamic database 
will be completed. The instrumentation sys- 
tem for the free-to-pitch apparatus is nearly 
complete. Analysis of the free-tumbling re- 
sults is continuing and the tumble simulation 
is being refined. 

Raymond D. Whipple 
Flight Dynamics Branch 
Langley Research Center 
(804)864-1194 
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7-5 F-106 \fortex Flap Flight Experiment 


Objective. To demonstrate the feasibility of 
the vortex flap concept in a free-flight envi- 
ronment to characterize the wing flow field 
and to provide a database for sharp leading- 
edge wings to be used for calibration of design 
and analysis tools. 

Approach. An F-106B aircraft was modified 
by the addition of wedge-shaped, leading- 
edge flaps to the existing wing. The wing flow 
field was characterized using surface pres- 
sure measurements and flow visualization 
data obtained through the use of flow cones, 
surface oil flows, and a light sheet/flow seed- 
ing system. A limited evaluation of the han- 
dling qualities and performance of the modi- 
fied aircraft relative to the baseline aircraft 
was conducted. 

Accomplishments. Flight activities for the 
modified aircraft with the vortex flap at 30° 
and 40° deflections have been completed. 
Aircraft performance data, detailed wing sur- 
face pressure measurements, and surface flow 
visualization were obtained at subsonic and 
transonic test conditions. Light-sheet and 
flow-seeding systems were installed to visual- 
ize the off-body flow and provided significant 
additional details of the wing flow-field char- 
acteristics at subsonic test conditions. The 
accompanying photograph shows the F-106 
aircraft in flight with the vortex flap deflected 
40° and a coating of oil defining surface flow 
patterns on the flap and wing upper surfaces. 

Significance. Flight testing has demon- 
strated the performance enhancement poten- 
tial of the vortex flap and provided consider- 
able insight concerning its flow field and load- 
ing distribution characteristics. 


Status/Plans. Testing of the aircraft was 
completed in April 1991. Analysis and publi- 
cation of the test results are continuing. 

James B. Hallissy, 

W. Elliott Schoonover, Jr., John E. Lamar 
Applied Aerodynamics Division 
Langley Research Center 
(804)864-2865 
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Figure 7-6. Generic Fighter Models Installed in the Unitary Plan Wind Tunnel 
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7-6 Wing Camber Effects on Flap Effectiveness at 

Supersonic Speeds 


Objective. To evaluate camber effects on 
leading- and trailing-edge flap effectiveness 
at supersonic speeds. 

Approach. A cooperative program was es- 
tablished between NASA and General Dy- 
namics to perform the research. Cambered 
and uncambered wing models were built with 
the same planform and flap geometry. The 
wing was designed to meet a broad spectrum 
of performance goals including efficient 
trimmed supersonic cruise performance and 
maneuver capability. The effects of camber 
on the flap performance were then studied 
through a comparison of force and moment 
measurements on both wings. The models 
were tested across the range from 1.6 to 2.16 
in test section I of the Unitary Plan Wind 
Tunnel. The wings were tested with leading- 
edge flap deflections of -4°, -2°, 0°, 5°, 10°, and 
15° and trailing-edge flap deflections of -30°, 
-20° -10°, 0°, and 10°. 

Accomplishments. The experimental re- 
sults show that leading-edge flap deflections 
of up to 5° reduced drag at high lift on the flat 
wing. Small negative leading-edge flap de- 
flections reduced drag on the cambered wing 
at low lift. Trailing-edge flap effectiveness 
decreased with increasing negative flap de- 
flection. Wing camber had little effect on 
trailing-edge flap effectiveness. 

Significance. This study shows that lead- 
ing-edge flap deflections, which are custom- 
arily used at subsonic speeds, can be used to 
obtain drag reduction benefits at supersonic 
speeds. The drag produced by wing camber at 
supersonic speeds can also be reduced by 
negative deflections of leading-edge flaps. The 
study also showed that wing camber has a 
negligible effect on the trim-producing trail- 
ing-edge flap effectiveness at supersonic 
speeds. 


Status/Plans. The study has been com- 
pleted, and results are being documented in a 
formal NASA publication. 

Gloria Hernandez 
Applied Aerodynamics Division 
Langley Research Center 
(804)864-5572 
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Figure 7-7. Outdoor Static Tests of the Full-Scale Ejector-Lift/V ectored-Thrust 

STOVL E-7A Configuration 
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7-7 Outdoor Static Tests of the Full-Scale Ejector-Lift/ 
Vectored-Thrust STOVL E-7A Configuration 


Objective. To test the full-scale STOVL E- 
7A model at Ames Outdoor Aerodynamic Re- 
search Facility (OARF) to measure the thrust 
of the lifting ejectors and vectoring ventral 
nozzle of the vehicle. The tests will also reveal 
the interactions of the plumes from the lifting 
jets beneath the vehicle as they impinge upon 
the ground. The E-7A research program is 
part of a joint technology development pro- 
gram involving Boeing Canada’s De Havilland 
Division and NASA Ames Research Center. 
The E-7A configuration emerged from a Gen- 
eral Dynamics design study of ejector-lift con- 
cepts. 

Approach. Static thrust calibrations of the 
propulsion system in the STOVL E-7A model 
were performed in the test sections of the 
Ames 40- by 80- and 80- by 120-foot wind 
tunnels as part of the forward-flight tests of 
the model. The static test of the model at the 
OARF will yield more accurate measurements 
of static thrust because jet-induced 
recirculation flow patterns in the wind tunnel 
tests due to the presence of tunnel wall bound- 
aries will not be a factor. The two wind tunnel 
tests have provided low- and high-speed aero- 
dynamic data. The outdoor test will complete 
the database on the full-scale model by pro- 
viding static thrust calibrations. The E-7A 
model has been mounted approximately 25 
feet above ground at the OARF. 

Accomplishments. Static tests of the E-7A 
model are complete. The measured thrust 
augmentation ratio of the lifting ejectors is 
approximately 1.6. The augmentation ratio is 
defined as the net ejector thrust divided by 
the thrust of the primary nozzles. This value 
is close to that measured in the test section of 
the 80- by 120-foot wind tunnel. Changes in 
the shape of the base of the fuselage between 
the ejector systems were shown to have an 


effect on the performace of the thrust-aug- 
menting ejectors. Laser light sheet, smoke, 
and infrared flow visualization techniques 
have revealed that the hot ventral nozzle jet 
penetrates forward in a thin layer beneath 
the lower energy ejector jets. 

Significance. Designers of ejector-lift/vec- 
tored-thrust STOVL configurations have been 
working under the assumption that an ejector 
thrust augmentation of 1.6 was achievable at 
full-scale. The level of ejector performance 
and forward-flight transition performance 
during the Ames E-7A tests demonstrates the 
viability of supersonic STOVL concepts which 
employ lifting ejectors. 

Status/Plans. The lateral/directional perfor- 
mance of the E-7A model in hover mode is 
scheduled to be tested at low speeds and high 
yaw angles in the 80- by 120-foot wind tunnel 
beginning in June of 1992. 

Tim Naumowicz, Brian E. Smith 
National Full-Scale Aerodynamics Complex 
Ames Research Center 
(415) 604-6674 
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7-8 Validation of Out-of-Ground Effect Prediction Capability 

for Powered Lift Aircraft 


Objective. To evaluate the capability of 
computational fluid dynamics (CFD) to pre- 
dict the critical flow features for a Short Take- 
Off Vertical Landing (STOVL) aircraft con- 
figuration through comparison of computed 
results with the measurements of a compan- 
ion experiment designed for CFD validation. 

Approach. Previous simulations of critical 
powered lift flow components such as a jet in 
crossflow, which isolates the lift jet/aerody- 
namic interaction during transitional flight, 
showed favorable agreement with experiment. 
To extend the application of the F3D code, a 
simplified but representative STOVL aircraft 
configuration was selected for testing and 
simulation. The model geometry, a Delta 
wing (E-7A) planform with two jets located in 
a blended fuselage, minimizes the geometric 
modeling complexity while retaining the im- 
portant physics of the lift jet/aerodynamic 
interaction. The unsteady Reynolds-aver- 
aged Navier-Stokes equations applied on 
patched or overlapped grids will be used for 
the computations. 

Accomplishments. The transitional flight 
characteristics of the STOVL aircraft configu- 
ration were measured in the NASA Ames 7- 
by 10-foot wind tunnel. The results include 
forces and moments, steady and unsteady 
surface pressures, and jet pressures and tem- 
peratures. Measurements of the flow were 
also made in the tunnel test section upstream 
and downstream of the model and at the jet 
exits to provide boundary conditions for the 
computations. Flow visualization and total 
pressure measurements in the jet plumes 
provide a description of the three-dimensional 
jet efflux flow field. 

Computational emphasis has been placed on 
establishing the numerical modeling require- 
ments for the propulsive flow field. Previous 
thin- Jayer Navier-Stokes computations using 
F3D for a jet in crossflow indicated that the 


computed vortices are less diffuse than in the 
experiment, and the jet wake region is not 
adequately simulated. By retaining viscous 
terms in all three directions, quantitative 
agreement with experiment has been im- 
proved. This modification is incorporated 
into the STOVL model simulation. 

Significance. A comprehensive and well- 
documented set of measurements is provided 
for CFD validation comparisons. 

Status/Plans. Computational grids are be- 
ing developed for the STOVL model geom- 
etry. Computations (out-of-ground effects) 
will be completed and compared with experi- 
ment. Flow field measurements in and out of 
ground effects will be made using laser 
velocimetry. 

Karlin R. Roth 

Fixed Wing Aerodynamics Branch 
Ames Research Center 
(415) 604-6678 
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7-9 Comparison of Shaved and Beveled Fins for Rolling Missile 

Applications 


Objective. To investigate the effectiveness 
of shallow shaving along the leading and 
trailing edges of wrap-around missile fins for 
producing roll on the configuration. 

Approach. An experimental investigation 
was conducted to investigate the aerodynam- 
ics of a tube-launched air-to-surface tactical 
missile configuration that employs planar 
canards and wrap-around tail fins arranged 
in opposing pairs. Fin shaping in the form of 
shallow shaving along the leading and trail- 
ing edges of the tail fins was investigated as a 
candidate for producing the rolling motion on 
this configuration that is needed for controlla- 
bility. A second set of fins having a more 
conventional 45 bevel along the tail fin lead- 
ing edges was investigated for comparison. 
Force and moment tests on this concept were 
conducted in the low Mach number test sec- 
tion of the Langley U nitary Plan Wind Tunnel 
at Mach 1.60, 1.90, and 2.16 at angles-of- 
attack up to 20 and roll angles up to 90°. 

Accomplishments. To be most effective, the 
rolling moments on this configuration should 
be constant over the anticipated flight condi- 
tions. The results of this experiment show 
that the shaved fins were much more effective 
than the beveled fins in producing these fa- 
vorable rolling moment characteristics. As 
shown in the figure, the shaved fins produced 
significantly larger rolling moments , and these 
moments were essentially invariantwith both 
Mach number and angle-of-attack. The bev- 
eled fins, however, produced large variations 
with both parameters, in some cases even 
producing negative rolling moments. 

Significance. Shallow shaving of the lead- 
ing and trailing edges of fins has been shown 
to be an effective roll-producing technique for 
missile configurations that require rolling 
motion for controllability . 


Status/Plans. Supersonic tests on this con- 
cept have been completed. Analysis of the 
data is continuing, and the results will be 
published in a NASA report. Subsonic/tran- 
sonic tests of this concept may be performed. 

Jerry M. Allen 

Applied Aerodynamics Division 
Langley Research Center 
(804)864-5592 
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Augmenter Inlet Flow Field Velocity Vectors 


Figure 7-10. STOL/STOVL Concepts for High-Performance Aircraft 
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7-10 STOL/STOVL Concepts for High-Performance Aircraft 


Objective. To provide advanced aerodynamic 
technology and advanced concepts applicable 
to STOL/STOVL operations of current and 
future high-performance aircraft. 

Approach. Through experimental and ana- 
lytical studies, using powered models in both 
the Subsonic Basic Research Tunnel (SBRT) 
and the 14- by 22-Foot Subsonic Tunnel (14- 
by 22-FST), key problems were examined and 
promising solution concepts for providing 
STOL/STOVL capabilities for advanced air- 
craft were identified. The SBRT and 14- by 
22-FST provide a unique combination of re- 
sources for the study of takeoff and landing 
aerodynamics in and out of ground effect. 

Accomplishments. In FY 1991, two major 
tests were conducted in the 14- by 22-FST and 
the SBRT, which were made operational with 
a moving belt ground plane and data system. 
A detailed investigation of augmenter perfor- 
mance, ground effects, stability and control, 
and a Laser Velocimetry_ (LV) map of the 
augmenter inlet flow field- was completed on 
an advanced VTOL concept shown in the 
accompanying figure. A study of the effects of 
high nozzle pressure ratio (NPR) exhaust 
flows; nozzle type, location, and number; vari- 
ous wing planforms; and moving belt ground 
plane on the aerodynamics of a generic VTOL 
configuration shown in the figure was also 
completed. The SBRT was used to assess the 
effect of a moving belt ground plane on the 
development of the ground vortex flow field 
from a single vertical jet in a cross flow. 

Significance. These three experimental in- 
vestigations have yielded the following sig- 
nificant results: (1) A method of venting the 
lower and upper surfaces to equalize pres- 
sure, which minimized the pitch-up tendency 
typical of configurations during transition, 
was demonstrated using the advanced VTOL 
configuration. (2) Comparisons of the compu- 
tational and LV measured augmenter inlet 


flow fields show very good agreement. (3) The 
data from the generic VTOL test and the 
isolated jet tests in SBRT are under analysis. 
Preliminary analysis of the effect of the mov- 
ing belt ground plane on the ground vortex 
flow field, both in 14- by 22-FST and SBRT, 
show good agreement with results obtained 
using a moving jet in the Vortex Research 
Facility (VRF). These results indicate signifi- 
cant changes in the extent of the ground 
vortex when the exhaust flow penetrated 
against the freestream (belt on or moving jet) 
or against the wall boundary layer (belt off). 

Status/Plans. The results of the ground 
vortex investigations will be analyzed to de- 
termine the effects of wing location, nozzle 
arrangements, and pressure ratio on configu- 
ration aerodynamics in and out of ground 
effect, A detailed quantitative determination 
of the actual extent of the ground vortex will 
be made and compared with the VRF moving 
model results. The SBRT will be completed 
with a model support system to allow model 
testing and an LV system for flow measure- 
ments. 

John W. Paulson, Jr. (804) 864-5071 
Guy T. Kemmerly (804) 864-5070 
Kevin J. Kjerstad (804) 864-5022 
Subsonic Aerodynamics Branch 
Applied Aerodynamics Division 
Langley Resaearch Center 
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7-11 Effect of a Close-Coupled Canard on Wing-Body 

Aerodynamics 


Objective. To accurately simulate the aero- 
dynamics of a canard-wing-body configura- 
tion and investigate the characteristics of the 
canard-wing vortex interaction and break- 
down. 

Approach. Using an extension of the Tran- 
sonic Navier-Stokes (TNS) code, Reynolds- 
averaged thin-layer Navier-Stokes equations 
were solved for the flow about a canard-wing- 
body configuration. Flow was considered fully 
turbulent and turbulence effects were mod- 
eled using the Baldwin-Lomax algebraic eddy 
viscosity model. The surface and flow-field 
grids were generated using the S3D code and 
the 3DGRAPE elliptic solver, respectively. A 
grid refinement study was performed and the 
total number of grid points ranges from 
250,000 to over 1.7 million. 

Accomplishments. Computations were 
made at a transonic Mach number of 0.90, 
angles-of-attack from 0° to 12°, and a Reynolds 
number based on mean aerodynamic wing 
chord of 1 .52 million. Very good agreement of 
computed wing surface pressures and inte- 
grated forces (e.g., lift, pitching moments, and 
drag) with experimental measurements has 
been obtained and indicates that the com- 
puted results are accurate. The accompany- 
ing figure illustrates the effect of the canard 
at a 12° angle-of-attack. Specifically, the 
effect of the canard on the characteristics of 
wing vortex breakdown was noted. A side-by- 
side comparison of the canard-on and canard- 
off cases showed that the wing exhibits a 
stable vortex in the presence of the canard. 
This result indicates the potential of the ca- 
nard to delay wing vortex breakdown and has 
been previously observed in numerous ex- 
perimental studies. 

The surface pressure map also shows the 
significant influence of the canard on the body 
aerodynamics. 


Significance. Utilizing canards for enhanced 
aircraft performance requires a thorough 
understanding of the associated complex aero- 
dynamics. Accurately predicting the nonlin- 
ear effects of the canard, including the ca- 
nard-wing vortex interaction and wing vortex 
breakdown characteristics, improves the ap- 
plicability of computational fluid dynamics 
towards the design and optimization of ca- 
nard-configured aircraft. 

Status/Plans. The effects of canard deflec- 
tion, positioning, and size on wing-body aero- 
dynamics are currently under investigation. 
A study of unsteady canard motion with ap- 
plications towards aircraft control and ma- 
neuverability is planned. 

Eugene L. Tu 

Applied Computational Fluids Branch 
Ames Research Center 
(415) 604-4486 
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Figure 7-12. Fundamental Research in Vortex Interactions 


144 


a (deg) Frequency (H*) 







7-12 Fundamental Research. in ¥ortex Interactions 



ske detail with UW and unsteady surface 
pressure "m»!yg i Ti*p~mp ntsi Thp figure shows 


breakdown moves ahe ad of the tails in the 
lower left of the figure. The power spectral 
densities (PSD) are shown on the right. Sev- 
eral of the buffet pressures peak near the 



Significance. This fundamental program is 
an important topic to the high-performance 
military aircraft manufacturers because of 
the problem of structural fatigue caused by 
va^tey -fim inter action on aircraft am pa nnag es 

Hus tost will broaden the database needed to 
validate the tail buffet laws currently being 
developed. This is the first experiment to 
corralate the unsteady off-body flow-fieldchar- 
aeterfsiics, i.e., velocity fluctuations, with the 
fin surface pressures. 

Statns/Plans. The unsteady surface pres- 
sure measurements are currently being ac- 
quired through the same angle-of-attack 
sweeps. The acquisition of the flow-field mean 
and fluctuating velocities is scheduled to be- 
gin in December 1991. 





Chapter 8 


Hypersonic Aerodynamics 

The Generic Hypersonics Aerodynamics program exists to institute a fundamental, permanent 
base of technological superiority in hypersonic aero/ aerothermodynamics and to provide NASA 
with the long-term capability to achieve the goals of future NASA/DoD/DOE hypersonic vehicle 
research and development programs. It focuses on increasing the understanding of hypersonic 
aero/aerothermodynamics phenomena associated with slender, air-breathing hypersonic ve- 
hicles that use highly integrated airframe/propulsion systems and increasing the capability to 
experimentally and computationally simulate and analyze these phenomena. The program 
emphasizes both enabling and enhancing technologies related to hypersonic flight such as 
numerical simulation of turbulence and boundary layer transition, real-gas effects, rarefaction 
effects, shock/boundary layer interactions, and validation of computational methods. It explores 
flight research opportunities for technology areas that cannot be adequately addressed by 
computational methods or ground facilities. A major objective of the program is to structure cost- 
effective partnerships with industry which will enhance the rapid and efficient transfer of 
research results into the design and development of aerospace vehicles. University participation 
involves sponsored research and a refocusing and rededication of the aeronautics/astronautics 
curricula to hypersonics. 

Program Manager: Jim Moss 

OAST/RF 

Washington, DC 20546 
(202)453-2820 
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Figure 8-1, Hypersonic Shock-Wave Boundary Layer Interaction 
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8-1 Experimental Study of Hypersonic Shock-Wave/Turbulent- 
Boundary Layer Interaction Flows 


Objective. To provide accurate new hyper- 
sonic data for turbulence modeling and code 
validation. 

Approach. Design, construct, and test a 
series of simple geometries that create three- 
dimensional (3-D) hypersonic shock- wave/tur- 
bulent-boundary layer interaction flows in 
the Ames 3.5-Foot Hypersonic Wind Tunnel. 

Accomplishments. A new testbed was con- 
structed for the 3.5-foot wind tunnel which 
enables one to test a series of 2-D and 3-D flow 
fields of varying complexity. The first test 
series was successfully completed. The test 
geometry consisted of a single sharp fin whose 
shock-wave interacts with a turbulent bound- 
ary layer. Test conditions were Mach number 
8.2, Reynolds number Re - 9 x 10® and fin 
angles from 5° to 15°. The measurements 
included surface flow direction, pressure and 
heat transfer, as well as detailed Pitot pres- 
sure and yaw angle surveys of the flow fields. 


Significance. These new data represent the 
first 3-D hypersonic shock-wave/boundary 
layer interaction data suitable for turbulence 
modeling and code validation. 

Status/Plans. More complex flow fields will 
be investigated with two fins and/or compres- 
sion ramps, which will begin to look like 
hypersonic inlets. 

C. C. Horstman, Marvin Kussoy 
Fluid Dynamics Division 
Ames Research Center 
(415) 604-5950 
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Figure 8-2. Numerical Performance Estimates for a Generic Hypersonic Forebody 
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8-2 Numerical Performance Estimates for a Generic Hypersonic 

Forebody 


Objective. To examine the effects of nose 
bluntness and real-gas models on the com- 
puted flow physics andintegrated performance 
parameters of a generic hypersonic forebody 
operating at NASP-like flight conditions. 

Approach. The flow over the McDonnell 
Douglas Generic Option #2 36-inch blended- 
wing-body was computed for Mach 16 flow at 
an altitude of 125,000 ft. The parabolized 
Navier-Stokes solver, UPS, was employed with 
each of three available gas models: (1) perfect 
gas, g = 1.4, (2) equilibrium air, and (3) finite- 
rate air. Initial conditions for the space- 
marching calculations were provided by the 
time-dependent Navier-Stokes solver, TUFF, 
which possesses the same real-gas capabili- 
ties as UPS. Bluntness effects were studied at 
wind tunnel conditions (Mach 11.4, ReL= 2.904 
x 107). 

Accomplishments. Solutions were obtained 
for sharp- and blunt-nosed versions of the 
Generic Option vehicle under wind tunnel 
conditions, assuming a perfect gas. At zero 
incidence, the effect of bluntness was ob- 
served to be most significant along the lower 
centerline because of the accumulation of en- 
tropy-layer fluid into a viscous core along this 
line. The magnitude of the bluntness effect on 
centerline pressure and heat transfer distri- 
butions predicted by the computations gener- 
ally agrees well with experimental data. 

Real-gas effects were studied for the Mach 16 
flight conditions by computing forebody flow 
fields using each of the three available gas 
models in the UPS and TUFF codes. A limited 
investigation of chemistry coupling and grid 
effects was also performed. These param- 
eters were studied to evaluate their effect on 
the lower centerline flow physics and inte- 
grated performance measures such as lift, 
drag, and inlet mass flux. Integrated inlet 


performance measures were generally more 
sensitive to real-gas effects than forebody 
performance measures such as lift, drag, and 
pitching moment. 

Significance. The bluntness study provided 
confidence in the ability of the present CFD 
codes to predict surface pressures accurately 
at high Mach numbers for relatively complex 
configurations. Confident heat transfer pre- 
diction remains elusive because of grid sensi- 
tivity and boundary layer transition uncer- 
tainties. The high energy results give an 
indication of the level of importance of real- 
gas effects at NASP-like cruise conditions on 
both forebody flow-field physics andintegrated 
performance parameters. 

Status/Plans: The Option 2 geometry will be 
employed again to address issues involved in 
chemistry modeling such as loose vs. tight 
coupling and conservative vs . nonconservative 
form of the species equations. 

Scott Lawrence 

Applied Computational Fluids Branch 
Ames Research Center 
(415)604-4050 
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Figure 8-3. Forebody Pressure Contours 
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8-3 Calculation of Forebody Flow Field for a Candidate 
Aero-Space Plane Configuration 


Objective. To demonstrate the application 
of the Langley Aerothermodynamic Upwind 
Relaxation Algorithm (LAURA) code to com- 
pute the flow over the forebody of a candidate 
Aero-Space Plane configuration, and to inves- 
tigate flow-field anomalies predicted with 
another state-of-the-art flow-field code and to 
compare the results. 

Approach. The LAURA code solves the time- 
dependent, thin-layer Navier-Stokes equa- 
tions using a point-implicit scheme to con- 
verge the numerical solution. A two-block 
grid was used to calculate the blunted nose 
region and the downstream portion of the 
forebody in sequential steps. For comparison 
purposes, the flow-field grids used in this 
work were identical to those used by the other 
computational methods. The calculation was 
performed for flight conditions at Mach 15 
and laminar, perfect gas flow was assumed. 

Accomplishments. A converged solution 
over the Aero-Space Plane configuration was 
obtained which yielded complete flow-field 
quantities and surface properties (i.e., heat- 
ing and pressure) for study. Computations 
were performed on the NAS Cray-2 computer 
and required approximately 22 hours of CPU 
time. Surface properties and the nature of the 
flow field at the engine inlet (forebody end) 
have been examined and comparisons with 
the other computational results have been 
made. 

Significance. Results show that the LAURA 
code is capable of flow-field calculations for 
this type of configuration and is able to pro- 
vide needed information about flow-field qual- 
ity and surface heating. For example, the 


calculation shows pockets of inflow and out- 
flow on the lower surface and a buildup of the 
viscous boundary layer along the lower 
centerline. Also, anomalies in profiles of the 
total enthalpy at the engine inlet face, which 
were predicted in the other computational 
results, were not confirmed by the present 
calculation. 

Status/Plans. Efforts are ongoing to resolve 
the anomalies in total enthalpy profiles. 
Knowledge gained in the present work will be 
applied in a forthcoming experimental/com- 
putational study in support of the National 
Aero-Space Plane program. 

Richard A. Thompson 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4367 
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Figure 8-4. Computational/Experimental Parametric Study of 3-D Scramjet 

Inlets at Mach 10 
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8-4 Computational/Experimental Parametric Study of 3-D 

Scramjet Inlets at Mach 10 


Objective. To calibrate current computa- 
tional fluid dynamics (CFD) computer codes 
for hypersonic forebody/inlet flows using ex- 
perimental data from the Generic Option #2 
program. 

Approach. Numerical simulation of three- 
dimensional hypersonic flow past a forebody/ 
inlet configuration is performed using a state- 
of-the-art family of upwind methodology CFD 
computer codes, CFL3DE and CFL3D. Nu- 
merical results, including surface pressure 
and heat transfer along the forebody and 
inside the inlet, are compared with experi- 
mental data for the same configuration and 
freestream conditions. 

Accomplishments. The forebody/inlet cal- 
culation has been completed. The computed 
forebody pressures and heat transfer exhibit 
reasonably good agreement with the experi- 
mental data as does the computed results in 
the inlet. In some regions of the inlet, the 
agreement is quite good considering the 
complex nature of the flow. In other regions, 
the agreement is poor indicating deficiencies 
in the numerical modeling. 


Significance. The results demonstrate the 
ability of this family of CFD computer codes to 
predict forebody pressures and heat transfer 
rates with reasonably good accuracy. The 
results in the inlet are encouraging, but also 
indicate areas of needed improvement!!! the 
numerical modeling of internal flows. In 
particular, the turbulence modeling needs to 
be improved. 

Status/Plans. Comparisons of the computed 
results and experimental data will be sent to 
the Generic Option #2 report authors for their 
comments and to learn of any update in the 
experimental data. After their response, the 
results will be documented in a NASP publi- 
cation. Also, new turbulence models, better 
suited to internal flows, will be incorporated 
into this family of CFD computer codes. 

Arthur D. Dilley, George F. Switzer, 

William M. Eppard 
Langley Research Center 
(804) 864-2288 
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8-5 Technique for Hypersonic Powered Tests of Airbreathing 

Configurations 


Objective. To experimentally evaluate an 
unproven technique of using metric model 
parts to conduct simulated exhaust powered 
tests of an airbreathing vehicle in a hyper- 
sonic wind tunnel. 

Approach. The method required separate 
tests of two identically scaled models. The 
external surface force/moment contributions 
from the nose to the engine cowl trailing edge 
were determined from tests of the sting 
mounted, unpowered forebody model that was 
fully metric. The aftbody external surface 
force/moment contributions from the cowl 
trailing edge rearward were determined from 
tests of an entire model that had a metric 
aftbody. This entire model was supported by 
a strut mounted to the forebody and had a 
closed inlet and a metric aftbody. The scramjet 
exhaust was simulated by using a 7 0/30 (per- 
cent mole fraction) mixture of tetrafluoro- 
methane and argon. 

Accomplishments. Both models were tested 
in the Langley 20-Inch Mach 6 Wind Tunnel 
at a Reynolds number of 1.0 x 10®/ft. The 
appropriate forebody and aftbody force ac- 
counting procedures for this metric model 
parts technique were developed and demon- 
strated. In addition, facility operating issues 
related to conducting exhaust simulation pow- 
ered testing with a metric aftbody model were 
examined. Also, the powered metric aftbody 
model was configured with various wing inci- 
dence angles (-1.5°, -3.0°, -4.5°) and elevon 
deflection angles (0°, +10°, and +20° ) to 
examine the trim characteristics with simu- 
lated exhaust. 


Significance. This test program produced 
the first hypersonic powered force/moment 
data on a National Aero-Space Plane-type 
configuration. The results of this effort have 
shown that (1) the effect of the exhaust plume 
on the measured forebody forces/moments is 
negligible; (2) the internal inlet contribution 
to the forebody drag is less than 10%; and (3) 
the interference effect on the measured aftbody 
forces/moments as a result of the forebody 
strut mount is negligible. 

Status/Plans. Additional tests of the pow- 
ered TTD model in the 20- Inch Mach 6 Wind 
Tunnel are planned to investigate the use of 
air, as opposed to the tetrafluoromethane/ 
argon mixture, as a simulated exhaust gas. 
This effort is part of Government Work Pack- 
age 5 in support of the National Aero-Space 
Plane Program 

David W. Witte 

Applied Aerodynamics Division 
Langley Research Center 
(804) 864-5589 
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8-6 Effect of Inlet Representation on Powered Hypersonic 

Aftbody Flows 


Objective. To assess the aftbody perfor- 
mance of a powered hypersonic airbreather 
(PHA) model by computationally examining 
the effects of a flow-through inlet and a faired- 
over inlet at a variety of freestream condi- 
tions. 

Approach. One method for obtaining pow- 
ered effects on PHA concepts is to fair over the 
inlet and then route alternate gases through 
the combustor section to simulate the com- 
bustion products. However, the fairing alters 
the external flow field and possibly affects the 
aftbody local flow field and surface pressures. 
The effect of the fairing was studied using the 
General Aerodynamic Simulation Program 
(GASP) to analyze two PHA configurations. 
The flow-through inlet model provided the 
necessary inlet geometry to simulate the flow 
associated with the cowl leading edge. The 
faired-over inlet model eliminated all cowl 
leading-edge effects by covering the inlet, 
resulting in a simple expansion over the cowl. 

Accomplishments. Two-dimensional CFD 
simulations were performed on the two con- 
figurations at various freestream conditions, 
including Mach numbers 3.0 to 14.0 at Re^ = 
2.0 x 10 vft. and Reynolds numbers 0.5 to 7.0 
x lOfyft. at Moo = 6.0. These conditions yielded 
static nozzle pressure ratio (SNPR) values 
from 12 to 256 because of the different 
freestream pressures. The internal condi- 
tions were identical for all solutions, namely 
M = 1.0 at the nozzle throat with a jet total 
pressure of 30 psi. 

Comparison of Mach number contours for the 
two inlet representations is shown in the 
figure for the conditions at = 14.0 and Re co 

= 2.0 x 10°/ft. Significant differences can be 
seen in the flow fields near the cowl for the two 
configurations. The cowl leading edge on the 
flow-through inlet creates a shock followed by 
a q ioderate expansion of the flow and then a 
recompression between this expansion region 


and the cowl. For the faired-over inlet, there 
is only an expansion. Also, the size of the 
boundary layer on the flow-through inlet case 
is significantly smaller because it develops 
only from the cowl leading edge instead of on 
the entire forebody. In the aftbody region, the 
flows are quite similar, including the location 
of the cowl tr ailing-edge shock, the shear 
layer, and the weak compression within the 
shear layer. The line plot compares aftbody 
surface pressures and integrated force and 
moment components for the two PHA configu- 
rations at this freestream condition. The 
result is that the difference in integrated 
aftbody pressures due to different represen- 
tations is less than 1 %. This comparison 
remained consistent for the entire set of 
freestream conditions that were simulated. 

Significance. Within the scope of this in- 
vestigation, two-dimensional CFD results in- 
dicate that there is no significant effect on 
powered aftbody performance due to fairing 
over the inlet compared to a flow-through 
inlet. Flow-field differences exist for the two 
inlet representations, but they are isolated to 
the flow region near the cowl and appear to 
have minimal impact on aftbaody surface 
pressures. 

Status/Plans. There is a need to extend the 
CFD solutions to three-dimensional to ac- 
count for 3-D relief effects as well as proper 
expansion of the plume onto the entire span of 
the aftbody and any wing or control surfaces 
that may be subject to influence by the ex- 
haust flow. 

Lawrence D. Huebner 
Applied Aerodynamics Division 
Langley Research Center 
(804) 864-5583 
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Figure 8-7. 2-D GASP PNS Turbulent Solution vs Experiment 






8-7 Generic Scramjet Nozzle/Aftbody Studies 


Objective. To determine optimum numeri- 
cal algorithms and codes for analyzing diffu- 
sive multiple-species, hypersonic 3D flow 
fields, including comparisons with appropri- 
ate experimental data. The goal is to develop 
procedures for the analysis and design of 
generic hypersonic aftbody geometries. 

Approach. Several Euler/N avier-Stokes 
codes, employing a variety of numerical algo- 
rithms, were exercised to analyze flow about 
a powered generic scramjet nozzle/aftbody 
model. The calculated results were compared 
with measured surface (static) and flow-field 
(pitot) pressures for several single- and mul- 
tiple-species flows. 

Accomplishments. Flow-field pitot pres- 
sures were measured about a powered Ge- 
neric Scramjet Nozzle/Aftbody (GSNA) model 
in the Langley 20-Inch Mach 6 Wind Tunnel. 
The GSNA has been modeled using computa- 
tional fluid dynamics (CFD) codes ranging 
from a 2-D Euler (inviscid) code to 3-D Navier- 
Stokes (viscous) codes. An explicit Parabolized 
Navier-Stokes (PNS) code with multiple gas- 
eous species capability was modified to allow 
solid wall boundary conditions internal to the 
grid to define a cowl separating a jet flow from 
a freestream flow. Calculations inEuler mode 
were made to determine the accuracy obtain- 
able with inviscid calculations. Results from 
both 2-D and 3-D calculations were compared 
with the experimental data. The implicit, 
multiple-zone General Aerodynamic Simula- 
tion Program (GASP) code was similarly exer- 
cised and shown to be more robust than the 
explicit code for high nozzle pressure ratios. 
The jet fluid was modeled as air (perfect gas), 
a mixture of nitrogen and oxygen, and two 
mixtures of argon and Freon to simulate the 
ratio of specific heats of actual hot scramjet 
exhaust gases. GASP 2-D PNS calculations, 
including turbulence modeling (as shown in 
the figure), showed improved agreement with 
experiment over Euler solutions, both in terms 


of surface pressures and flow-field details 
around the cowl trailing edge. A 3-D explicit, 
elliptic full Navier-Stokes code was modified 
to accept block-structured grids and two mul- 
tiple-species, half-span nozzle calculations 
were performed. 

Significance. Proposed hypersonic vehicles 
will operate in flight regimes beyond the ac- 
cessibility of ground-based facilities. The 
ability to numerically analyze and design 
scramjet/aftbody components for such vehicles 
is crucial for efficient, interdisciplinary pro- 
pulsion/airframe integration. 

Status/Plans. A NASA Technical Paper has 
been written documenting the flow-field pres- 
sure measurement experiment (to be released) 
and an AIAA paper documenting the CFD 
results and data comparisons (surface and 
flow field) has been published. A full 3-D 
hypersonic vehicle will be numerically mod- 
eled with powered scramj et exhaust flow simu- 
lation and comparisons made between the 
CFD results and ongoing wind tunnel mea- 
surements. In addition, the GASP code will 
be coupled with an iterative design/optimiza- 
tion code to permit automated design of 
scramjet aftbody geometries for optimum lift 
and thrust constrained by pitching moment 
requirements. 

Kenneth E. Tatum, William J. Monta 
Applied Aerodynamics Division 
Langley Research Center 
(804) 864-5587 
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8-8 Advanced Aero-Propulsion Performance Design Tool 

(PEMACH) 


Objective. To improve the current capabil- 
ity to rapidly predict local pressure and ther- 
mal loads on hypersonic conceptual designs 
for structural analysis and stability and con- 
trol analysis. 

Approach. An engineering technique based 
on the method of characteristics (MOC) was 
chosen as a starting point for this study, 
which began last year. Although this particu- 
lar technique was originally developed for 
applications to supersonic wings and conical 
shapes, it has shown very good agreement 
with experimental and flight data up to a 
Mach number of 4. The method is being 
adapted to complete hypersonic configura- 
tions and results compared with experimen- 
tal data, Navier-Stokes (NS) calculations, and 
engineering methods currently in use. 

Accomplishments. The code, PEMACH, (a 
Practical Engineering Method for Aero-pro- 
pulsion Characteristics at Hypersonic Mach 
Numbers) can now solve the flow over the 
fuselage and wings, including capture of en- 
gine mass flow, in a quasi-3-D fashion at 
hypersonic Mach numbers. The code includes 
solution capabilities for either sharp or blunt 
leading and trailing edges, and for completely 
arbitrary planform. 

Part of the method development was calibra- 
tion to experimental results, as illustrated in 
the figure. The illustrated comparison used 
the experimental results from the National 
Aero-Space Plane (NASP) Generic Option #2 
3-D forebody-inlet model. Excellent agree- 
ment with wall pressures, both on and off 
centerline, and shock locations are illustrated 
in the figure. Comparison with CFD solutions 
verify accurate prediction of mass capture, 
and comparisons with unpowered aerody- 
namic test results from the NASP Test Tech- 
nique Demonstrator (TTD) model verify accu- 
rate prediction of C m a over a large variation 
in angle-of-attack. 


Significance. This capability brings the 
ability to quickly predict configuration loads 
suitable for structural, thermal, and stability 
and control analysis. 

Status/Plans. A preliminary version of 
PEMACH has been released and is currently 
being used for aerodynamic (lift, drag, mo- 
ments) analysis and design of hypersonic ve- 
hicles. Work is continuing to efficiently model 
the inlet and nozzle flows, to fully couple the 
boundary layer to the 3-D effects, and include 
boundary layer transition prediction. 

Suresh H. Goradia, 

VIGYAN Associates 

Abel O. Torres, Charles R. McClinton 
Hypersonic Technology Office 
Langley Research Center 
(804) 864-6253 
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Figure 8-9. Wing Glove for Pegasus Crossflow Transition Experiment 
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8-9 Wing Glove for Pegasus Crossflow Transition Experiment 


Objective. To design and fabricate a wing 
glove that provides a smooth non-ablating 
test surface for the measurement of hyper- 
sonic 3-D crossflow boundary layer transi- 
tion. 

Approach. Ablation products from the cur- 
rent thermal protection system on the Pega- 
sus delta wing will interfere with sensitive 
aerodynamic measurements like boundary 
layer transition. To minimize this interfer- 
ence, an aerodynamically shaped glove with a 
smooth nonablating surface was installed on 
the existing wing. A ceramic glove fairing was 
used to provide a smooth transition between 
the wing and glove surface. All of the flight 
instrumentation was contained within the 
glove so that the integration with Pegasus 
was simply a cable connection. 

Accomplishments. A prototype wing glove 
was fabricated to identify and demonstrate 
concepts for meeting the requirements of the 
crossflow transition experiment. The glove 
consists of a copper skin attached to the ce- 
ramic blanket known as TABI (Tailorable 
Advanced Blanket Insulation). Since the 
aerothermodynamic heating at hypersonic 
speeds is substantial, a copper skin was cho- 
sen to quickly conduct heat away from the 
leading edge and therefore reduce the ther- 
mal gradient and stresses. TABI both sup- 
ports the copper skin and insulates it from the 
composite wing. Recent work with the rigid 
ceramic insulation known as TUFI (Tough 
Uni-Piece Fibrous Insulation) shows good 
progress toward fabrication of the ceramic 
glove fairing. 


Significance. Based on this work, a FX-0 
glove will be fabricated and flown on the 
Pegasus flight, scheduled for June 1992. This 
FX-0 glove is the first of three gloves planned 
for this hypersonic flight research program. 
The goal of FX-0 is to validate the glove 
concept for measuring crossflow transition on 
a swept delta wing at hypersonic speed. Fu- 
ture gloves will be optimized for making more 
precise measurements of this phenomena. 

Status/Plans. Further work is required to 
finalize the fairing design and glove assem- 
bly. Fabrication of the FX-0 glove and its 
fairing is scheduled to begin January 1992. 

Paul Kolodziej 

Thermal Protection Materials Branch 
Ames Research Center 
(415) 604-5377 
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Chapter 9 

Aeroacoustics Research and Technology 


Aeroacoustics is concerned with the fundamental understanding, prediction, and control of noise 
and acoustic loads produced by the motion of fluids and bodies moving through the atmosphere. 
To achieve this objective, research is being conducted in three main areas: 1) Computational 
Aeroacoustics - establish CFD methodology as a discipline for acoustics technology development 
focusing on computational models of high-speed flow noise and shock/vortex interactions. 2) 
Supersonic Jets/Loads - understand/predict/reduce noise and acoustic loads generated by high- 
speed aircraft and, in particular, establish innovative concepts for supersonic jet noise suppres- 
sion and develop a detailed acoustics loads database on airframe structures. 3) Long Range 
Propagation - model atmospheric propagation effects supported by fundamental verification/ 
validation experiments. 

Program Manager: Benjamin Neumann 

OAST/RF 

Washington, DC 20546 
(202) 453-2795 
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100% thrust, standard lift configuration 


Sideline 116.3 
Centerline 116.3 


Sideline 112.3 
Centerline 112.3 


80% thrust, 30% lift increase 




Figure 9-1. Contours of Effective Perceived Noise Levels (EPNdB) for Standard 
Takeoff and Reduced Thrust/Increased Lift Takeoff for a Jet Aircraft 


168 






9-1 Jet Noise Predictions for Reduced Thrust Takeoff Study 


Objective. To demonstrate the use of the 
Aircraft Noise Prediction Program (ANOPP) 
to evaluate the effects of reduced jet thrust on 
community noise levels. Reduced jet thrust 
levels at takeoff increase the required dis- 
tance on the runway (ground roll) before lift- 
off. To compensate for the increased ground 
roll, several levels of increased lift were in- 
cluded in the study. 

Approach. Predictions were made using the 
ANOPP dual stream coannular jet noise mod- 
ule for the AST-205- 1 aircraft powered by 
GE21/J11B14A scaled engines. A power set- 
ting of 100% thrust and a standard lift con- 
figuration were used as a baseline for this 
study. Power settings varied from 100% to 
80% maximum thrust, while lift was increased 
in increments of 15% to a maximum of 60%. 
An angular rotation rate of 3 per second and 
a constant climb angle of 8° were used. Effec- 
tive Perceived Noise Level (EPNL) values 
were calculated at the FAA prescribed noise 
measuring points. Noise levels were also 
calculated on a l-by-6 mile grid of 65 observer 
positions, with each observer located 4 feet 
above the ground. Noise contours were gener- 
ated from these observer positions. 

Accomplishments. Contours for both the 
baseline case and the 80% thrust/30% in- 
creased lift case are shown in the figure. The 
reduced thrust case decreases the necessary 
ground roll distance by 646 feet and lowers 
the sideline EPNL level by 4 EPNdB from the 
maximum thrust baseline case. Although the 
centerline noise level is not of major concern, 
it too benefits by 4 EPNdB from the reduction 
in thrust and increase in lift. Additionally, as 
the figure shows, reducing thrust and in- 
creasing lift alter the shapes of the contours. 
As thrust is decreased and lift is increased, 
the contours tend to shorten in length and 
become more compact. 


Significance. The use of the ANOPP predic- 
tion code shows that reduced thrust takeoff 
procedures can be beneficial. Sideline noise 
levels are of greatest concern during takeoff. 
Since the FAA sideline noise measurement 
point location depends on the point of lift-off, 
only a slight benefit can be derived from in- 
creasing the amount of lift while keeping 
thrust at a maximum. The key to reducing jet 
noise lies in thrust reduction. 

Status/Plans. It has been shown that re- 
duced thrust, used in conjunction with in- 
creased lift to compensate for increased run- 
way distance, can reduce Effective Perceived 
Noise levels at the sideline and centerline 
FAA noise measurement locations. Better jet 
noise prediction modules being designed for 
ANOPP will allow the noise impact from a 
high-speed civil transport (HSCT) aircraft on 
community noise to be evaluated. 

Robert A. Golub 
Acoustics Division 
Langley Research Center 
(804) 864-5281 
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Figure 9-2. Scale Model of McAir 279-3C Installed for Hot-Gas Ingestion Testing (top); 
feffect of Jet Temperature on Measured Overall Sound Pressure Levels Measured at Three 

Levels on the Model (bottom) 
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9-2 ASTOVL Acoustic Loads Test 


Objective. To characterize the acoustics and 
dynamic structural loading of an advanced 
short takeofFvertical landing (ASTOVL) air- 
craft operating in ground effect and at el- 
evated jet exit temperatures. This research 
was part of a Lewis/Langley cooperative pro- 
gram, complementary to in-house research 
aimed at improving the understanding of the 
impingement acoustics of multiple heated jets. 

Approach. A 9.2% scale model of the McAir 
279-3C aircraft was tested in the Lewis 9- by 
15-Foot Low-Speed Wind Tunnel. The model 
(see accompanying figure) was built to evalu- 
ate the hot gas ingestion (HGI) characteris- 
tics of the design at temperatures up to 1000° 
F, at nozzle pressure ratios (NPR) of 4. The 
model geometry parameters varied included 
forward nozzle splay angle (from 12° inward 
to 18° outward) and Lift Improvement Device 
(LID) configuration (on, off). Acoustic data 
were obtained from a linear array of five 
freestream microphones located parallel to 
the model centerline. Model dynamic loads 
were measured with eight watereooled dy- 
namic pressure transducers located on the 
underside of the canard, wing, and aft fuse- 
lage. 

Accomplishments. Data were taken for 
over 700 combinations of model geometry and 
test conditions. The data shown illustrate the 
minimal effect of jet temperature on the noise 
level (OASPL) as measured on the canard, the 
wind, and in the freestream, while the jets 
were operating at the design condition. As 
shown by the data, levels up to 160 dB were 
observed on the canard (which was closer to 
the nozzles than the wind) and up to 140 dB in 
the freestream. At other conditions, levels as 
high as 180 dB with discrete tones 160 dB in 
amplitude were measured. Tones were gen- 
erated at a wide variety of conditions (model 
height, nozzle pressure ratios, jet tempera- 
turep) and geometries. Both free-jet tones 


and impingement tones were observed, and 
their frequencies were well predicted by cur- 
rent jet noise theories. 

Significance. Acoustic loads on ASTOVL 
configurations can potentially attain suffi- 
ciently high levels to cause structural failure. 
No satisfactory theory exists that can predict 
the acoustics of multiple heated non- 
axisymmetric jets operating in ground effect 
in an off-design condition. The acquired data 
will be valuable in understanding both the 
fundamental physics of jet-impingement and 
the configuration issues in ASTOVL aircraft 
acoustics, 

Status/Plans. This test served as Phase II of 
Lewis’ hot gas ingestion work with this de- 
sign. Additional acoustic testing is planned 
for Phase III, including the placement of trans- 
ducers on the underside of the fuselage be- 
tween the nozzles. This will allow for the 
investigation of the dynamics of fountain im- 
pingement on the model, and the propagation 
of noise from the fountain to the farfield. 

L. Kerry Mitchell 
Acoustics Division 
Langley Research Center 
(804) 864-3625 
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Figure 9-3. Aeroacoustic Loads on Thrust Vectoring Vanes of a Model F-18 High Alpha 

Research Vehicle (HARV) 
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9-3 Aeroacoustic Loads on F-18 High Alpha Research 

Vehicle (HARV) 


Objective. To provide a benchmark data- 
base on the nearfield acoustic loads generated 
by supersonic jet plumes of high-performance 
aircraft. High aeroacoustic load levels due to 
twin jet resonance have caused component 
failures on F-15 and B-1B aircraft. A model 
and full-scale static study was conducted to 
evaluate aeroacoustic loads associated with 
the thrust vectoring control system (TVCS) of 
the F-18 High Alpha Research Vehicle 
(HARV). 

Approach. A 7.4% aft-end model of the F-18 
HARV/TVCS was designed, constructed, and 
instrumented with water-cooled dynamic pres- 
sure sensors to measure aeroacoustic loads on 
the three thrust vectoring vanes. Dynamic 
loads were measured on the model for vane 
angles from stowed (-10°) to fully deployed 
(+20°). Only military power nozzles were 
studied with nozzle pressure ratio’s ranging 
from 2 to 5 and jet total temperatures from 
104° to 680°F. In a parallel effort, researchers 
at NASA Ames/Dryden measured the near- 
field acoustic characteristics of the F-18 HARV 
with TVCS installed on the twin F-404 en- 
gines and on an unmodified F-18. The diver- 
gent nozzle flaps on the F-404 engines were 
removed on the F-18 HARV, leaving essen- 
tially convergent nozzles. Near-field sensors 
were positioned near the aircraft in a similar 
fashion to one used in the model scale study. 

Accomplishments. Overall dynamic pres- 
sure loads varied considerably with vane angle 
position and nozzle pressure ratio, as shown 
in the figure. The peak pressure value of 1.5 
psi at a nozzle pressure ratio of 2.8, recorded 
for thrusting vanes at 20°, translates to an 
estimated full- scale dynamic force load on a 
vane of between 300 and 400 lbs. Stowed 
vanes (-10°) experience much smaller dynamic 
loads, 0.4 psi rms. At higher pressure ratios, 
the initial spread rate of the jet is less, and 
less jet/vane interaction occurs. These levels 
are significantly higher than previously re- 


corded on previous F-15 and B- IB outer diver- 
gent flap studies. The F-15 and B-1B loads 
were dominatedby a twin jet resonance mecha- 
nism, which is absent in the present research 
when the vanes are fully deployed. Ground 
tests of a full scale F-18 HARV at Dryden also 
confirmed high levels of near-field loads. The 
removal of the divergent nozzle flaps on the 
baseline F-18 to accommodate the TVCS re- 
sulted in a convergent nozzle configuration. 
Convergent nozzles produce significantly 
greater shock noise. At the maximum after- 
burner power setting, a condition not tested 
in the model study, an intense tone was present 
in the near-field spectrum. This tone lies near 
the predicted frequency for twin supersonic 
plume resonance — a phenomenon known to 
contribute to the high acoustic loads for the 
convergent nozzle configuration. 

Significance. The database obtained in the 
present research effort provides structural 
design guidelines for future application of a 
TVCS. The removal of the engine nozzle flaps 
and installation of the TVCS confirms prior 
predictions that significant increases in plume 
shock noise will interact with the aircraft 
structure, which could possibly reduce the 
expected life cycle for components on the F-18 
aircraft. 

Status/Plans. The results of research will be 
reported with NASA Ames/Dryden. 

John M. Seiner 
Acoustics Division 
Langley Research Center 
(804) 864-6276 
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Figure 9-4. Direct Computation of Aerodynamic Sound Generation. Figure 1: Scattering of 
Plane Sound Waves by a Compressible Vortex; Figure 2: Sound Generated by Co-Rotating 

Vortices 
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9-4 Direct Computation of Aerodynamic Sound Generation 


Objective. The prediction of the far-field 
sound produced by a turbulent flow requires 
a detailed knowledge of acoustic source terms. 
In many theories these terms are crudely 
approximated, based on empirical correla- 
tions. Computation ofboth the acoustic sources 
and far-field sound using the unsteady N avier- 
Stokes equations allows the physics of the 
sound generation process to be studied di- 
rectly. Such a knowledge will allow direct 
validation of aeroacoustic theories and help 
develop noise control strategies. 

Approach. The Navier-Stokes equations are 
solved numerically. The infinite domain is 
truncated and newly developed obliquely 
nonreflecting boundary conditions are em- 
ployed at the computational boundaries. 

Accomplishments. To address the feasibil- 
ity and accuracy of direct computation of 
aerodynamic sound generation, several model 
problems have been considered, including: 
the scattering of sound waves by a two- dimen- 
sional (2-D) compressible vortex, and sound 
generation by a 2-D co-rotating vortex pair. 
Figure 1 shows (1) the amplitude of the acous- 
tic field generated by the interaction of inci- 
dent planar waves with the 2-D compressible 
viscous vortex, and (2) a comparison of scat- 
tering amplitude directly measured in the 
computations, with a theoretical prediction 
based on an acoustic analogy. Their agree- 
ment verifies the accurate representation of 
the acoustic sources by the numerics: 
discretization does not act as a significant 
source of sound in the computations. Figure 
2 shows a preliminary computation of the far- 
field pressure fluctuations of the sound field, 
which is generated by a co-rotating vortex 
pair, and vorticity contours in the vicinity of 
the vortex pair. 


Significance. Successful simulation of the 
simple test problems mentioned above, includ- 
ing excellent agreement with theories and 
experiments, shows that accurate direct 
computation of fluid flow with sound sources 
and far-field sound, whose energy is 10 orders 
of magnitude smaller than the flow energy, is 
feasible. 

Status/Plans. The codes developed will be 
used to compute the sound generation in two 
and three dimensional turbulent shear flows. 

Sanjiva K. Lele 

Department of Aeronautics and Astronautics 
Stanford, University 
(415) 723-7721 

Parviz Moin 

Ames Research Center 

(415) 604-5127 

Tim Colonius, Brian E. Mitchell 
Department of Mechanical Engineering 
Stanford, University 
(415) 723-9602 
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Figure 9-5. Comparison of Boundary Element Method with Experimental Results 



9-5 Computational Model for Long-Range Acoustic Propagation 


Objective. To develop a method capable of 
predicting the collective influence of the im- 
portant physical features oflong-range propa- 
gation. 

Approach. Two-dimensional sound propa- 
gation over an arbitrarily shaped ridge situ- 
ated on a locally reacting infinite boundary in 
a homogeneous medium is treated using the 
Boundary Element Method. The discretiza- 
tion procedure provides a set of linear alge- 
braic equations which are solved for the un- 
knowns over the boundary surface ofthe obsta- 
cle. The solution of the boundary unknowns 
is used to compute the sound field at any point 
of interest in the domain. 

Accomplishments. A computational method 
capable of studying long-range propagation 
was developed. An effective scheme to reduce 
computational time and storage while retain- 
ing reasonable solution accuracy using cubic 
isoparametric elements, approximate polyno- 
mials, and asymptotic expansions of the 
Hankel function was implemented. An out- 
of-core equation solver for the large, fully 
populated unsymmetric algebraic equations 
produced by the Boundary Element Method 
was developed, and computed results com- 
pared with measured data. 


Significance. This research has shown that 
the Boundary Element Method can be used to 
study long-range propagation. Several inno- 
vations have made it possible to apply the 
method to frequencies and distances never 
before possible with numerical methods. Re- 
sults also show that the method compares 
well with measured results. 

Status/Plans. The method will be used to 
help analyze data from field experiments. 

Willie R. Watson 
Applied Acoustics Branch 
Langley Research Center 
(804) 864-5290 
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Figure 9-6. Quiet Struts and Sensors for Wind Tunnel Acoustic Studies 


178 







9-6 Quiet Struts and Sensors for Wind Tunnel Acoustic Studies 


Objective. To develop the technology and 
design methodology for quiet struts and sen- 
sors to be used in subsonic airflow: primarily 
closed-return wind tunnels. 

Approach. Computational fluid dynamics 
(CFD) will be used to analyze flow fields 
around struts and junctions so as to identify 
low-drag configurations, which should thereby 
have low levels of self noise. Selected shapes 
and mounting configurations will be tested in 
the Ames 7- by 10-foot wind tunnel. Acoustic 
sensors will be tested with nose cones de- 
signed to generate low levels of wind-induced 
noise. Investigations will be made of signal 
processing methods for combining signals from 
microphone pairs such that uncorrelated 
(wind) noise will be attenuated. 

Accomplishments. A 40 x 80 wind tunnel 
test of special acoustic sensors developed by 
Northrop Corporation has shown that nose 
cones shaped like Pitot probes round noses, 
and with the sensor openings several probe 
diameters downstream from the nose shoul- 
der, generate 6 dB to 10 dB less wind noise 
than do the commercial microphones with 
nose cones used at Ames. Further reductions 
in turbulence-induced noise were achieved by 
summing the signals from two sensors in the 
same probe or by summing signals from two 
sensors in adjacent probes. Uncorrelated 
sound was thereby cancelled. 


CFD simulations have identified several strut/ 
brace junctions that have well-behaved flow 
fields and are potentially quiet components 
for microphone struts. A wide range of impor- 
tant geometric parameters have been identi- 
fied for aeroacoustic testing in the 7- by 10- 
foot tunnel. 

Significance. Improved acoustic sensors 
and support struts will result in significantly 
lower background noise in closed-return wind 
tunnels, without any modification to the facil- 
ity as long as wind noise is the dominant noise 
source, as it is in the 40- by 80-foot tunnel. 
This is equivalent to adding expensive screens 
to the circuit for turbulence control. The 
quality of acoustic data acquired from low- 
noise Advanced Ducted Propellers, HSR jet 
suppressors, and similar devices will be sig- 
nificantly improved. 

Status/Plans. Quiet struts have been de- 
signed for testing in the 7- by 10- foot wind 
tunnel. Fabrication of the struts and associ- 
ated hardware will commence in the near 
future. Discussions concerning future coop- 
erative research on sensor development with 
Northrop Corporation are being planned. 

Paul T. Soderman 
Fixed-Wing Aerodynamics Branch 
Ames Research Center 
(415) 604-6675 
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Chapter 10 

High-Speed Research 


The High-Speed Research Program seeks to provide the technology answers for an environmen- 
tally safe design of a supersonic transport. In aerodynamics the technologies focus on efficiencies 
needed to provide a low noise configuration, both in climb out and cruise. The specific objectives 
of the program are to provide enabling technologies for supersonic transports to achieve in a 
practical configuration: acceptable overland sonic booms, efficient second stage climb-out for low 
noise, and supersonic laminar flow control for peak cruise efficiency. The laminar flow control 
will significantly reduce the aircraft gross weight, hence the noise signatures. The sonic boom 
element will be evaluated in FY 1992 to determine whether the promise of low boom is great 
enough to continue research. 

There is a great body of data from the previous supersonic cruise research program that ended 
in the early 1980s. However, prediction methods, optimization techniques, and experimental 
testing have improved for complex configurations. In addition, new ideas for high lift devices and 
subsonic laminar flow control are now candidates for investigation for supersonic configurations. 
These developments have given hope that supersonic aircraft efficiencies can be improved to help 
the environmental and economic realities of the coming age. Based upon systems studies these 
efficiencies translate into goals for lift-to-drag ratios of 10 in climb, 15 in transonic cruise, and 
10 in supersonic cruise. 

The High-Speed Research Program is focused on ( 1 ) CFD modeling and validation of aeroacoustic 
analysis including measurement techniques and atmospheric modeling; (2) development of 
innovative high lift devices to provide increased lift-to-drag ratios on climb out; (3) design, wind 
tunnel test and analysis of low sonic boom configurations; (4) large-scale flight evaluation of 
supersonic laminar flow control including design codes for transition prediction, and the practical 
installation of suction systems. 

Program Manager: Benjamin Neumann 

OAST/RF 

Washington, DC 20546 

(202) 453-2795 
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0.045 Scale AST-105 High Speed Civil Transport Model 
30- by 60-Foot Tunnel 

Complete Configuration, LEF=30 deg, TEF=20 deg, Q=7.0 psf 
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Figure 10-1. Effect of Forebody Length on Directional Stability 







10-1 Effect of Fuselage Forebody Fineness Ratio on HSCT High- 

Lift Directional Stabilty 


Objective. To determine the effect of in- 
creased forebody fineness ratio (length/diam- 
eter) on the directional stability characteris- 
tics of high-speed civil transport (HSCT) con- 
figurations under high-lift conditions. 

Approach. An existing model, representa- 
tive of current HSCT designs and having a 
forebody fineness ratio of 5.0, was modified to 
permit testing with alternate forebodies hav- 
ing fineness ratios of 6.5 and 8.0. Tests were 
conducted in the 30- by 60-foot tunnel to 
identify potential vehicle integration effects 
which may possibly limit the utility of HSCT 
high-lift concepts. 

Accomplishments. Results of the investiga- 
tion showed only minor effects of forebody 
fineness ratio on configuration directional sta- 
bility for the angle-of-attack range correspond- 
ing to normal low-speed operating conditions. 
For higher angles-of-attack, increasing fine- 
ness ratio was found to have a positive stabi- 
lizing effect on directional stability. 

Significance. Projected advances in HSCT 
high-lift system effectiveness will permit sig- 
nificant reductions in vehicle wing size. Since 
the fuselage is sized by passenger/volume 
requirements (which will remain unchanged) 


reductions in wing area result in correspond- 
ing increases in relative fuselage forebody 
length and consequently increases in forebody 
fineness ratio. Previous experience with mili- 
tary vehicles has indicated that high-fineness 
ratio forebodies can significantly affect ve- 
hicle directional stability characteristics. 

These results indicate that advances in HSCT 
high-lift systems may be used to accomplish 
reduction in wing size without adversely af- 
fecting directional stability characteristics. 

Status/Plans. No future work is planned in 
this area. 

E. Richard White 
Flight Dynamics Branch 
Langley Research Center 
(804)864-1147 
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Figure 10-2. Computed Streamlines for Various Fence Heights on a NACA 4412 Airfoil 
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10-2 High-Lift Systems Research for HSCT Application 


Objective. To develop advanced high-lift 
systems that will help reduce the community 
noise impact and contribute to the economic 
success of new high-speed civil transport air- 
craft. 

Approach. A combined computational and 
experimental program is being pursued. Based 
on limited small-scale testing, theoretical and 
computational analyses, and a literature 
search for suitable high-lift devices, the 
trapped vortex was selected as the most prom- 
ising technology. 

Accomplishments. Potential flow analyses 
have demonstrated the stability conditions 
required for trapping a vortex on a wing. The 
need for two fences to trap a vortex with very 
little suction at the vortex core was a signifi- 
cant discovery. Preliminary water-channel 
and wind-tunnel tests have shown the effec- 
tiveness of the two-fence configuration for 
simple geometries. These results also vali- 
dated the computations. Additional Navier- 
Stokes computations have shown that the 
vortex trapping process is not sensitive to 
Reynolds number. 

Cooperative wind tunnel tests were performed 
with Boeing Commercial Airplane Company 
at the University ofW ashington Aeronautical 
Laboratory. These tests showed that for a 
representative high-speed civil transport 
model, a double-fence configuration gener- 
ated a tighter vortex than did a single-fence. 


Significance. The theoretical and water 
ch ann el results are very encouraging. A prac- 
tical system based on the trapped vortex should 
eliminate the need for complex leading and 
trailing edge flaps, while providing sufficiently 
high lift coefficients to allow safe, controllable 
landings. This concept can generate signifi- 
cant lift coefficients at low angles-of-attack, 
reducing the tail strike problem at takeoff and 
touchdown. Reducing the lift-off angle-of- 
attack by 4° saves at least 2500 pounds in 
landing-gear weight. In addition to lower 
landing gear weight, the wing area required 
for takeoff is reduced. Overall, a significant 
aircraft weight savings can be achieved which 
will result in lower engine thrust require- 
ments and safer, quieter operation. 

Status/Plans. The computational studies 
have shown the trapped vortex to be a viable 
technique for generating high-lift at low 
angles-of-attack. A small-scale test is under 
way in the Ames 7- by 10-foot wind tunnel 
which will define the effect that various geo- 
metric parameters (fence spacing and height, 
etc.) have on the vortex trapping process. 
These results will be applied immediately in 
another cooperative test of the Boeing low- 
speed model at the Ames 7- by 10- foot wind 
tunnel. 

James C. Ross 

Fixed-Wing Aerodynamics Branch 
Ames Research Center 
(415) 604-6722 
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10-3 Integration of Flight Simulation With Aircraft Noise 

Prediction 


Objective. To provide a tool to evaluate the 
effects of high-lift aerodynamics and aircraft 
operating procedures upon high-speed civil 
transport (HSCT) community noise. 

Approach. The approach combines the flight 
simulation capabilities of the Langley Visual 
Motion Simulator (VMS) and the noise predi- 
cation capabilities of the Aircraft Noise Pre- 
diction Program (ANOPP) to predict ground 
contours for HSCT community noise evalua- 
tion. The VMS is a ground-based flight mo- 
tion simulator with six degrees-of-freedom. It 
features a variable atmosphere model and 
computer-generated out-the-window visual 
scenery. It has an advanced transport-type 
cockpit equipped with a full complement of 
aircraft flight and engine-thrust controls as 
well as flight instrumentation displays. Aero- 
dynamic control forces on the sidestick con- 
troller and rudder pedals are provided by a 
hydraulic system coupled to an analog com- 
puter. The simulator provides the pilot with 
realistic variable-feel characteristics of air- 
craft stiffness, damping, coulomb friction, 
breakout forces, and inertia. The ANOPP 
code consists of dedicated noise prediction 
modules forjet powered aircraft. This predic- 
tion program requires a detailed description 
of the aircraft position and engine operating 
conditions as a function of time. With the 
flight trajectories and engine propulsion re- 
quirements provided from the VMS, ANOPP 
is capable of predicting the ground contours of 
the noise levels associated with simulated 
takeoff, landing, or level flight. The accompa- 
nying figure conceptualizes the relationship 
between the two software-based prediction 
systems. 

Accomplishments. In the initial evaluation 
study, the VMS has been reprogrammed to 
simulate the NASA designed AST- 105 air- 
craft. This aircraft has a cruise Mach number 
of 2.6, carries 273 passengers, and has a range 
of 4500 nautical miles. It is powered by four 


Pratt & Whitney VSCE-516 engines. The 
complete set of aerodynamic and propulsion 
data available for this supersonic aircraft 
provides a validation for the simulator. The 
ANOPP code has had upgrades performed to 
the flight dynamics, atmospheric absorption 
and jet mixing noise modules for better super- 
sonic aircraft noise prediction. The VMS 
simulator has been operated for trial flights 
and the output has been utilized in the ANOPP 
code to predict contours. 

Significance. A successful high-speed civil 
transport must be environmentally compat- 
ible. Advanced operating procedures, ad- 
vanced propulsion systems, high-lift aerody- 
namic concepts, and acoustic suppression have 
all been proposed as techniques for commu- 
nity noise reduction. The effects of using 
these new technologies must be quantified to 
ensure that a proposed high-speed civil trans- 
port can meet acceptable community noise 
limits with minimum economic penalty. 

Status/Plans. A baseline community noise 
evaluation of the AST-105 configuration will 
be completed to assess the impact of the ad- 
vanced engines and advanced piloting proce- 
dures. The simulator and noise prediction 
code will be updated to simulate a more cur- 
rent proposed HSCT aircraft. 

Robert A. Golub 
Acoustics Division 
Langley Research Center 
(804) 864-5281 
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10-4 Design System for Low Sonic Boom Configurations 


Objective. To facilitate and accelerate the 
design and analysis of low sonic boom blended 
wing-body configurations. 

Approach. Existing codes were modified so 
that all input and output geometry and array 
formats were compatible. New codes were 
written and incorporated as required. One 
new code automatically adjusts the thickness 
distribution to obtain the geometry required 
by low-boom constraints. Computerized loft- 
ing was incorporated to reduce drawing board 
work. Use of a blended wing-body concept 
eliminates the tedious problem of accounting 
for the wing-fuselage gap that occurs when 
they are treated as separate components. 
Graphics codes were incorporated for rapid 
display of relevant geometries and distribu- 
tions. 

Accomplishments. Some sample designs 
have been completed. Ground-level signa- 
ture results for one case are shown in the 
accompanying figure. 


Significance. High-speed civil transport 
(HSCT) milestones for assessment of the sonic 
boom problem can be met more easily if model 
design can be accelerated. Configurations 
with acceptable sonic boom characteristics 
would open the possibility of overland super- 
sonic flight. 

Status/Plans. Experience in the use of the 
system is being gained. Some minor tailoring 
of the codes occurs as a result of this experi- 
ence. 

Raymond L. Barger, Mary S. Adams 
Fluid Mechanics Division 
Langley Research Center 
(804) 864-2315 
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Figure 10-5. Benefit of Boom Shaping 
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10-5 Sonic Boom Shaping For Reduced Loudness 


Objective. To quantify the effects of sonic 
boom signature shaping on subjective judg- 
ments of boom loudness. 

Approach. A new Langley Research Center 
sonic boom simulator was used to obtain sub- 
jective loudness judgments from a group of 
local residents of a wide range of candidate 
boom shapes. The shapes were front-shock 
minimized and covered a range of initial to 
peak overpressures (0.25 to 1.00), initial rise 
times (1 to 4 msec), and secondary rise times 
(20 to 50 msec). Duration for all booms was 
held constant at 300 msec. 

Accomplishments. Subjective loudness rat- 
ings were obtained for the range of shapes 
shown at the top of the figure. These shapes 
were characterized by an initial rapid rise to 
an intermediate overpressure (level B) fol- 
lowed by a more gradual rise to peak 
overpressure (level A). The mean loudness 
ratings, as a function of initial rise time for 
four ratios of intermediate to peak 
overpressure, are shown in the bottom figure. 
These data show that loudness, for constant 
peak overpressure, decreased with increasing 
rise time and, for a specific rise time, in- 
creased with intermediate to peak pressure 
ratio. 


Significance. These results indicate that 
significant reductions in boom loudness can 
be achieved by detailed shaping of the boom 
signatures. They also provide data that can 
be used in the selection of a preferred boom 
shape. 

Status/Plans. Future studies will introduce 
a “house filter” into the boom simulation pro- 
cess in order to assess indoor effects. Asym- 
metrical boom signatures will also be investi- 
gated to determine whether asymmetry may 
provide additional reductions in perceived 
loudness. 

Jack D. Leatherwood 
Acoustics Division 
Langley Research Center 
(804) 864-3591 
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Figure 10-6. (1) Sonic Boom Simulator (L-90-5755); (2) Subjective Loudness of Sonic Booms 
as a Function of Perceived Level (PL). 




10-6 Prediction of Subjective Response to Sonic Booms 


Objective. A noise metric that adequately 
predicts the subjective loudness of sonic booms 
that have different physical characteristics is 
needed for the assessment of sonic boom sig- 
natures predicted for candidate high-speed 
transport configurations. To address this 
need, the sonic boom simulator shown in the 
upper figure was developed to study subjec- 
tive response to sonic booms. This investiga- 
tion was aimed at assessing the effects on 
subjective loudness of the overall duration, 
rise time, overpressure, and wave shape of 
sonic booms. 

Approach. Loudness judgments of selected 
sonic boom signatures were obtained from 36 
test subjects using the method of numerical 
category scaling. Each test subject judged the 
loudness of each of 150 simulated sonic booms 
on a scale from 0 to 10. The set of 150 sonic 
booms represented different combinations of 
overall duration, rise time, overpressure, and 
wave shape. Overall duration ranged from 25 
to 425 msec with rise times of 1, 2, 4, and 
8 msec. These ranges of overall duration and 
rise time easily encompass the values pre- 
dicted for high-speed civil transports. The 
overpressures for the different wave shapes, 
N- wave and “minimized,” varied from approx- 
imately 0.6 to 1.6 psf. The loudness judg- 
ments were averaged across test subjects and 
converted into decibel-like subjective loud- 
ness levels for each of the 150 sonic booms. 

Accomplishments. Analyses of several dif- 
ferent noise metrics indicated that perceived 
level (PL) was the best predictor of subjective 
response to sonic booms. The lower figure 
shows the subjective loudness levels plotted 
against PL for each of the four rise times. 
Perceived level adequately accounted for the 
effects of rise time, overpressure, and wave 
shape. However, a small effect of duration on 
subjective response (~1 dB) over the very wide 
duration range considered was found, which 
is not accounted for by PL. 


Significance. The results of this study indi- 
cate that a noise metric exists that effectively 
predicts the subjective response to sonic booms 
as heard outdoors. Thus, a validated method 
is available for comparative assessment of the 
subjective acceptability of sonic boom signa- 
tures predicted for candidate high-speed trans- 
port configurations. 

Status/Plans. The sonic boom simulator will 
continue to support efforts to develop a high- 
speed transport having sonic boom signa- 
tures that permit overland supersonic flight. 
Relative acceptability of various proposed 
“minimized” signatures will be investigated 
in future studies. 

David A. McCurdy 
Acoustics Division 
Langley Research Center 
(804) 864-3596 
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Figure 10-7. Indoor/Outdoor Sonic Boom Simulation Facility 
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Comparison with other noise sources; July - August, 1991 
Building response; January - March, 1992 



10-7 Indoor/Outdoor Sonic Boom Simulation Facility 


Objective. To determine an acceptable level 
of sonic boom exposure for a high-speed civil 
transport, human response to booms heard 
both indoors and outdoors is to be examined. 
The first step in this study is the development 
of a suitable sonic boom simulator. 

Approach. A simulator is required which 
can provide realistic sonic boom signatures at 
the facade of a typical residential structure. 
Since the frequency content of sonic booms 
extends to very low frequencies, conventional 
loudspeaker arrangements are unable to gen- 
erate realistic sonic boom waveforms. The 
Georgia Tech Research Institute (GTRI) was 
given the task to design and fabricate an 
array of acoustic drivers that could cover the 
frequency range from 3 Hz to 4 kHz, and 
generate an N-wave having an amplitude of 
approximately 2 pounds per square foot. 

Accomplishments. A sonic boom simulator 
has been designed that consists of 25 acoustic 
horns and loudspeakers. The custom-de- 
signed, very low frequency horns are servo- 
motor driven and radiate sound over the range 
3 Hz to 30 Hz. The other, conventional loud- 
speakers radiate sound from 30 Hz to 4 kHz. 
The loudspeakers are arranged in a vertical 
array adjacent to an uninhabited house lo- 
cated at the GTRI laboratory. 

Significance. When fully operational, this 


sonic boom simulation facility will provide the 
unique capability of exposing a house to re- 
alistic sonic boom waveforms. It is known 
that human response to sonic booms heard 
indoors is quite different from outdoor re- 
sponse. This is due to a variety of factors such 
as acoustic transmission loss and noise-in- 
duced vibration and rattle. This facility will 
enable a systematic investigation to be made 
of these indoor/outdoor differences. 

Status/Plans. The simulator is currently 
operational and human response testing, both 
indoors and outdoors, will begin in fall 1991. 

Kevin P. Shepherd 
Acoustics Division 
Langley Research Center 
(804) 864-3583 
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Surveyed areas 


• Objectives: 

• Develop sonic boom response 
questionnaire 



• Provide preliminary data on 
extent of sonic boom annoyance 

Sonic boom exposure: 

• Long term SR-71, 0.5 to 1.0 psf, 

~ 1 per week 

• Exposure ceased 6 months prior 
to study 

Findings: 

• Little to moderate annoyance 

• Startle reaction frequently noted 


• Vibration frequently noted, some 
damage attributed to sonic booms 


Figure 10-8. Preliminary Sonic Boom Survey 
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10-8 Preliminary Sonic Boom Survey 


Objective. To develop a response question- 
naire for an extensive survey and to provide 
some data on the extent of annoyance to sonic 
booms in areas which have been exposed for 
many years to sonic booms with intensities 
comparable to those from proposed future 
supersonic civil transports. 

Approach. The present study developed a 
questionnaire for subsequent and more ex- 
tensive surveys through interviews in two 
communities that had experienced supersonic 
overflights ofSR-71 airplanes for several years. 
It was estimated, based on analysis of USAF 
files of supersonic flight operations, that the 
SR-71 overflights occurred on average about 
once per week and produced sonic booms with 
peak overpressures of 0.5 to 1.0 psi, although 
exposures had ceased about 6 months prior to 
the survey. A total of 23 respondents living in 
the border area of Washington and Idaho and 
22 respondents living in central Utah were 
interviewed during July 1990. Questions 
that addressed the following topics were in- 
cluded: terminology for sonic booms , extent of 
vibration, activity interference, links to previ- 
ous surveys (e.g. , Oklahoma City, 1 964), startle 
reaction, and extent of annoyance with links 
to community annoyance by conventional air- 
craft noise. 

Accomplishments. This effort has produced 
a draft questionnaire that appears to meet 
most of the goals established for this prelimi- 
nary survey. In addition, the preliminary 
survey provides a number of findings related 
to the extent of annoyance from and percep- 
tion of sonic booms in the survey areas. People 
were aware of sonic booms and on average 
gave ratings of “a little annoyed” to “moder- 
ately annoyed” and about equal to some other 
noise source. Almost all of the respondents 
reported that they had been startled by sonic 
booms and nearly half reported that their 
sleep had on occasion been disturbed. A high 
proportion of people reported that things vi- 


brated, shook, or rattled during a boom. The 
most frequently mentioned items were the 
whole house, dishes, pictures, and knick- 
knacks. About one-quarter of the respon- 
dents associated sonic booms with damage to 
their house, although only eight people in the 
entire sample thought that sonic booms were 
an issue in their community. 

Significance. This preliminary survey has 
provided information for final development of 
a questionnaire for a major sonic boom annoy- 
ance survey. In addition, the insight on the 
extent of annoyance due to long- term sonic 
boom exposure has provided valuable but 
statistically limited information as to whether 
people will accept sonic booms with intensity 
on the order of those predicted for future 
supersonic transports. 

Status/Plans. A major survey of human 
annoyance is planned to be conducted in con- 
junction with a USAF investigation of the 
extent and prediction of sonic boom exposure 
in the Nellis Military Operating Area of south- 
ern Nevada. The NASA portion of this study 
will provide estimates of the relationship be- 
tween sonic boon exposure and community 
response. 

Kevin P. Shepherd 
Acoustics Division 
Langley Research Center 
(804) 864-3583 
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Figure 10-9. Sonic Boom Propagation - Effect of Atmosphere on Boom Shape 
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10-9 Minimum Sonic Boom Rise lime Determined by Absorption 


Objective. To quantify the effects of atmo- 
spheric absorption and turbulence on sonic 
boom wave forms. This research is in support 
of the High-Speed Research Program, which 
is addressing the feasibility of supersonic flight 
over land. 

Approach. Predictions of a recently devel- 
oped theory describing the effects of molecu- 
lar absorption and viscous dissipation on sonic 
boom rise times were compared with rise 
times determined from measured sonic boom 
signatures generated by F-16 and SR-71 air- 
craft. A total of 116 sonic boom signatures 
were analyzed in this study. The Mach num- 
ber range was from a low of 1.14 to a high of 
3.0. Flight altitudes varied between 3.84 km 
and 22.3 km, thus providing a wide range of 
peak overpressure, which is the variable of 
primary importance for predicting sonic boom 
rise time. This work was carried out by Dr. A. 
Pierce of Pennsylvania State University un- 
der a grant to consider atmospheric effects on 
sonic boom wave forms. 

Accomplishments. The accompanying fig- 
ure provides comparisons of two theoretical 
rise time curves as a function of shock 
overpressure, and rise times as determined 
from measured sonic booms. The compari- 
sons show that the rise time predicted by the 
new theory is consistently less than that de- 
termined from the measured sonic boom sig- 
natures. Since atmospheric turbulence could 
not be controlled during the measurement 
period, and atmospheric turbulence is thought 
to be a major factor in sonic boom rise times, 
most of the difference is currently attributed 
to atmospheric turbulence. This is a tentative 
conclusion, however, as the two data points 
below the predicted curve on the right-hand 
figure attest to the possibility of exceptional 
cases. 


Significance. The prediction of rise time 
provided by the new theory is considerably 
better than that given by previous theory 
which neglected molecular absorption. Fur- 
ther, the new theory provides a conservative 
estimate of sonic boom rise time. 

Status/Plans. The new absorption theory 
will be integrated into currently available 
sonic boom prediction codes. Additional com- 
parisons with other sonic boom databases are 
planned. Further work on atmospheric ef- 
fects on sonic booms will concentrate on the 
development of the tools required to predict 
the effects of turbulence. 

Gerry L. McAninch 
Applied Acoustics Branch 
Langley Research Center 
(804)-864-5269 
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Figure 10-10. Supersonic Laminar Flow Control Program 
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10-10 Supersonic Laminar Flow Control Program 


Objective. To develop and validate technol- 
ogy for practical, reliable, and maintainable 
supersonic laminar flow control concepts for 
future high-speed civil transports (HSCTs). 

Approach. A balanced program involving 
NASA and industry has been structured to 
carry out the supersonic laminar flow control 
program. The program utilizes a mix of com- 
putational efforts, ground facility experiments, 
and flight testing. Advanced computational 
fluid dynamics (CFD) methods and boundary 
layer stability codes are being developed to 
guide the design and for code validation. 
Swept- wing model experiments are under way 
in low disturbance supersonic wind tunnels to 
provide data on leading-edge transition phys- 
ics. Ongoing and planned precursor flight 
tests on the F-16XL-1 and -2 will obtain de- 
sign criteria and code calibration data to 
reduce the risk for the suction panel tests on 
theF-16XL-2. Flight tests of suction panels in 
FY 1995 will achieve extensive laminar flow 
over a range of flight conditions and provide 
data for code validation and development of 
design methodologies. 

Accomplishments. The RFP for the design, 
fabrication, and installation of suction panels 
and a suction system on the F-16XL-2 was 
released on July 29. Data from the Rockwell 
suction panel on the F-16XL-1 are being used 
to calibrate CFD codes. Laminar flow to about 
8.6% chord was obtained in initial flights at 
M=1.4 and an altitude of 50,000 feet both 
with low suction and without suction. Recent 
flights have demonstrated laminar flow to 
about 25% chord with full suction for similar 
flight conditions. The F-16XL-2 supersonic 
laminar flow control feasibility study was 
completed. Technology studies with industry 
were initiated. A non-suction swept wing 
pressure model was tested in the LaRC low 
disturbance supersonic tunnel. 


Significance. Industry is preparing propos- 
als in response to the RFP. Data obtained 
from the F-16XL-1 flight tests are providing 
inputs for code calibration and transition pre- 
diction methodology, as well as improving our 
understanding of laminar flow/transition 
physics over a highly swept supersonic wing. 
Data from ground facility tests, to date, are 
verifying the design methods and will provide 
data for CFD calibration and comparison with 
flight tests. The F-16XL-2 feasibility study 
verified that achieving extensive laminar flow 
with a suction glove on the F-16XL-2 is fea- 
sible, but illustrated that special treatment of 
the wing/fuselage area would be required to 
achieve laminar flow while not producing 
stability and control problems. 

Status/Plans. Award of the F-16XL-2 con- 
tract is anticipated for March 1992. Flight 
testingofthe F-16XL-1 is scheduled for comple- 
tion at the end of calendar year 1991 . Precur- 
sor tests with the F-16XL-2 will be initiated in 
about March 1992. 

Michael C. Fischer 
Flight Applications Division 
Langley Research Center 
(804) 864-1921 
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10-11 Flow Quality Study in the Unitary Plan \Wnd Tunnel 


Objective. To determine the flow quality in 
the NASA Langley Unitary Plan Wind Tun- 
nel (UPWT) based on measurement of cone 
transition Reynolds number. 

Approach. The potential deterioration of 
the UPWT nozzle and test section wall sur- 
face finish over the past several years has led 
to a reevaluation of the tunnel flow quality 
using the Arnold Engineering Development 
Center (AEDC) 10° Transition Cone model. 
This model was last tested in the UPWT in 
1974. The present study was conducted at the 
same Mach number and Reynolds number 
test conditions as the previous test as well as 
several additional Mach numbers concen- 
trated in the region of interest for high-speed 
civil transport (HSCT) studies. 

Accomplishments. The AEDC 10° Transi- 
tion Cone model was tested in test section two 
of the UPWT at Mach numbers of 2.36, 2.86, 
3.51, and 4.60 and over a Reynolds number 
range of 1.0 x 10® to 5.0 x 10® per foot. 
Generally the current levels of transition 
Reynolds number were slightly lower than 
those obtained in the 1974 tests. This reduc- 
tion in transition Reynolds number may be a 
result of increased tunnel turbulence levels 
caused by deterioration of the nozzle/test sec- 
tion wall finish or perhaps other flow distur- 
bance producing mechanisms (e.g., turning 
vane induced disturbances). However, the 
AEDC 10° Transition Cone model surface 
finish was found to have numerous imperfec- 
tions and may have caused the relative reduc- 
tion in transition Reynolds number. 


Significance. The current emphasis on 
HSCT studies has reiterated the need to have 
well-defined wind tunnel flow quality levels. 
The AEDC 10° Transition Cone, which has 
been tested in numerous wind tunnels around 
the world, provides a common reference to 
compare flow quality between facilities. The 
present results indicate the effects of poten- 
tial deteriorations in the UPWT flow quality 
and/or the test model surface finish. 

Status/Plans. Requests have been made to 
AEDC to permit Langley Research Center to 
refinish the surface of the AEDC 10° Transi- 
tion Cone. Once refinished the model will be 
retested in the UPWT. 

Jeffrey D. Flamm 
Applied Aerodynamics Division 
Langley Research Center 
(804) 864-5955 
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Chapter 11 

Aerothermodynamics Research and Technology 


The Aerothermodynamics Research and Technology program focuses on advancing both our 
understanding of and ability to address the issues associated with high-temperature gas effects 
as they impact the aerodynamic and heating environments encountered by vehicles and 
spacecraft for both Earth and planetary missions. The problems of high Mach number flight 
throughout the continuum, transitional, and free molecular flow regimes are of primary 
importance to the program. For the Earth-to-Orbit (ETO) vehicles, which encounter Mach 
numbers as great as 25, a significant amount of dissociation and chemical nonequilibrium exists 
at the high altitudes (i.e., above 50 km). Also, because many ETO vehicle concepts must return 
to Earth for landing, the relationship between high-speed aerothermodynamic efficiency and low- 
speed flight performance must be investigated for each configuration. For more energetic 
missions that may involve probes or aeroassist space transfer vehicles (ASTVs), flight Mach 
numbers as large as 50 are encountered, and ionization, radiation, and thermochemical 
nonequilibrium can be significant to the vehicle designs. 

The Aerothermodynamics R&T program addresses the issues through emphasized research in 
four primary areas: computational tool development; experimental research and computational 
validation; facilities research and development; and configuration assessment. The computa- 
tional tools developed include not only computational fluid dynamics methods but also computa- 
tional chemistry methods which provide a unique ability to calculate high-temperature gas 
properties. The vehicle synthesis engineering tools developed are becoming industry standards 
for preliminary design and analysis. Application of aerothermodynamics technology to vehicle 
designs will result in reduced flight environment uncertainty, optimized configurations, and 
improved performance margins. 

Program Manager: Jim Moss 

OAST/RF 

Washington, DC 20546 

(202) 453-2820 
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Figure 11-1. Advancement in Radiative Transport 
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11-1 Advancements in Radiation Transport 


Objective. To develop improved methods to 
predict radiative heating over three-dimen- 
sional configurations in thermochemical 
nonequilibrium flows. 

Approach. A radiation transport code was 
developed that minimizes the spectral array 
required for adequate resolution of the radia- 
tive flux. Techniques were defined to opti- 
mize the coupling of the radiative transport to 
a flow-field code. A method was developed to 
efficiently account for 3-D radiation effects 
over complex configurations. 

Accomplishments. A nonequilibrium ra- 
diative heating prediction method, the Lan- 
gley Optimized Radiative Nonequilibrium 
(LORAN) code has been developed. Numeri- 
cal studies have been performed to determine 
the minimum spectral points required for 
adequate predictions of the radiative flux. 
The results have been compared with ground- 
based and flight measurements, and with the 
detailed Park NEQAIR code, with good agree- 
ment. Initial studies have demonstrated the 
feasibility of using the LORAN code coupled 
with the flow-field code LAURA. A Modified 
Differential Approximation (MDA) method 
has been adopted for 3-D radiative transport, 
and has been applied successfully for 
axisymmetric, blunt-body cases. 

Significance. High-temperature radiative 
transport at equilibrium conditions over 
axisymmetric blunt bodies at zero angle-of- 
attack with the one-dimensional tangent-slab 
approximation is reasonably advanced over 
earlier studies on Earth entries and the work 
on the Pioneer Venus and Project Galileo 
probes. However, for aerobraking applica- 
tions, the radiative transport is expected to be 
at thermochemical nonequilibrium conditions 
over complex 3-D configurations at non-zero 
angles-of-attack. Current advancements in 
the development of the nonequilibrium LO- 
RAN code, coupled to the LAURA flow field, 


and the development of the 3-D MDA code are 
an excellent start to predicting and interpret- 
ing Aero-Assist Flight Experiment results, 
and in the application to aerobrakes in the 
Space Exploration Initiative. 

Status/Plans. Efforts are continuing to im- 
prove the accuracy and efficiency of the re- 
cently developed 3-D radiative transport meth- 
ods. The methods will be further demon- 
strated though applications at additional flow 
conditions and vehicle configurations. 


Lin C. Hartung 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4371 


207 



VELOCITY DISTRIBUTION FOR N 


O 



Figure 11-2. Energetics of Gas-Surface Interactions in Transitional Flows at 
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11-2 Energetics of Gas-Surface Interactions in Hansitional Flows 

at Entry Velocities 


Objective. To provide fundamental infor- 
mation on the velocities and energies of mol- 
ecules incident on a surface behind a shock 
layer in transitional flows at entry velocities 
and to provide data to support experiments on 
proposed tethered satellites. 

Approach. The direct simulation Monte 
Carlo (DSMC) method was used to model the 
flow about a 1.6-m diameter sphere in the 
altitude range of 130 to 90 km for Earth entry 
at 7.5 km/s. Samples of the velocity and 
translational energy distributions were taken 
for gas molecules striking the surface at vari- 
ous points on the sphere and for each gas 
species. 

Accomplishments. Velocity distributions 
were computed for each of five gas species at 
five different altitudes as a function of the 
circumferential location on the sphere. The 
computations required the sampling of more 
than one million molecules at each condition. 
The results will be published in AIAA Paper 
91-1338. 

Significance. Proposed experiments in tran- 
sitional, entry flows will require sophisticated 
instrumentation to study the details of gas- 
surface interactions. These DSMC results 
show the wide range of incident velocities and 
energies that must be considered in designing 
such experiments. The results also demon- 
strate the high degree of translational 
nonequilibrium present behind the shock 
layer. 


Status/Plans. Additional studies to fully 
characterize the translational velocity distri- 
butions throughout the flow field are planned 
for the purpose of evaluating possible means 
of coupling DSMC and Navier-Stokes solu- 
tion methods. 


Richard G. Wilmoth, James N. Moss 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4367 

Virendra K. Dogra 
Vigyan Research Associates 
(804) 864-2957 


209 




11-3 Asynchronous Macro tasked Relaxation Strategies for the 
Solution of Viscous Hypersonic Flows 


Objective. To exploit parallel, multitasking 
capabilities of CRAY-class supercomputers to 
improve efficiency of numerical simulations 
of viscous, hypersonic flows in thermochemi- 
cal nonequilibrium. 

Approach. An upwind-biased, point-implicit 
relaxation of the governing equations was 
implemented. Relaxation factors were inde- 
pendently defined for the viscous and in viscid 
contributions to the solution. Computational 
work was divided into several, independently 
executed tasks through partitioning of the 
computational domain using CRAY 
macrotasking. The point-implicit relaxation, 
which limits the computational stencil, in 
combination with the appropriate relaxation 
factors eliminated the need for synchronized 
communication between the tasks for steady 
flow problems. 

Accomplishments. Solutions for hypersonic 
flows over aerobrakes, including the wake 
flow, have been obtained using as many as 
eight tasks executing asynchronously. Sacri- 
ficing synchronization on coarse grain ma- 
chines has shown no adverse effects on total 
CPU time required to obtain a solution, and 
reduces turnaround time for a given job. 
Adaptive partitioning provides further, mod- 
erate improvement in required CPU time by 
moving partition boundaries to concentrate 
task activity in areas where the contribution 
to the solution error is the highest. The 
accompanying figure shows the convergence 
history for each of six tasks applied simulta- 
neously and asynchronously to the solution of 
hypersonic, reacting flow over an aerobrake. 
The pattern is somewhat chaotic, because of 
the lack of synchronization and moving parti- 
tion boundaries, but the overall trend is de- 
creasing by approximately a factor of six times 
faster in turnaround time than a single task 
application. 


Significance. Simulations of moderately 
complex entry vehicles (i.e., HL-20, aerobrake- 
payload assemblies) have large memory re- 
quirements, sometimes requiring dedicated 
use of a supercomputer. The present asyn- 
chronous, macrotasking strategy enables com- 
plete utilization of all processors in the dedi- 
cated mode with speedup proportional to the 
number of processors for coarse grain ma- 
chines like the Cray-2. Furthermore, while a 
single Cray Y-MP does not have sufficient 
memory to execute these simulations, net- 
working several such machines to provide 
even faster turnaround is theoretically pos- 
sible. 

Status/Plans. These strategies will be ap- 
plied to STS and HL-20 flow-field simula- 
tions. 

Peter A. Gnoffo 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4380 
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Figure 11-4. Particle Simulation in a Multiprocessor Environment 
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11-4 Particle Simulation in a Multiprocessor Environment 


Objective. To develop a particle simulation 
code that is portable to a wide class of vector 
and multiprocessor supercomputers and ca- 
pable of addressing large-scale rarefied flow 
problems. 

Approach. Direct particle simulation meth- 
ods model flows as a collection of discrete 
particles that interact with each other and 
with boundaries through collisions. The colli- 
sion-selection rule used is highly compatible 
with the requirements for vectorization. These 
methods have been implemented on the 128 
processor Intel iPSC/860, and a dynamic do- 
main decomposition is utilized to distribute 
the simulation across processors and to bal- 
ance the load. 

Accomplishments. A new code has been 
developed that can take advantage of me- 
dium-grain parallelism through a domain 
decomposition while retaining the vector com- 
patibility of previous codes. This provides a 
capability to examine three-dimensional flow 
problems involving gas mixtures in thermo- 
chemical non-equilibrium. The code is por- 
table over a wide class of Multiple Instruction 
Multiple Data (MIMD) machine architectures. 
Speedup as additional processors are devoted 
to the simulation remains linear to the full 
128 processors of the iPSC/860 tested to date. 

Significance. With continuing interest in 
hypersonic flight, a need exists for a computa- 
tional capability to predict the aerodynamic 
and thermal environment found around ve- 
hicles such as the National Aero-space Plane 
(NASP) or the Aero- Assist Flight Experiment 
(AFE). Because of the difficulties in applying 
the continuum equations to the very low- 
density hypersonic regimes these vehicles will 
encounter, application of direct particle simu- 
lation methods is an appropriate consider- 
ation. The numerical efficiency of this ap- 


proach and the scalability to large numbers of 
processors enables larger-scale applications 
of direct particle simulations than previously 
possible. 

Status/Plans. The current version of the 
developed particle simulation code will port to 
several parallel MIMD and vector machine 
architectures. Grid adaption techniques will 
be employed to increase resolution in areas of 
the flow having larger gradients. More ad- 
vanced models for thermal and chemical re- 
laxation, now under development, will be in- 
corporated in the current code. 


Jeffrey D. McDonald 
Eloret Institute 
Aerothermodynamics Branch 
Ames Research Center 
(415) 604-1140 
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11-5 Numerical and Experimental Analyses of Small Rocket Flows 


Objective. To provide experimental and nu- 
merical data for accurate prediction of rar- 
efied flow of small, low-thrust rockets usedfor 
control of spacecraft. 

Approach. Expansion of diatomic nitrogen 
through a low thrust rocket was considered. 
A small nozzle flow provided conditions that 
pass from continuum to rarefied. To assess 
rarefaction effects, the nozzle flow was com- 
puted with both a traditional continuum CFD 
technique and a stochastic particle simula- 
tion method. The plume was also calculated 
with the particle method. The computations 
were compared with new experimental data 
for Pitot pressure and flow angle. 

Accomplishments. Comparison of pitot pres- 
sure at the nozzle exit reveals that the par- 
ticle method predictions agree very well with 
the experimental data, whereas the continuum 
solutions overpredict the pressure by 25%. 
The particle method solutions continue to 
agree with experimental data further out into 
the expansion plume. The numerical efficiency 
of the particle method has been improved 
substantially by restructuring the algorithms 
to suit the architecture of vector supercom- 
puters. A typical run of the particle method 
code required 3 CPU hours and 24 Mwords of 
central memory on a Cray Y-MP. 

Significance. The differences in the con- 
tinuum and particle method solutions will 
increase as the plume expands. The present 
study reveals that significant errors will be 
incurred through the prediction of such flows 
with continuum CFD. To obtain accurate so- 
lutions of small plumes, and hence accurate 
assessment of the interaction between the 
plume and the spacecraft, it is necessary to 
employ the particle simulation technique. 


Status/Plans. The particle method computa- 
tions will be extended further in the trans- 
verse direction and into the backflow region 
behind the rocket. Furthermore, solutions 
will be generated through the same nozzle for 
the expansion of argon for comparison with 
additional experimental measurements. 

Iain D. Boyd 
Eloret Institute 
Aerothermodynamics Branch 
Ames Research Center 
(415) 604-4907 

Paul F. Penko 
Lewis Research Center 
(216) 433-2404 
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11-6 Hypersonic Rarefied Flow About a Compression Comer 


Objective. To numerically simulate the phe- 
nomena of shock/boundary layer interactions 
and compare the numerical results with the 
experimentally measured data obtained in 
the DLR wind tunnels at Gottingen, Ger- 
many. 

Approach. The focus is on the merged-layer 
interactions that occur for hypersonic (Moo * 
20) nitrogen flow over a two-dimensional model 
consisting of a flat plate followed by a com- 
pression ramp. The ramp corner is located 
71.4 mm from the plate leading edge, and 
ramp angles of 0°, 15°, 25°, and 35° are consid- 
ered. Freestream Reynolds number based on 
model length ranges from 6,000 to 39,000. 
The calculations are made with a 2-D version 
of the direct simulation Monte Carlo (DSMC) 
method. 

Accomplishments. Calculations have been 
made for four test conditions using the DSMC 
method. Results of these calculations have 
been compared with various experimental 
measurements to provide both quantitative 
and qualitative comparisons. The experi- 
mental and computational findings agree on 
the combination of flow conditions and ramp 
angles necessary to produce separation. 


Significance. The capability to accurately 
calculate the phenomena associated with 
shock/boundary layer interactions is critical 
for the design of future space transportation 
vehicles. The present study gives added con- 
fidence in the ability of the DSMC method to 
simulate such flows for hypersonic rarefied 
conditions. 

Status/Plans. Results of the study have 
been published in AIAA Paper 91-1313. 

James N. Moss, Joseph M. Price 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4379 

Ch.-H. Chun 

DLR Institute for Experimental 
Fluid Mechanics 
Gottingen, Germany 
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Figure 11-7. Rarefied Flow Regime Hypersonic Waverider Configuration, Mach No. 25 


11-7 Advancement in Developing an Efficient 3-D DSMC Code 


Objective. To significantly improve our 
present capability to calculate rarefied flows 
over complex configurations at low altitudes. 

Approach. The capability of a new direct 
simulation Monte Carlo (DSMC) code, de- 
noted as F3, was further developed, demon- 
strated, and validated and flow fields were 
simulated in the rarefied flow regime about 
complex configurations in three dimensions. 

Accomplishments. Significant code im- 
provements have been made in the areas of 
surface definition and grid generation for code 
setup, graphical displays for diagnosis proce- 
dures, and post-processing and analyzing re- 
sults, and in data compressing techniques to 
reduce memory size. F3 was used to simulate 
the flow field over a delta wing at angles-of- 
attack from 0° to 30° for a Mach 20 flow at a 
Knudsen number of 0.016. The aerodynamic 
coefficients were in excellent agreement with 
ground-based experimental data, thus pro- 
viding partial validation of the code. Solu- 
tions had been obtained over a hypersonic 
waverider configuration proposed by the Uni- 
versity of Maryland that was optimized for 
Mach 25 at an altitude of 90 km. Also, solu- 
tions were obtained using an approximate 
bridging technique to correct for rarefied flow 
effects at higher altitudes. Finally, solutions 
were obtained with the current F3 DSMC 
code at altitudes from 97 to 145 km. The F3 
lift-to-drag ratios were much lower than those 
predicted by the University of Maryland. 


Significance. The computational efficiency, 
code setup procedures, and memory manage- 
ment of the F3 code are far superior to our 
earlier 3-D DSMC code, with no loss in accu- 
racy. Our present simulation for the delta 
wing with the F3 code was set up in a few days 
and required 1 to 2 days on a 32-bit worksta- 
tion as contrasted with several months to 
setup and one Cray CPU day with the earlier 
code. The present results with the F3 code 
show that the high values of lift-to-drag pre- 
dicted for a waverider concept are signifi- 
cantly reduced at high-altitude, rarefied flow 
conditions. 

Status/Plans. Future improvements to the 
code will include adaptive gridding, a vari- 
able inner mesh resolution, improved memory 
management, and optimization for vector and 
parallel computing. Code will be demon- 
strated for selected 3-D configurations. 

Didier F. G. Rault 
Aerothermodynamics Branch 
Langley Research Center 
(804)864-4388 
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Figure 11-8. Shuttle Orbiter Pressure 
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11-8 Flow-Field Computations for Shuttle Orbiter 


Objective. To obtain flow-field solutions for 
the shuttle orbiter adequate for comparison 
with the flight data from the Shuttle Infrared 
Leeside Temperature Sensor (SILTS) experi- 
ment. 

Approach. Improved gridding techniques 
were used to define the surface and volume 
grids. The Langley Aerothermodynamic Up- 
wind Relaxation Algorithm (LAURA) code 
was used to obtain the flow-field solutions. 

Accomplishments. A grid and flow-field 
solution has been obtained for inviscid flow at 
a 40° angle-of-attack for Mach 7.4 flow, a wind 
tunnel case. The computed surface pressures 
are in excellent agreement with the experi- 
mental data on both the windward and lee- 
ward surfaces. The results identified a small 
region of non-smooth grid surface near the aft 
end of the vehicle. This small grid problem 
has been resolved, and a grid for a viscous 
solution has been completed. 

Significance. The first computational solu- 
tion has been obtained for the complete or- 
biter configuration at a hypersonic flow condi- 
tion. The analysis of the results for the wind 
tunnel case provides confidence that the grid 
technique and the LAURA code can provide 
excellent results; thus, we can proceed with 
more computer-intensive cases and analyze 
viscous flow with nonequilibrium chemistry 
at actual shuttle flight conditions. This solu- 
tion also adds significantly to our database of 
computational results for complex configura- 
tions such as the HL-20, National Aero-Space 
Plane, and the Aero-Assist Flight Experi- 
ment vehicle. 


StatualPlans. Viscous flow solutions with 
non equilibrium chemistry will be obtained at 
actual shuttle flight conditions. The compu- 
tational results will be compared with SILTS 
flight data and data from the Developmental 
Flight Instrumentation flights. 

William L. Kleb, K James Weilmuenster 
Aerothermodynamics Branch 
Langley Research Center 
(804) 864-4364 
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Figures 11-9. Approximate Flow Field and Heating Analysis 
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11-9 Approximate Heating Analysis 


Objective. To develop rapid but reliable 
approximate flow-field analysis capability for 
application in conceptual design studies and 
preliminary aerothermal design. 

Approach. Boundary layer techniques were 
incorporated in a versatile invisdd method 
for heat-transfer predictions over 3-D sur- 
faces at various angles-of-attack. 

Accomplishments. An inviscid method has 
been developed incorporating the normal pres- 
sure approximation of Maslen. The method 
has been demonstrated to yield rapid but 
reliable inviscid solutions over axisymmetric 
and 3-D surfaces at angle-of-attack. The code 
requires 10 seconds of computer time on a 
workstation to obtain a solution over a 3-D 
blunt nose. Boundary layer approximations 
have been included for axisymmetric vehicles 
at angle-of-attack for (1) perfect gas and equi- 
librium chemistry and (2) laminar and turbu- 
lent heating calculations. 


Significance. Successful demonstration of 
the approximate flow-field technique for 3-D 
configurations will significantly advance the 
useability of engineering methods. The analy- 
sis would provide the aerothermal commu- 
nity with a unique approximate method for 
computing the heating rates over 3-D blunt 
noses needed in thermal design studies. Cur- 
rently, more detailed methods require many 
hours of computation. 

Status/Plans. The boundary layer methods 
are presently being included in the inviscid 
method for application to 3-D vehicles. Com- 
parison of results will be made with a detailed 
Navier-Stokes method. 

Christopher J. Riley 
Aerothermodynamic Branch 
Langley Research Center 
(804) 864-4387 
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11-10 Comparison of Heating Rate Calculations with Experimen- 
tal Data on a Modified Shuttle Orbiter at Mach 6 


Objective. To provide detailed information 
on the ability of an “engineering code” based 
on combined in viscid and boundary theory 
AA3DBL/HALIS (^asymmetric Analog for 
Xbree-Dimensional Boundary Layer/Kigh 
ALpha Inviscid Solution) and a “benchmark” 
Navier-Stokes code LAURA (Langley 
Aerothermodynamic Upwind Relaxation Al- 
gorithm) to predict heating rates from the 
windward symmetry plane and in the vicinity 
of wing leading edges on winged vehicles at 
large angles-of-attack. 

Approach. Heating rate calculations from 
AA3DBIVHALIS and LAURA were compared 
with experimental data obtained at Mach 6 on 
a modified Shuttle Orbiter model. This is the 
first step in a process designed to validate the 
codes for application to winged vehicles in the 
flight environment. 

Accomplishments. Comparisons have been 
made at 30° and 40° angle-of-attack for condi- 
tions where the flow over the model was 
completely laminar. Both codes compare well 
with the experimental data along the wind- 


ward symmetry plane, in a circumferential 
direction around the model, and in the high 
heating region along the fuselage or wing 
leading edge. 

Significance. These results give confidence 
that both AA3DBL/HALIS and LAURA can 
be used to predict more accurately the heat- 
ing environment of future entry vehicles and 
should lead to less conservatism in design. In 
addition, since AA3DBL/HALIS require only 
1% of the time required for similar LAURA 
calculations, it should be very useful for per- 
forming parametric studies. 

Status/Plans. AA3DBL/HALIS and LAURA 
will be compared with wind tunnel data at 
other flow conditions and with available flight 
data. 

H. Harris Hamilton II, Francis A. Greene, 

K. James Weilmuenster 
Aerothermodynamics Branch 
Langley Research Center 
(804)864-4365 
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Figure 11- 11. Computed Meridional Pressure Contours from Mach 5.04 Flow of Premixed 
Hydrogen- Air Over a Blunt Body Using a 181x127 Grid (4 Frames of an Oscillation Cycle). 
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11-11 Numerical Simulation of Unsteady Shock-Induced 

Combustion 


Objective. To develop a computational 
method that will lead to the simulation and 
understanding of the unsteady shock-induced 
combustion observed in ballistic range ex- 
periments. 

Approach. Ballistic range experiments pro- 
duce a periodic, unsteady flow when a blunt 
projectile is fired supersonically into a 
premixed hydrogen-air mixture (Frame 1 in 
the figure). This unsteady behavior is caused 
by the combustion of the mixture, which has 
been heated by the bow shock in front of the 
projectile. Observations show that the heat 
release due to combustion does not occur im- 
mediately behind the bow shock but is de- 
layed by the finite rate of the chemical reac- 
tions in the high-speed flow. The region 
between the shock and the heat release is 
called the induction zone. There is little or no 
heat release in this zone, but the concentra- 
tions of chemical species important for igni- 
tion increase exponentially. Previous work 
has shown that the chemistry in the induction 
zone must be modeled accurately to simulate 
these flows. The grid spacing requirements 
for sufficient accuracy are so stringent that 
simulating these unsteady flows has not been 
feasible. A trans-formation of the species 
continuity equations has been found to reduce 
grid requirements and make such simula- 
tions practical. This trans-formation is com- 
patible with current numerical methods. 

Accomplishments. This work represents 
the first known numerical simulation to in- 
corporate a detailed chemical reaction model 
capable of simulating the unsteady shock- 
induced combustion phenomenon seen in bal- 
listic range experiments (Frames 2 and 3). A 
new formulation of the species continuity equa- 
tions has been developed and incorporated in 
a proven numerical scheme. This new formu- 


lation increases the accuracy of numerical 
calculations in flows with exponentially vary- 
ing species concentrations. 

Significance. The physical mechanism caus- 
ing unsteady combustion in supersonic ballis- 
tic range experiments is now better under- 
stood. This will lead to improved simulations 
of combustion processes for hypersonic sys- 
tems. 

Status/Plans: The current formulation will 
be applied to other problems in which shock- 
induced combustion or ignition is important. 
Examples are flows in shock-tunnels or ram- 
accelerators. Application of this work to other 
types of flows with exponentially changing 
species concentrations, such as ionization 
processes, may also be valuable. 


Gregory J. Wilson, Myles A. Sussman 
Aerothermodynamics Branch 
Ames Research Center 
(415)604-4228 
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11-12 Aerodynamic/Aerothermodynamic Characteristics 
of HL-20/HL-20 A/B Lifting Body Configurations 


Objective. To experimentally determine 
aerodynamic/aerothermodynamic character- 
istics of an HL-20 lifting body configuration at 
hypersonic to subsonic flow conditions (i. e., 
entry to landing), and to evolve optimum 
lifting body configurations from aerodynamic 
perspective via modification of aerolines and 
control surfaces. 

Approach. Force and moment measure- 
ments were performed on the HL-20 configu- 
ration at subsonic (Low-Turbulence Pressure 
Tunnel [LTPT], 7- by 10-ft. Tunnel), transonic 
(Calspan Transonic Tunnel, 8-foot Transonic 
Tunnel), supersonic Unitary Plan Wind Tun- 
nel [UPWT], and hypersonic flow conditions 
(HFC) to determine the effects of compress- 
ibility and viscosity and to simulate real-gas 
effects. Thermal mappings and flow visual- 
izations were performed at hypersonic condi- 
tions to determine aerothermodynamic char- 
acteristics. HL-20 aerolines, including con- 
trol surfaces, were modified to enhance aero- 
dynamic performance with no sacrifice in 
aerothermodynamic performance. 

Accomplishments. HL-20, HL-20A series 
(four different configurations having a stream- 
lined upper surface including canopy, wind- 
ward surface chamber, and reduced base area) 
configurations were tested extensively at sub- 
sonic, supersonic, and hypersonic conditions. 
Tests of the HL-20B series (second iteration; 
flat-like upper surface and wings faired into 
body) and HL-20/Titan 4 ascent configuration 
were initiated in the MSFC 14-inch tunnel. 

Significance. Unprecedented ground-based 
data have been established for HL-20/HL- 
20A lifting body configurations, clearly sur- 
passing previous phase A studies. Aerody- 
namic characteristics were determined from 
Mach 20 to 0.1 (entry to landing). Aerody- 
namic performance of HL-20 at subsonic 
speeds has improved significantly via modifi- 
cation! of the wing from cylinder-slab to airfoil 


(lift-to-drag ratio increased approximately 
20%); at subsonic conditions, the effect of 
Reynolds number vanished for values in ex- 
cess of 8 million per foot in the LTPT. Aerody- 
namic performance of HL-20s free of real-gas 
effects primarily because of flat windward 
surface (unlike the shuttle orbiter, which ex- 
perienced significant real -gas effects on pitch- 
ing moment). The HL-20A concept provides 
improved subsonic aerodynamic performance, 
but must reenter at high incidence to match 
lift/drag (i.e., cross range) of the HL-20 and 
the will experience significantly higher heat- 
ingin the nose region as a result of steamlining 
of the nose. 

Status/Plans. We plan to complete wind 
tunnel testing of HL-20, HL-20A, and HL- 
20B concepts; reduce and analyze data; and 
disseminate results. Will initiate CFD code 
calibration study at hypersonic conditions for 
HL-20 configuration. 

Bernard Spencer, Jr., George M. Ware, N. 
Ronald Merski 

Experimental Hypersonics Branch 
Langley Research Center 
(804) 864-5245 
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Chapter 12 
Aerobraking 

Aerobraking can provide a substantial increase in payload or decrease in initial mission mass 
for missions that require a trajectory deceleration in the vicinity of an atmosphere-bearing 
planet. The objective of the Aerobraking Technology program is to develop and validate key 
supporting technologies which allow the effective use of aerobraking for maimed lunar missions, 
robotic and manned round trip missions to Mars, and other planetary exploration missions. The 
term aerobraking in this context includes 1) capture into a planetary orbit, 2) reduction of orbit 
altitude, 3) direct entry from a hyperbolic trajectory, and 4) entry from orbit. The key supporting 
technologies include aerothermodynamic modeling; high temperature thermal protection sys- 
tems; adaptive guidance, navigation, and control; and lightweight structures that minimize or 
eliminate in-space construction and/or assembly. Aerobraking systems are strongly affected by 
each of these supporting technologies. System analysis that addresses missions, vehicle 
concepts, and operations is used to integrate the discipline areas, to assess technology require- 
ments and trades, and to ensure that system requirements are satisfied. 

Program Manager: Jim Moss 

OAST/RF 

Washington, DC 20546 
(202)453-2820 
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Figure 12-1. Additional Post-Aerocapture AV Requirements Resulting From Multiple 

Off-Nominal Effects 
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12-1 Effect of Atmospheric and Aerodynamic Uncertainties on 
Manned-Mars Aerobrake Feasibility 


Objective. To develop a state-of-the-art Mara 
aerobraking guidance algorithm capable of 
managing off-nominal atmospheric and aero- 
dynamic conditions. 

Approach. A predictor-corrector guidance 
algorithm with deceleration feedback was 
developed and implemented into a three-de- 
gree-of-freedom simulation to assess the im- 
pact of uncertainties in the Mars atmospheric 
density profile and vehicle aerodynamics on 
mission feasibility. Additionally, because the 
mission is manned, a load-relief capability 
was incorporated such that the maximum 
sensed deceleration did not exceed 5 g’s (1 g is 
the acceleration of a body at rest on the Earth’s 
surface). 

Accomplishment. The predictor-corrector 
algorithm has demonstrated the capability to 
accurately guide a Mars aerobraking vehicle 
through a wide range of off-nominal condi- 
tions (while providing load-relief) to the proper 
exit orbit. The off-nominal conditions simu- 
lated include errors in the atmospheric den- 
sity profile (both vertical and horizontal struc- 
ture) as well as a mispredicted trim angle-of- 
attack. The figure shows the amount of addi- 
tional AV required to insert the vehicle into 
the proper parking orbit as a result of various 
combinations of off-nominal effects. The cases 
illustrated in this figure are representative of 
several worst-case scenarios; however, in each 
case, the mission was completed successfully. 
Note that a contingency AV budget on the 
order of the nominal AV (75 m/sec) would be 
more than adequate to account for post- 
aerocapture orbital trim after atmospheric 
flight through an off-nominal environment. 
In comparison to the Mara departure propul- 
sive requirements (on the order of 2-3 km/ 
sec), this is a relatively minor increase in the 


total AV budget. The algorithm was also used 
in support of the Mars Atmospheric Knowl- 
edge Working Group (MAKWG) to help char- 
acterize the significance of off-nominal atmo- 
spheric conditions. 

Significance. With the use of an adaptive 
guidance algorithm, atmospheric and aerody- 
namic uncertainties can be adequately man- 
aged and do not pose mission failure concerns 
that would preclude the use of aerobraking in 
a manned Mars mission. 

Status/Plans. The results of this study have 
been documented in AIAA Papers 91-0058 
and 91-2873 as well as the final report of the 
MAKWG. Further refinement of this analysis 
includes incorporation of the algorithm into a 
six-degree-of-freedom Monte Carlo simula- 
tion to better understand the coupled effect of 
various off-nominal conditions and identify 
the vehicle control implications. 

R. D. Braun, R. W. Powell 
Vehicle Analysis Branch 
Langley Research Center 
(804) 864-4507 
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MANNED VEHICLE 



Figure 12-2. Heat Shield Erosion in a Dusty Martian Atmosphere 


234 








12-2 Heat Shield Erosion in a Dusty MartianAtmosphere 


Objective. To determine the extent of heat 
shield erosion from dust particle impacts that 
occur on a large, blunt aerobrake configura- 
tion during aerocapture maneuvers at Mars 
and the Mars Environmental Survey probes 
entering the atmosphere along ballistic flight 
paths. 

Approach. An explicit Navier-Stokes code 
was used to compute the flow field about the 
vehicle shapes shown in Figures la and 2a. 
Calculations using a particle size distribution 
ranging from 1 to 19 pm were made. The 
paths, deceleration, heating, and sublimation 
of silicate (SiC>2) dust spherules which enter 
the forebody shock-layer at various locations 
were calculated and traced to surface impact. 
Two heat shield materials were considered: a 
glass surface that is representative of Shuttle 
ceramic tiles and a low-density ablator, 
AVCOAT, that was used on the Apollo cap- 
sule. 

Accomplishments. The dustless flow field 
was computed at a number of points along 
several entry trajectories for the manned ve- 
hicle and the MESUR probe. Subsequently, 
dust particles were introduced into the shock 
layer and their trajectories were traced. The 
approximate surface erosion is shown in Fig- 
ures lb and lc for the manned vehicle as a 
function of radial distance from the stagna- 
tion point. For a vehicle with m/CDA = 200 kg/ 
m 2 , over 9 mm of glass can be destroyed, while 
for ablating AVCOAT the erosion exceeds 10 
mm in the stagnation region. 

This research was extended to examine the 
heat shield erosion caused by atmospheric 
dust for small, unmanned, probes entering 
the Martian atmosphere along ballistic flight 
paths. Heat shield erosion, energy and fre- 


quency of particle impacts were calculated at 
several points along atmospheric trajectories 
with entry velocities of 7 and 9 km/sec and 
vehicles with ballistic coefficients of 30 (see 
Figures. 2b, 2c) and 40 kg/m 2 . 

Significance. The presently available ero- 
sion results for the manned vehicle indicate 
that the design of the heat shield can be 
significantly affected by atmospheric dust. 
For example, it appears that glass-coated 
tiles are unsuitable heat shield materials for 
high-speed Mars entries unless coated with a 
protective material. For AVCOAT, the TPS 
mass may have to be increased by 25% to 29%, 
or about 1% of the vehicle mass, for entry at 
8.6 km/sec. For the MESUR probe the results 
clearly show that the surface insulation glassy 
tile, which has a 0.35 mm glass coating, re- 
quires a protective layer of TUFI (matrix of 
silica and alumina fibers). For the AVCOAT 
ablator, the heat shield mass must be in- 
creased by about 30% to account for erosion 
during or shortly after a dust storm. There- 
fore, the vehicle’s heat shield must be de- 
signed to account for dust particle impact 
damage. 

Status/Plans. The previous research will be 
extended to implement a coupled flow model 
to assess gas-solid chemical coupling and ex- 
amine the effect of micron-sized particles on 
the vehicle heat shield for different insulation 
materials. The latest results will be reported 
at the AIAA Aerospace Science Meeting in 
January 1992 . 

Periklis Papadopoulos, Michael Tauber 
Aerothermodynamics Branch 
Ames Research Center 
(415) 604-1146 
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Figure 12-3. Aerobraked Lunar Transfer Vehicle (LTV) Cost and Operations Study 





12-3 Aerobraked Lunar Transfer Vehicle Cost and 

Operations Study 


Objectives. To use cost and risk analyses to 
compare all-propulsive versus aerobraked 
concepts for a lunar transfer vehicle (LTV) 
and to use resultant cost and risk knowledge 
to define technology development roadmaps. 

Approach. Develop life cycle cost (LCC) 
models with sufficient depth of detail in all 
phases of LTV operations and in technology 
development programs to capture differences 
between aerobraked and all-propulsive LTV 
concepts. Use risk analysis to identify critical 
technology test requirements and risk-reduc- 
ing design changes. 

Accomplishments. Relative life cycle cost 
estimates show a 7%-10 % (inside uncertainty 
band) cost advantage for an aerobraked LTV 
when compared with an all-propulsive con- 
cept based on the 90-Day Report Option 5 
lunar mission (25 years, one mission per year, 
originating in low-Earth orbit). Cost savings 
stem from reduced Earth-to-orbit launch costs, 
reduced lunar orbit propellant boiloff, reduced 
gravity losses, and fewer on-orbit operations 
tasks (marginal advantage; assumes three- 
piece erectable aerobrake). Cost increases 
stem from Design Development Test & Evalu- 
ation (DDT&E) programs, aeropass software 
production and maintenance, and mission 
operations. Test and production facilities and 
space station (or free flyer) support facilities 
have been considered. Cost and risk analysis 
results have led to several design modifica- 
tions (relative to the baselined MSFC con- 
cept), and technology development programs 
have been defined (for costing purposes) in 
major discipline areas. 


Significance. Viability of the aerobrake for 
the lunar mission has been examined in terms 
of life cycle costs. To date, estimates had 
shown that the mass savings from aerohrakes 
were significant, but the question remained 
whether resources invested into aerobrake 
development programs exceeded these sav- 
ings. After in-depth estimates of major 
aerobrake DDT&E and operations issues, 
aerobrake options have been shown to be 
competitive with all-propulsive approaches 
for the lunar missions. New, larger launch 
vehicles or fewer missions could alter the 
results. 

Status/Plans. Additional analysis and cost- 
ing studies will address single-use ablator 
designs and aerobrake disposal, and identifi- 
cation oflunar aerobrake development spinoffs 
for future Mars missions will be made. 

Lawrence F. Rowell 
Space Systems Division 
Langley Research Center 
(804) 864-44502 
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V E = 8.6 km/sec, 5g limit 



Figure 12-4. Mass Fraction of Aerobrakes for Manned Mars Mission 


238 


12-4 Aerobrake Design Studies for Manned Mars Missions 


Objective. To systematically access the im- 
pact of aerobrake 4D and ballistic coefficient 
on aerobrake mass fraction and mission tra- 
jectory and verify the advantages of using 
aerobrakes over using chemical propellants 
for Mars entry. 

Approach. The vehicles’ mass fractions that 
must be devoted to the aerobrakes, including 
the heat shields, have been computed for a 
high-speed manned Mars entry of 8.6 km/sec 
and with an Earth 5 g deceleration limit. 
Blunt, low lift-drag ratio (L/D), configura- 
tions with ballistic coefficients (m/CjxA) of 
100 and 200 kg/m^ were studied. In addition, 
a delta-winged vehicle, with a L/D of 0.85 and 
a ballistic coefficient of 375 kg/m2, was stud- 
ied. The convective and equilibrium radiative 
heating and the pressures were computed at 
selected body locations for the Martian atmo- 
spheric gases along the overshoot and under- 
shoot trajectories. Bothinsulative, heat shields 
and ablators were considered. Polyimide/ 
graphite (PI/Gr) primary structural members 
were assumed for structural analyses. How- 
ever, the basic shell structure consisted of an 
aluminum honeycomb sandwich core between 
PI/Gr face sheets; the shell was stiffened by I- 
beam frames. Multi-variable optimization was 
used to minimize the structural weights. 

Accomplishments. After adding heat shield- 
ing and optimizing the structure, the 
aerobrakes’ mass fractions (defined as heat 
shield plus aerobrake mass divided by the 
vehicle’s total mass at entry) varied from 
about 15% to 13% for ballistic coefficients of 
100 and 200 kg/m2, respectively, for the blunt 
shapes and was slightly under 17% for the 
winged vehicle (see figure). The higher mass 
fraction for the winged-vehicle resulted from 
its much higher ballistic coefficient which 
created more intense heating and a higher 
mass fraction for thermal protection. How- 
ever, the increase in L/D from 0.3 to 0.85 for 
the winged vehicle nearly doubles the entry 


flightpath angle corridor width and provides 
an order of magnitude increase in cross-range 
capability during descent from orbit. Also, the 
winged vehicle can be launched into Earth 
orbit fully assembled. The blunt shapes will 
probably have to be assembled in orbit, which 
is complex and potentially hazardous. 

Significance. The aerobrakes’ mass frac- 
tions range from less than to slightly over the 
15% value that is considered to make 
aerobraking indisputably superior to chemi- 
cal propulsive braking. In fact, using the best 
current chemical propellants yielded a pro- 
pulsion system mass fraction that exceeded 
the aerobrakes’ values by factors of 4 to 5. 

Future Plans. The results of the study were 
reported at an AIAA meeting. In addition, the 
analysis procedures that were developed will 
be applied to minimize the aerobrake masses 
for the unmanned probe/landers being stud- 
ied for the Mars Environmental Survey 
(MESUR) missions. 

M. Tauber, M. Chargin, W. Henline, 

K. R. Hamm, Jr., H. Miura 
Ames Research Center 
(415) 604-6086 

A. Chiu and L. Yang 
Sterling Software 
Palo Alto, CA 
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12-5 Heating Rates forAeroassist Flight Experiment 


Objective. To define the heating rates and 
flow-field characteristics over the Aeroassist 
Flight Experiment (AFE) vehicle for the pur- 
pose of vehicle design, experiments’ design, 
and scientific pre-flight predictions. 

Approach. The convective and radiative 
heating rates and the flow-field characteris- 
tics were calculated with several different 
computational codes at flight conditions from 
currently considered trajectories. 

Accomplishments. The stagnation-point 
heating rates have been calculated at flight 
conditions along the current trajectory. The 
heating rates over the complete vehicle, 
forebody and afterbody, have been calculated 
at two flight conditions - time of peak heating 
and beginning of prime data period. The 
calculations are based on the best estimate of 
physical input data and also consider uncer- 
tainties in critical values such as catalytic 
wall coefficients and variations in radiation 
methods. The results have been supplied to 
key organizations within the AFE project. 
The results are being used in the heat shield 
design and in the design of the experiments. 
The stagnation-point heating rates are used 
to define correction factors to the heating 
method used in the heat shield design, and 
the heating rates over the complete configu- 
ration are used as anchor points in the design 
process. 


Significance. The heating rates, supplied 
by Langley, are based on the best techniques 
available, and are the most complete data set 
of accurate calculations. The AFE project has 
relied heavily on these results and our assess- 
ment of the heating environment in the ve- 
hicle design. Our results have shown that the 
heating environment for the current trajec- 
tory yield heat shield surface temperatures 
that are near the limit of the tile material. 
Any changes to the trajectory that result in 
increased heating rates could seriously jeop- 
ardize the survival of several experiments 
and acquisition of quality data. 

Status/Plans. The plan is to remain continu- 
ously involved in the calculation and assess- 
ment of the heating for the AFE. The impact 
in uncertainties of several critical parameters 
on the heating rates is being assessed. 

H. Harris Hamilton II 
Aerothermodynamics Branch 
Langley Research Center 
(804)864-4365 



REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of '".formation *s estimated to average 1 hour oer response, including the time for reviewing instructions, searching existing data sources, 
gathering and maintaining the data needed, and con-oietmg and reviewing the collection of information Send’comments regarding this burden estimate or any other aspect Of this 
collection c-f information, including suggestions tor reducing this ouraen to Wasmngton Headauarters Services. Directorate for information Operations and Reports, 12 IS Jefferson 
Davis Highway. Suite 12C4. Arlington. / L 222 02-4302. and tome Office o* Management and Budget. Paperwork Reduction Project (0704-0188). Washington, DC 20503. 


1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

April 1992 Technical Memorandum 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 

NASA Aerodynamics Program Annual Report 1991 NASW 4430 


6. AUTHOR(S) 

Louis J. Williams, Kristin A- Hessenius, Victor 
R. Corsiglia, Gary Hicks, Pamela F. Richardson, 
George Unger, Beniamin Neumann, and Jim Moss 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Office of Aeronautics and Space Technology 
Aerodynamics Division (OAST/RF) 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546 


10. SPONSORING /MONITORING 
AGENCY REPORT NUMBER 


NASA TM-4368 


11. SUPPLEMENTARY NOTES 

This report is an annual accomplishments review for the Aerodynamics 
Division (Code RF) .. The information contained herein covers FY 91. 


12a. DISTRIBUTION /AVAILABILITY STATEMENT 

Unclassified - Unlimited 
Subject Category 4^ 


12b. DISTRIBUTION CODE 


13. ABSTRACT (Maximum 200 words) 

This report is the annual accomplishments review for the Aerodynamics 
Division (Code RF) during FY 1991. The program includes both 
fundamental and applied research directed at the full spectrum of 
aerospace vehicles, from rotorcraft to planetary entry probes. This 
report contains a comprehensive review of the following aerodynamics 
elements: computational methods and applications, CFD validation, 

transition and turbulence physics, numerical aerodynamic simulation, 
test techniques and instrumentation, configuration aerodynamics, 
aeroacoustics , aerothermodynamics , hypersonics , subsonics , 
fighter/attack aircraft and rotorcraft. 



17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 
Unclassified 


19. SECURITY CLASSIFICATION 
OF ABSTRACT . . 

Unclassified 


15. NUMBER OF PAGES 
260 


16. PRICE CODE 
A12 


20. LIMITATION OF ABSTRACT 

Unlimited 


' MSN 7540-01-280-5500 


Standard Form 298 (Rev 2-89) 

prescribed by ANSI Std 739-18 
298-102 



















